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PROJECT  ABSTRACT 


The  objective  of  the  proposed  work  was  to  inform,  refine,  and  validate  a  recently  developed 
mechanism  based  predictive  model  of  Hydrogen  Environment  Cracking  (HEAC)  in  high  strength 
Alloys  with  a  focus  on  Monel  K-500.  In  this  proposal,  laboratory  based  material  hydrogen 
environment  assisted  cracking  (HEAC)  data  and  hydrogen  interaction  data  will  be  developed 
based  on  testing  of  material  obtained  from  field  returned  heats  of  aged  Monel  K-500  bolts  or 
studs  documented  to  have  HEAC  issues  compared  to  such  data  on  unfailed  or  commercially 
procured  material  all  satisfying  the  QQ-N-286  specification.  This  comparison  will  serve  to 
determine  whether  a  recently  developed  and  emerging  predictive  model  of  embrittlement  ( under 
Award  #:  N00014-10-1-0552  focusing  on  development  of  the  scientific  concepts  and  inputs  for 
the  stress  corrosion  life  prediction  program,  SCCrack,  developed  to  predict  crack  growth  rate 
and  life  time  to  failure )  for  aged  Monel  K-500  fasteners  subjected  to  cathodic  polarization  in 
marine  environments  must  be  modified  to  consider  strong  dependencies  of  the  heat-to-heat  or 
lot-to-lot  variations  in  material  when  attempting  to  assess  SCC  lifetimes  or  whether  the  model  is 
robust  in  this  regard  and  instead  should  focus  on  general  material,  environmental  and  geometric 
factors  such  as  applied  electrochemical  potential,  remote  stress,  geometry  and  surface  roughness. 
In  addition,  the  need  for  and  treatment  of  statistical  variations  in  input  parameters  will  become 
better  defined.  The  approach  would  be  to  test  failed  bolt  material  from  the  field  in  comparison  to 
unfailed  bolt  or  commercial  material  using  cutting  edge  world-class  crack  growth  and  hydrogen- 
metal  interaction  characterization  methods.  In  addition,  grain  boundary  microchemistry  would  be 
examined  via  Auger  or  other  surface  science  approaches  to  determine  whether  material 
microchemistry  and  intrinsic  susceptibility  vary  from  lot  to  lot.  This  work  was  accomplished  via 
a  research  team  led  by  Professors  J.  R.  Scully  and  James  T.  Bums,  students,  one  graduate 
students,  post-doctoral  researchers  as  well  as  UVA  characterization  facility  microscope  experts. 
The  payoff  from  this  project  was  improved  insight  into  on-going  questions  regarding  whether 
predictive  computational  models  of  component  SCC  lifetime  must  strongly  consider  heat  to  heat 
variations.  Ultimately  this  will  contribute  to  the  development  and  refinement  of  a  predictive  tool 
that  could  be  transitioned  to  DoD  labs  to  anticipate  and  manage  SCC  lifetimes  as  a  function  of 
many  user  controlled  inputs  such  as  cathodic  protection  design,  use  of  low  voltage  anodes,  bolt 
load,  bolt  roughness  etc.  by  providing  significant  input  regarding  whether  heat-to-heat  or  lot-to- 
lot  variability  in  material  must  be  strongly  considered  and  incorporated  into  such  models.  An 
unprecedented  opportunity  was  exploited  which  combined  field  information  with  cutting  edge 
computational  modeling  to  refine  models  of  stress  corrosion  cracking  optimize  model  outputs. 
The  work  described  herein  lead  to  two  peer  reviewed  publications,  numerous  talks,  and  the 
production  of  one  M.S.  in  Materials  Science  and  Engineering  with  subsequent  employment  in 
the  DOD  sector 
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ONR  TECHNICAL  BACKGROUND 


Introduction 

A  range  of  hydrogen  assisted  cracking  problems  currently  limits  use  of  high  strength  alloys  in 
the  marine  environment.  These  interactive  processes  are  highly  localized  at  the  crack  tip  and 
apply  to  Monel  k-500. 

The  alloy  of  focus  is  Monel  K-500  (63Ni-27Cu-3Al-0.5Ti).‘ Monel  K-500  (UNS  N05500)  is  a 
precipitation-hardenable  nickel-copper  alloy  that  combines  the  corrosion  resistance  of  MONEL 
Alloy  400  with  greater  strength  and  hardness.  It  is  strengthened  by  a  combination  of  cold  work 
and  aging  to  form  fine  Ni3(Ti,Al)  phases  The  alloy  is  used  for  pump  shafts,  oil-well  tools  and 
instruments,  surgical  blades  and  scrapers,  springs,  valve  trim,  fasteners,  and  marine  propeller 
shafts.  These  Ni-based  alloys  differ  considerably  from  quenched  and  tempered  steels  which 
derive  strength  from  a  body-centered  martensite  crystal  structure  and  secondary  age  hardening. 
A  range  of  properties  differ  from  H-metal  interactions  to  physical  properties.  However,  a 
common  feature  is  that  cracking  is  dominated  by  crack  tip  hydrogen  concentration  and  that 
HEAC  leads  to  intergranular  fracture  that  can  be  developed  in  cold  worked  and  aged  Monel  K- 
500  in  near  neutral  NaCl  solution  polarized  to  -1.0  V  vs.  SCE  and  subjected  to  slow-rising  crack 
mouth  opening  displacement. 

Seawater  HAC  susceptibility  of  Monel  K-500 

In  the  early  1960s,  Monel  Alloy  K-500  (UNS  N05500)  was  believed  to  be  relatively  immune  to 
hydrogen  environment  embrittlement  and  stress  corrosion  cracking".  Stress  corrosion  cracking  in 
H2S  environments  under  high  stress,  high  pressure  and  temperature  steam  and  HF  vapors  and  Hg 
were  noted  (for  Monel  alloy  400  or  Monel  Alloy  K-500  (UNS  N05500)).* 2 3 4 5 6  Reports  in  the  60s 
through  80s  began  to  indicate  susceptibility  to  room  temperature  hydrogen  environment 
embrittlement 4  5"  mainly  seen  as  a  low  in  tensile  ductility  and  notch  strength  ratio.  Scully  noted 
that  embrittlement  on  SHT  and  aged  Monel  Alloy  K-500  (UNS  N05500)  (Rockwell  C  25-35) 
occurred  in  seawater  at  room  temperature  when  coupled  to  zinc  and  that  hydrogen  levels 
increased  over  an  8  week  period  to  50  wppm.  A  strain  rate  effect  was  also  observed  even  in  pre¬ 
charged  material. 


1.  Special  Metals  product  data  sheet  on  Monel  Alloy  K-500  (UNS  N05500),  www.specialmetals.com 

2 .  R.  Copson  and  C.  F.  Cheng,  Corrosion,  vol.  12  (1956):  p.  647. 

3  K.  D.  Efird,  Materials  Performance,  vol.  24  (1985):  p.  37-40. 

4  J.  A.  Harris,  C.  D.  Stephens  and  Scarberr.Rc,  Corrosion,  vol.  28  (1972):  p.  57. 

5  C.  E.  Price  and  R.  S.  Fredell,  Metallurgical  Transactions  a-Physical  Metallurgy  and  Materials  Science,  vol.  17  (1986): 
p.  889-898. 

6J.  R.  Scully  and  M.  G.  Vassilaros,  DTNSRDC  SME  84-69,  (1984):  p.  1-28. 
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Since  these  initial  reports,  nickel-copper  alloy  bolt  failures  have  also  been  observed  numerous 
times  in  age  hardened  Monel  Alloy  K-500  (UNS  N05500)  subjected  to  normal  cathodic 
polarization  when  coupled  to  aluminum  anodes  in  seawater.  Failures  by  intergranular  cracking 
have  been  attributed  to  hydrogen  production,  uptake  and  embrittlement  as  a  result  of  cathodic 
polarization  and  was  first  attributed  to  high  thread  root  hardness  of  HRC  39  due  to  age  hardening 
after  threading.  It  was  recommended  that  first  annealing  and  then  age  hardening  be  conducted 
after  threading  to  limit  hardness  below  the  acceptable  hardness  (HRC  35)  of  Monel  Alloy  K-500 
(UNS  N05500).  This  was  a  suggested  hardness  limit  in  sour  systems.  This  hardness  limit  was 
also  suggested  by  a  Military  specification.  However,  additional  failures  have  occurred  in  roll- 
threaded  Monel  Alloy  K-500  (UNS  N05500)  bolts  annealed  at  980-1050°C,  water  quenched,  and 
precipitation  age  hardened  at  500-600°C  for  16  h  producing  a  Rockwell  C  hardness  of  only  25.7 8 9 
Embrittlement  failures  occurred  in  bolts  after  about  one  year  under  load  to  about  59%  of  the 
tensile  yield  strength  and  continuously  cathodically  protected  with  anode  grade  aluminum.  These 
failures  have  continued.  Intergranular  cracking  could  be  reproduced  under  slow  strain  rate  testing 
in  air  after  8  day  pre-charging  when  coupled  to  A1  or  charged  potentiostatically  (annealed  and 
aged  Monel  Alloy  K-500  (UNS  N05500))  and  to  a  lesser  extent  when  coupled  to  steel  or 
potentiostatically  polarized  to  -800  mV  vs.  Ag/AgCl  (only  seen  for  aged  Monel  Alloy  K-500 
(UNS  N05500))  in  ASTM  artificial  ocean  water. 

The  hydrogen  embrittlement  assisted  cracking  (HEAC)  behavior  of  other  Ni-based  superalloys  in 
gaseous  H2  and  aqueous  exposure  has  been  reviewed  in  detail  elsewhere.'  The  severity  of  this 
cracking  problem  is  apparent  in  many  nickel  base  alloys  such  as  Alloy  718."  However,  few 
studies  have  documented  the  threshold  stress  intensity  for  hydrogen  environment  embrittlement 
cracking  in  Monel  Alloy  K-500  (UNS  N05500)  and  there  were  no  reports  on  crack  growth  rate  in 
controlled  testing  prior  to  ONR  funded  work  under  the  prior  projects  N000 14-06- 1-0366  and 
N00014-10-1-0552  and  work  at  LUNA.  Raymond  noted  using  the  step  rising  load  method  that 
the  threshold  stress  intensity  of  aged  Monel  Alloy  K-500  (UNS  N05500)  (Rockwell  C  hardness 
=  30)  decreased  from  >69  MPaVm  at  -0.7  Vsce  to  <18  MPaVm  at  -1 .0  Vsce  in  room  temperature 
NaCl  solution/  Moreover,  little  was  known  about  hydrogen  uptake,  transport  and  trapping 
behavior  in  this  alloy  until  projects  N00014-06-1-0366  and  N00014-10-1-0552.  The  effective 
hydrogen  diffusion  coefficient  in  pure  single  crystal  Ni  was  8.4  x  10'1"  cm2/sec  with  an  activation 
energy  of  39.5  kJ/mole.  In  contrast,  diffusion  was  faster  in  copper  at  1.29  to  1.95  xl0'y  cm2/sec 
with  an  activation  energy  from  36-38  kJ/mole.  Harris  as  well  as  Scully  suggested  from  the 
limited  depth  of  intergranular  cracking  in  cylindrical  specimens  charged  with  hydrogen  that 
hydrogen  diffusion  rates  must  be  on  the  order  of  10'  cm  /sec.  However,  these  estimates  are 
indirect.  A  Ni-30  at.  %  Cu  alloy  had  a  hydrogen  diffusivity  of  3x1 0'1"  cm2/sec  and  an  effective 
activation  energy  for  diffusion  of  41  kJ/mole.v  Only  one  paper  reports  data  for  Monel  Alloy  K- 

7  L.  H.  Wolfe  and  M.  W.  Joosten,  SPE  Production  Engineering,  1988):  p.  382-386. 

8  L.  Raymond,  LRA  report  #  CTC"071024,(2008).  p.  1-16. 

9  D.  J.  Mitchell  and  E.  M.  Edge,  Journal  of  Applied  Physics,  vol.  57  (1985):  p.  5226-5235. 
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500  (UNS  N05500).  The  material  was  SHT  at  900°C  and  radiatively  cooled.  The  microstructure 
and  hardness  were  not  reported  but  minimal  aging  is  suggested  by  this  treatment.  Deff  was 
reported  to  be  7.6x10'"  cnT/sec.1"  In  recently  completed  work  under  our  previous  ONR  grant, 
Defrat  25  °C  was  suggested  to  be  1  to  3.3xlO'10  cm2/sec  for  aged  Monel  K-500.  However,  it  was 
felt  that  this  value  was  uncertain  by  a  factor  of  2  to  4.  This  difference  is  crucial  to  predictions  of 
stage  II  crack  growth  as  will  be  discussed  below  [da/dt]u. 

Moreover,  hydrogen  assisted  cracking  susceptibility  in  aqueous  environments  is  controlled  by 
and  increases  (to  a  theoretical  upper  bound)  with  the  crack  tip  diffusible  hydrogen  concentration 
which  is,  in  turn,  controlled  by  the  hydrogen  overpotential  and  production  rate.10  Metallurgical 
trap  states  are  important  and  either  serve  as  benign  sinks  for  hydrogen,  function  as  the  crack 
initiation  sites  and  the  crack  propagation  path  itself  or  function  as  a  reservoir  of  hydrogen  that 
can  be  supplied  to  the  tensile  triaxial  stress  field  of  the  stationary  or  moving  crack  tip. 


Additionally,  stage  II  environmental  crack  velocities  at  the  upper  bounds  of  hydrogen  uptake 
rates  where  cracking  is  not  surface  uptake  limited  is  often  controlled  by  the  effective  diffusion 
rate  of  hydrogen  at  the  crack  tip.  1(1  These  parameters  may  differ  between  SHT,  cold  worked  and 
aged  Monel  K-500  and  may  help  explain  observed  differences  in  susceptibility  of  aged  Monel  K- 
500  from  heat  to  heat.  For  all  these  reasons,  CH,difr,  the  diffusible  hydrogen  concentration 
absorbed  in  the  nickel  base  alloy  as  a  function  of  cathodic  potential,  and  the  effective  hydrogen 
diffusion  coefficient  Deff  are  critical  parameters  of  interest  in  order  to  predict  [da/dt]n  and  Kth  in 
UNS  N05500.  However,  Cn.diff  behavior  as  a  function  of  hydrogen  overpotential  or  cathodic 
hydrogen  production  rate  has  only  recently  been  established  in  Monel  Alloy  K-500  (UNS 
N05500)  in  previous  ONR  grants.  These  values  were  uncertain  by  a  factor  of  2-3.  Morever,  Defr 
and  trapping  behavior  have  just  been  explored  over  a  range  of  technologically  significant 
conditions  such  as  the  solution  heat  treated,  cold  worked  as  well  as  the  age  hardened  conditions. 
Refinement  of  these  parameters  is  crucial  to  not  only  the  accurate  prediction  of  Kth  and  [da/dt]n 
but  assessment  of  heat  to  heat  variations. 


Prediction  of  Stage  II  crack  velocity  [da/dfln  and  Kth  in  recent  work  under  ONR  N00014-06-1- 
0366  and  N00014-10-1-0552 


In  addition  to  the  threshold  stress  intensity  factor  for  HEAC,  is  it  important  to  understand  and 
model  predict  the  Stage  II  crack  growth  rate  as  mentioned  above.  This  information  will  be 
crucial  to  the  prediction  of  component  lifetimes.  Figure  1  shows  [da/dt]  vs  Kapp  for  aged  Monel 
K-500  (Allvac).  The  upper  bound  crack  growth  rate  is  ~2xl0'1  um/sec  or  2x1  O'4  mm/sec. 


10  B.  A.  Kehler  and  J.  R.  Scully,  Corrosion,  vol.  64  (2008):  p.  465-477. 
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A  substantial  correlation  of  the  fastest  rates  of  HEAC  were  obtained  from  past  experiments  at 
UVa,  using  a  common  experimental  method,  and  from  literature  results  for  high  strength  alloys 
stressed  during  exposure  in  chloride  solutions  HEAC  experiments  are  challenging  and  costly, 
and  should  be  augmented  by  modeling.  Diffusion-based  models  were  developed  to  describe 
da/dtn  for  HEAC  based  on  FPZ  supply  from  the  concentration  of  H  on  the  crack  tip  surface  and 
in  equilibrium  with  the  local  environment  (Cs),  as  reviewed  and  implemented  for  steel.  211  In 
general  models  of  this  type  are  of  varying  complexity  depending  on  whether  the  H-diffusion 
field  is:  (a)  concentration  and/or  stress  driven,  (b)  transient  or  steady  state,  (c)  ahead  of  a 
stationary  or  moving  crack,  (d)  modeled  in  one  or  two  dimensions  (e)  emanating  from 
environment-sensitive  Cs,  (f)  microstructure-trap  affected,  (g)  plastic  strain-trap  affected,  and  (h) 
coupled  with  a  H-failure  criterion.  The  results  of  diffusion  modeling  are  of  the  general  form: 
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Where:  Xcrit  is  the  critical  distance  for  H  diffusion  ahead  of  the  crack  tip,  Dneff  is  trap  modified 
H  diffusivity,  Co  is  the  concentration  of  H  and  governed  by  crack  electrochemistry  (equilibrium 
or  surface  reaction  rate  limited)  ays  is  alloy  tensile  strength,  t  is  time,  and  %  is  a  function  of  the 
indicated  variables  (equaling  between  4.0  and  0.01  depending  on  diffusion  boundary  condition). 
The  observed  and  predicted  proportionality  between  Stage  II  crack  growth  rate  and  Dnefr  in  Fig. 
2  supports  H-diffusion  control  of  HEAC.  An  upper  bound  xcrit  can  be  determined  from  these 
data.  The  highest-reasonable  value  of  ^  is  3  and  constant,  provided  that  the  da/dtn  are  not 
sensitive  to  Co  for  the  various  chloride  solutions  and  if  cracking  is  governed  by  constant 
Ccrit/Cs  of  order  2.  These  assumptions  result  in  Xcrit  =  0.7  pm  for  the  trend  line  in  Fig.  4, 
independent  of  alloy.  The  data  for  aged  Monel  K-500  are  consistent  with  the  trend  curve 
assuming  a  Deff  of~2xl0'1"  cm2/sec. 


uYongwon  Lee  and  Richard.  P.  Gangloff,  "Measurement  and  Modeling  of  Hydrogen  Environment  Assisted  Cracking  of  Ultra- 
High  Strength  Steel",  Metallurgical  and  Materials  Transactions.  A.  vol.  38,  pp.  2174-2190  (2007). 
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Figure  1.  Slow  rates  of  IG  SCC  in  Monel  K-500  (ST-A)  Figure  2.  Diffusion-based  upper  bound  level  of  da/dtn 
accurately  measured  in  6  day  experiment  with  slow-  for  HEAC  based  on  FPZ  supply  from  the  concentration 
rising  CMOD;  consistent  with  H-diffusion  limited  Stage  of  H  on  the  crack  tip  surface  and  in  equilibrium  with  the 
II  SCC.  (Allvac  Monel  K-500  tested  at  FTA  in  201 1).  local  environment  (Cs),  as  reviewed  and  implemented 

for  steel.  The  data  for  aged  Monel  K-500  are  consistent 
with  the  trend  curve  assuming  a  Dcn-of  2x1  O'10  cmVsec 

Application  of  Approach  to  Monel  K-500  under  cathodic  polarization  in  seawater 

Under  a  previously  funded  project,  ONR  N000 14-06- 1-0366,  [da/dt]n  and  KiH  were  measured 
and  modeled.  Examples  of  the  application  of  this  proposed  approach  to  Monel  K-500  are  shown 
in  Figure  2  where  hydrogen  uptake  law  was  explored  on  planar  electrodes  as  a  function  of 
hydrogen  overpotential  and  charging  current  density  in  an  alkaline  solution  designed  to  represent 
a  cathodically  polarized  crack  tip.  This  hydrogen  concentration  can  be  connected  to  crack 
geometry  through  Figure  6  where  hydrogen  overpotential  has  been  predicted  as  a  function  of  the 
crack  geometry  parameter  x2/G  and  applied  potential  at  the  crack  mouth.  These  data  were 
obtained  for  a  rescaled  crevice  and  based  on  cathodic  reaction  kinetics  measured  for  steel  and 
will  be  redone  for  Monel  K-500  in  the  proposed  work.  Crack  tip  hydrogen  concentration  must  be 
predicted  and  measurements  were  calibrated  against  LECO  measurements.  A  refined  assessment 
of  the  diffusible  hydrogen  content  as  a  function  of  applied  external  potential  and  crack  depth  was 
undertaken  for  Refinements  include  updated  determination  of  CH,difr  (probably  lowered)  which 
should  raise  predicted  values  of  Kjh  towards  experimentally  determined  values.  Additional 
crucial  pieces  of  information  required  as  input  to  da/dt  modeling  include  Dneff,  Xcrn  and  Co/,. 
D/zetf-was  determined  via  the  TDS  method  in  N00014-06-1-0366.  To  refine  estimates  of  Defr, 
these  measurements  are  corroborated  using  some  other  independent  method  in  the  current  work. 
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Figure  3.  CHdiff  and  CTOt  for  hot  finished  and  age- 
hardened,  Monel  K-500  ALLVAC  (63%Ni  -  33%Cu 
-  3.2%A1  -  2%Fe)  [UNS  N05500]. 


Figure  4.  CH,diff  versus  applied  potential  1  for  various 
x2/G  values.  The  data  applies  to  Allvac  aged  Monel  K- 
500  analyzed  after  300  sec  with  hydrogen  analysis  from 
TDS. 


Figure  5.  Predicted  C^difr  from  TDS  versus  applied  Figure  6.  Experimental  data  for  the  potential 
potential  for  SHT/aged  Monel  K-500  (ALLVAC)  at  dependence  of  [da/dt]n  for  hot  finished  and  aged 
various  X2/G  values.  Monel  K-500  (ALLVAC)  stressed  in  3.5%  NaCl 

solution.  Hydrogen  levels  were  obtained  from  TDS 
measurements. 

In  crack  growth  studies  under  work  under  ONR  N00014-06-1-0366  and  N000 14-1 0-1  -0552,  the 
following  findings  were  obtained.  The  crack  growth  rate  was  found  to  be  a  strong  function  of 


electrochemical  potential  (Figure  6)  and  therefore  crack  tip  hydrogen  concentration.  In  fact,  a 
substantial  increase  in  crack  growth  rate  is  predicted  as  a  function  of  Eapp  and  thus  crack  tip 
hydrogen  concentration  as  seen  in  Figures  7  and  8  and  seen  experimentally  in  aged  Monel  K-500 
(Figure  8).  This  means  that: 

•  HEAC  of  aged  Monel  K-500  will  be  a  strong  function  of  applied  potential  and  will  likely 
be  extremely  sensitive  to  applied  cathodic  potential,  local  potential,  geometry,  materials 
and  stress  state. 

•  Any  metallurgical  variable  or  heat-to-heat  variable  not  captured  by  gross  changes  in 
material  hardness  that  might  change  the  trajectory  of  [da/dt]II  with  Eapp  as  seen  in  Figures 


6-8  could  significant  change  the 
embrittlement. 


Figure  7.  Prediction  of  crack  growth  rate  in  SHT-aged 
K500  as  a  function  of  reduced  hydrogen  concentration. 
Ys  =  783  M  Pa  aged  Monel  K-500.  Vh  -  1.73  cmVmole 
H  from  Moody,  Baskes  for  Nickel  see  Ha  spreadsheet. 
ays=  783  MPa  aH  =  6  ays  VH  =  1.73  cmVmole  II 


lifetime  of  components  exposed  to  hydrogen 


Figure  8.  Crack  Growth  rate  in  SHT-aged  K-500: 
Experiment  vs.  Calculation  Using  TDS  data(ALLVAC) 
to  calculate  CH,difrat  Eappiied  at  x2/G  -  10  cm.  The  first 
experimental  data  point  ALLVAC)  is  at  -700m V  vs. 
SCE.  The  last  is  at  -1100  mV  vs.  SCE.  These  data 
indicate  the  effect  of  potential. 
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SUMMARY  BASED  ON  FINDINGS  UNDER  ONR  N00014-06-1-0366  and  N00014-10-1- 

0552 


The  crack  propagation  perspective  provides  a  framework  to  quantify  subcritical  cracking 
phenomena  for  component  life  assessment  in  a  marine  application.  Improvements  to  H  cracking 
models  were  achieved  in  prior  work  by  developing  new  understanding  of  crack  tip 
electrochemistry  and  mechanics,  as  guided  by  the  existing  interdisciplinary  experience  base.  The 
Monel  K-500  alloys  system  provides  an  ideal  case,  of  both  practical  and  scientific  relevance,  to 
extend  understanding  of  HEAC  modeling.  By  coupling  with  external  definition  of  the  mechanics 
condition,  the  lifetime  of  the  component  can  be  estimated  as  long  as  the  state  of  stress,  material, 
geometry,  electrochemical  conditions,  and  temperature  are  defined.  However,  this  work  fills  a 
crucial  gap.  Specifically,  assessing  how  heat  to  heat  variations  affect  model  inputs  and  how 

statistical  variations  should  be  treated  w  ith  a  rational  technical  basis  inside  of  the  SCCcrack 

model. 


OBJECTIVES 


The  overall  objective  of  this  interdisciplinary  research  program  is  to  implement  quantitative- 
predictive  models  of  hydrogen  environment  assisted  cracking  (HEAC)  of  Monel  K-500  in  a 
pseudo-application  specific  situation. 

The  objective  of  this  work  was  to  inform,  refine,  and  validate  recently  developed  mechanism 
based  predictive  model  of  Hydrogen  Environment  Cracking  (HEAC)  in  High  Strength  Alloys 
with  a  focus  on  Monel  K-500.  Lab  based  material  cracking  data  and  hydrogen  interaction  data 
were  developed  based  on  testing  of  material  obtained  from  failed  heats  of  Monel  K-500  bolts  or 
studs  as  compared  to  such  data  on  unfailed  material  all  satisfying  the  QQ-N-286  specification. 
This  comparison  served  to  determine  whether  the  UVA  recently  developed  and  emerging 
predictive  model  of  embrittlement  must  be  modified  to  consider  strong  dependencies  of  the  heat- 
to-heat  variations  in  material  on  SCC  lifetimes  or  whether  the  model  is  robust  in  this  regard  and 
instead  should  focus  on  general  material,  environmental  and  geometric  factors  such  as  applied 
electrochemical  potential,  remote  stress,  geometry  and  surface  roughness. 

The  specific  objectives  of  this  work  were  to  ascertain  from  field  failed  as  well  as  commercial 
purchased  Monel  K-500: 


Quantitative  characterization  of  the  HEAC  susceptibility  of  several  heats  of  aged 
Monel  K-500  over  a  range  of  cathodic  potentials  relevant  to  cathodic  protection  in 
seawater  via  state-of-the-art  crack  growth  assessment  techniques. 


10 


•  Characterization  of  basic  metallurgical  properties  of  several  heats  of  aged  Monel 
K-500  material  as  well  as  the  grain  boundary  microchemistry  that  may  control 
variations  in  intergranular  HEAC  susceptibility. 

•  Characterization  of  the  number  of  traps  and  the  trap  binding  energies  and  other 
hydrogen  related  properties  of  several  heats  of  aged  Monel  k-500  material. 

•  Determination  of  the  heat  to  heat  and  within  heat  statistical  variations  in  hydrogen 
concentrations  and  [da/dt]n  rates  over  a  range  of  pertinent  cathodic  potentials. 


TECHNICAL  APPROACH 


The  objectives  were  met  by  study  of  six  separate  lots  of  aged  Monel  K-500  each  of  which  meet 
the  QQ-N-286  specification  that  governs  procurement  of  this  material.  One  of  the  material  heats 
was  obtained  from  a  commercial  supplier  (Allvac)  and  5  separate  material  lots  were  supplied  by 
NSWC-CD  from  field  returned  aged  Monel  K-500  components.  The  research  effort  was  divided 
into  three  main  thrusts. 

First  a  detailed  characterization  of  the  metallurgy  and  baseline  material  properties  of 
each  lot  was  performed  to  identity  the  microstructural  differences  between  each  lot  of 
material. 

Second,  the  HEAC  susceptibility  of  each  lot  of  material  was  established  for  various 
environment  via  both  slow-strain  rate  tensile  testing  and  state-of-the-art  fracture 
mechanics  based  characterization  of  the  crack  growth  behavior  as  a  function  of  the 
applied  stress  intensity  (K). 

Third,  the  H-metal  interactions  were  assessed  for  each  material  lot.  The  data  from 
each  of  these  efforts  were  interpreted  in  the  context  of  previously  developed  and 
validated  micro-mechanical  models  to  quantitatively  evaluate  how  the  observed 
changes  in  microstructure  would  result  in  the  measured  changes  in  HEAC 
susceptibility. 


Detailed  motivation,  experimental  methods,  results,  discussion  and  conclusions  for  the  efforts  in 
each  of  the  three  thrust  areas  are  reported  in  the  following  two  chapters.  Figures  are  numbered  in 
each  section  chronologically. 

Chapter  1  details  the  materials  used,  the  microstructural  characterization  relevant  to 
the  HEAC  behavior,  the  LEFM  fracture  process  and  results,  and  the  micro¬ 
mechanical  modeling.  . 

Chapter  2  details  the  H-metal  interaction  analysis,  slow-strain  rate  results,  Auger 
analysis,  and  further  micro-mechanical  modeling. 
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SUMMARY  OF  FINDINGS 


The  detailed  findings  of  this  work  are  incorporated  in  each  chapter.  However,  three  at-large 
findings  traced  to  the  above  primary  objectives  were  met.  Approaches  for  incorporating  heat  to 
heat  variations  in  fracture  modeling  from  either  the  fracture  mechanics  perspective  and  notched 
tensile  bar  perspectives  were  proposed.  Possible  additions  to  the  federal  specification  based  on 
beneficial  and  detrimental  elements  and  a  possible  equation  for  trace  element  composition  were 
proposed.  Lastly  a  few  test  methods  were  proposed  that  can  be  used  to  elucidate  relative 
susceptibility.  Specifically,  the  slow  strain  rate  method  may  be  utilized  to  rapidly  assess  alloys 
prone  to  intergranular  hydrogen  cracking  at  intermediate  cathodic  protection  levels. 


TRANSITIONS 


The  long  range  payoff  of  this  research  for  the  USN  will  be  the  ability  to  predict  the 
environmental  cracking  threshold  stress  intensity  and  crack  growth  rate  for  any  nickel  based 
alloy  such  as  Monel  K-500  fasteners  in  a  marine  environment,  as  a  function  of  applied 
electrochemical  potential,  defect  geometry,  alloy  metallurgy  and  state  of  stress  with  an 
assessment  of  heat  to  heat  variations.  The  question  of  whether  and  how  predictive  computational 
models  of  component  SCC  lifetime  can  consider  heat  to  heat  variations  was  addressed. 
Ultimately  this  work  will  enhance  the  development,  refinement  and  the  use  of  a  predictive  tool 
that  could  be  transitioned  to  Navy  labs  to  anticipate  and  manage  SCC  lifetimes  as  a  function  of 
cathodic  protection  design,  low  voltage  anodes,  bolt  loads,  bolt  roughness  etc.  with  significant 
input  on  whether  heat-to-heat  or  lot-to-lot  variability  in  material  must  be  strongly  considered  and 
incorporated  into  such  models.  An  unprecedented  opportunity  combine  field  information  with 
cutting  edge  computational  modeling  was  undertaken. 

For  instance  these  alloys  are  used  in  flooded  spaces  and  experience  various  degrees  of  cathodic 
protection.  Such  predictions  will  be  sufficiently  quantitative  for  direct  use  in  component-level 
analysis  including  time  to  failure.  The  output  will  be  suitable  for  a  range  of  components  but 
including  high  strength  fasteners,  rods  and  studs  in  wetted  marine  environments.  The  models 
produced  to  achieve  this  capability  will  ultimately  be  sufficiently  simple  to  be  implemented  on  a 
desktop  computer  in  the  engineering  environment. 


Insight  might  also  be  provided  leading  to  the  further  refinement  of  QQ-N-286  standard  leading  to 
an  enhanced  ability  to  weed  out  HEAC  susceptible  tempers  of  Monel  K-500  when  subjected  to 
cathodic  protection  in  seawater. 
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CHAPTER  1 


The  Effect  of  Microstructure  Variation  on  the  Hydrogen  Environment-Assisted  Cracking  of  Monel 

K-500 

Zachary  D.  Harris,  Justin  D.  Dolph,  Greger  L.  Pioszak,  Brendy  C.  Rincon  Troconis,  John  R.  Scully,  and 

James  T.  Bums 

Center  for  Electrochemical  Science  and  Engineering 
Department  of  Materials  Science  and  Engineering 
University  of  Virginia 
Charlottesville,  VA  22904 


Abstract 

The  influence  of  microstructure  variation  on  hydrogen  environment-assisted  cracking  (HEAC)  of 
Monel  K-500  was  evaluated  using  five  nominally  peak-aged  lots  of  material  tested  under  slow-rising 
stress  intensity  loading  while  immersed  in  NaCl  solution  under  cathodic  polarizations.  Minimal  variation 
in  HEAC  resistance  amongst  material  lots  was  observed  for  an  applied  potential  of  -950  mVScE  (Eapp,  vs. 
saturated  calomel),  while  significant  differences  in  both  fracture  morphology  and  HEAC  resistance  were 
observed  at  the  less  negative  potential  of  -850  mVScE>  suggesting  that  relatively  severe  H  environments 
produce  sufficient  crack-tip  H  to  minimize  the  impact  of  metallurgical  differences.  Sensitivity  analyses 
accomplished  by  varying  the  inputs  used  in  decohesion-based,  micromechanical  models  imply  significant 
variations  in  HEAC  resistance  are  possible  for  realistic  changes  in  grain  boundary  toughness,  hydrogen 
uptake  behavior,  and  yield  strength.  Grain  size,  impurity  segregation  (including  the  effects  of  gettering 
elements),  grain  boundary  character/connectivity,  and  crack  path  tortuosity  are  also  considered  in  the 
context  of  HEAC  susceptibility.  Yield  strength,  global  hydrogen  content,  as  well  as  impurity  segregation 
to  grain  boundaries,  especially  boron  and  sulfur,  are  speculatively  considered  to  be  the  dominant 
contributions  in  determining  HEAC  resistance.  Modifications  that  would  incorporate  the  effects  of  grain 
boundary  segregation  are  proposed  for  the  KTH  model;  detailed  validation  of  such  changes  require  high 
fidelity  and  quantitative  inputs  for  the  degree  of  grain  boundary  segregation.  Regardless,  fracture 
mechanics-based  HEAC  results,  detailed  microstructure  characterization,  and  micromechanical  modeling 
were  successfully  coupled  to  gain  insights  into  the  influences  governing  the  microstructure-dependent 
HEAC  susceptibility  of  Monel  K-500. 
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Introduction 


Monel  K-500  is  a  Ni-Cu-Al  superalloy  commonly  used  in  marine  applications  due  to  its  excellent 
combination  of  corrosion  resistance,  strength,  and  fracture  toughness.  For  such  applications, 
environment-assisted  cracking  (EAC)  is  historically  evaluated  via  a  wide-  and  diverse-range  of  testing 
approaches  (e.g.  NACE  MR0175  and  ASTM  G30,  G39,  G129,  E1681,  and  F1624  [1-6]).  Material 
selection  engineers  will  often  use  data  from  such  tests  in  a  “go/no-go”  paradigm  where  the  long-term 
EAC  resistance  is  inferred  by  empirical  material  testing.  Specifically,  if  testing  indicates  that  the  service 
environment  and  loading  conditions  are  below  a  threshold  level  of  severity,  it  is  then  assumed  that 
cracking  will  not  occur  over  the  lifetime  of  the  component.  Such  an  approach  is  codified  in  NACE 
MR0175  which  outlines  guidelines  for  the  selection  of  “cracking  resistant  materials’'  in  H2S  environments 
(of  note  is  that  this  specification  does  state  that  materials  meeting  this  standard  are  “not  necessarily” 
immune  under  all  service  conditions).  The  inadequacy  of  this  “go/no-go”  approach  is  demonstrated  by 
isolated  long-time  (~10  year)  service  failures  of  Monel  K-500  components  on  North  Sea  oil  and  gas 
platforms,  deep  gas  wells,  and  US  Navy  sea  vessels  [7-13].  In  each  of  these  cases,  the  Monel  K-500 
components  were  immersed  in  seawater  and  cathodically  polarized  to  protect  adjacent  carbon  and  low 
alloy  steel  structures.  This  service  environment,  combined  with  the  observation  of  fracture  surfaces  that 
exhibited  intergranular  (IG)  with  small  amounts  of  transgranular  (TG)  crack  growth  [7,8,11-13], 
suggests  that  the  failure  of  these  components  can  be  attributed  to  EAC  -  an  assessment  corroborated  by 
laboratory  testing  of  Ni-based  superalloys  [14-18]. 

These  service  failures  expose  critical,  non-conservative  flaws  in  the  existing  approach  to  EAC 
management  for  current-  and  next-generation  materials.  Specifically,  current  EAC  characterization 
techniques,  material  selection  criteria,  and  structural  integrity  prognosis  approaches  fail  to  fully  account 
for  conditions  where  slow  sub-critical  crack  propagation  occurs  over  a  significant  portion  of  component 
life.  This  deficiency  demonstrates  the  need  for  high-fidelity  linear  elastic  fracture  mechanics-based 
(LEFM)  material  characterization  of  both  the  threshold  stress  intensity  for  EAC  (KTh)  and  the  steady  state 
(stage  II)  sub-critical  crack  growth  kinetics  (da/dtu)  for  the  pertinent  service  environment  [19]. 
Additionally,  if  slow  crack  propagation  dominates  the  total  service  life  of  the  component,  then  sub-critical 
crack  growth  kinetics  can  be  coupled  with  LEFM-based  life  prediction  methods  (e.g.  SCCrack ,  recently 
developed  for  EAC  modeling  [20])  to  provide  a  powerful  tool  for  structural  design, 
material/coating/cathodic  protection  selection,  and/or  structural  life  management.  An  analogous  approach 
is  well  established  in  the  fatigue  community  where  LEFM-based  programs  are  commonly  used  to  aid 
structural  life  management  [21-23].  Ultimately,  the  rigor  of  the  LEFM  approach  is  critically  dependent 
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on  the  applicability  of  the  crack  growth  kinetics  data  to  the  relevant  environment,  loading  conditions,  and 
material  of  interest.  An  initial  study  evaluated  the  influence  of  the  former  two  variables  on  EAC  behavior 
in  Monel  K-500  [18],  but  the  effect  of  lot-to-lot  metallurgical  variation  on  HEAC  susceptibility  is  not 
known. 


Crack  growth  kinetics  for  environment/cathodic  polarization  combinations  pertinent  to  marine 
service  conditions  were  established  in  prior  work  using  direct  current  potential  difference  (dcPD)  [24]  to 
monitor  crack  progression  during  slow-rising  displacement  testing  of  single-edge  notched  (SEN) 
specimens  from  a  single  lot  of  Monel  K-500  [18].  The  resulting  data  quantitatively  demonstrated  that 
decreasing  the  applied  potential  from  -750  mVScE  to  -1 100  mVScE  in  0.6  M  NaCl  solution  will  decrease 
the  EAC  threshold  (Kth)  and  increase  the  stage  II  crack  growth  rate  (da/dtn).  In  addition  to  their 
engineering  relevance,  these  data  provided  the  basis  for  micromechanical  modeling  of  potential- 
dependent  EAC  in  the  context  of  a  decohesion-based  hydrogen  (H)  embrittlement  mechanism.  Estimates 
of  the  overpotential  for  H  production  at  the  crack-tip  were  established  as  a  function  of  bulk  applied 
potential  via  rescaled  crevice  experiments  [18].  These  data  were  coupled  with  effective  H  solubility 
measurements  versus  H  overpotential  to  establish  that  the  crack-tip  diffusible  H  concentration  (CH-Diff) 
rises  with  decreasing  applied  potential  from  -750  to  -1100  mVScE  [25].  As  the  level  of  cathodic 
polarization  increases,  the  increase  in  EAC  susceptibility  directly  scales  with  increasing  crack-tip 
diffusible  H  concentration,  thereby  suggesting  that  H  environment-assisted  cracking  (HEAC)  is  the 
governing  mechanism  for  sub-critical  cracking  in  Monel  K-500.  This  finding  is  consistent  with  a  wide 
literature  that  supports  the  H-embrittlement  phenomenon  in  Ni-based  alloys  [26-31] 


While  several  viable  mechanisms  for  hydrogen  embrittlement  in  Ni-based  alloys  exist  [32-34], 
hydrogen-enhanced  decohesion-based  (HEDE)  models  provide  a  means  to  quantitatively  relate  the  stress- 
enhanced  H  concentration  (CH-CT)  in  the  fracture  process  zone  (FPZ)  to  the  relevant  HEAC  metrics  (KTH 
and  da/dtn)  [35-38].  Critically,  CH.a  can  be  calculated  from  the  empirically  determined  Cn-Diff,  thereby 
establishing  CH.aas  the  key  environmental  variable  in  the  HEDE  models  for  KTH  and  da/dtn.  Prior  work 
argued  that  this  relationship  (Eqn  1)  can  be  described  with  Fermi-Dirac  statistics  by  reasonably  assuming 
a  single  dominant  reversible  trap  site  and  low  lattice  H  concentrations  [18]. 


Ch-o~  CH.Diff  exp 


PhVh' 
.  RT  . 


(Eqn  1) 


where  VH  is  the  partial  molar  volume  of  hydrogen  in  a  nickel  lattice  (1.73  cm3/mol  H)  [18],  R  is  the  gas 
constant,  T  is  temperature,  and  aH  is  the  hydrostatic  stress  located  a  distance  of  xcrit  ahead  of  the  crack- 
tip.  Depending  on  the  approach  (discrete  dislocation  models,  blunted  crack  finite  element  models,  or 
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strain  gradient  plasticity  models),  crack-tip  mechanics  estimates  of  aH  vary  from  3  to  20  times  the  yield 
strength  of  the  alloy  [39]. 


Kth  for  HEAC  is  predicted  using  the  decohesion-based,  micromechanical  model  proposed  by 
Gerberich  et  al.  [32,35,38]  (Eqn  2)  to  describe  the  near  crack-tip  stress  state: 


Km  =^exP 
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(Eqn  2) 


where  P’  ((MPaVm)'1)  and  a”  (MPa-m)  are  constants  from  modeling  the  dislocation  shielding  of  the 
crack-tip  elastic  stress,  k]G  (MPaVm)  is  the  intrinsic  Griffith  toughness  of  the  grain  boundary,  a 
(MPaVm(atom  fraction  of  hydrogen)’1)  is  a  weighting  factor  which  defines  the  decrease  in  kiG  per  unit  H 
concentration,  and  oYs  (MPa)  is  the  yield  strength  [18,32,35].  The  steady-state  (stage  II)  sub-critical 
HEAC  growth  behavior  is  limited  by  the  diffusion  of  H  from  the  crack-tip  surface  to  the  location  of 
maximum  stress  within  the  FPZ  [18,37].  As  such,  a  HEDE-based  model  for  H-diffusion  ahead  of  a 
discontinuously  propagating  crack  results  in  the  relationship  given  in  Eqn  3. 
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where  Deff  is  the  effective  diffusivity  of  hydrogen  which  accounts  for  the  presence  of  trap  sites  in  the  alloy 
(reported  values  are  bounded  by  5.0xl0'n  cm2/s  and  1 .3xl0"lu  cm2/s  for  Monel  K-500  [25]),  x^t  is  the 
critical  distance  ahead  of  the  crack-tip  where  discontinuous  H-assisted  crack  advance  nucleates,  and  CHa- 
cnt  is  the  critical  concentration  of  hydrogen  necessary  to  activate  decohesion  at  xcrit  [18].  The  parameter 
Xcnt  is  on  the  order  of  1  pm  for  many  high  strength  alloys  and  is  related  to  the  position  of  maximum 
hydrostatic  stress  in  front  of  the  crack  tip  [40].  Prior  work  demonstrated  that  using  a  single 
experimentally  calibrated  parameter  and  the  empirically  derived  applied  potential  versus  CH-Difr 
relationship  gave  predictions  from  Eqn  2  and  3  that  showed  excellent  agreement  with  measured  HEAC 
parameters  for  applied  potentials  ranging  from  -700  to  -1 100  mVsce  in  a  single  lot  of  Monel  K-500  [18]. 
The  success  of  these  HEDE-based  micromechanical  modeling  paradigms,  along  with  the  high  fidelity  of 
the  dcPD  technique  for  measuring  crack  growth,  justify  the  extension  of  these  approaches  to  understand 
the  effect  of  metallurgical  variation  on  HEAC  behavior. 

The  effect  of  lot-to-lot  variations  in  metallurgy  on  HEAC  susceptibility  is  especially  relevant  for 

Monel  K-500.  The  specification  governing  material  procurement  (Federal  Specification  QQ-N-286G)  is  a 
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property-driven  document  that  specifies  benchmarks  for  composition,  mechanical  properties,  IG 
susceptibility,  and  grain  size  uniformity  [41].  However,  while  there  are  suggested  heat  treatments,  the 
specification  allows  for  age  hardening  to  “be  accomplished  by  any  process  proved  adequate  to  meet  the 
mechanical  property  requirements"  [41].  This  flexibility  in  production  can  lead  to  vast  variations  in 
strength  and  metallurgical  features  and  their  effects  are  unlikely  to  be  captured  by  the  rudimentary  testing 
mandated  by  QQ-N-286G  for  determining  IG  failure  susceptibility.  Changes  in  the  grain  size,  the 
distribution  of  grain  boundary  misorientation  angles,  the  bulk  impurity  concentration,  the  degree  of 
impurity  segregation,  the  precipitate  morphology,  and  the  post-processing  dislocation  density  will  vary 
with  material  lot/processing  route.  Critically,  these  changes  can  influence  the  HEAC  behavior  by 
modifying  the  material  strength,  slip  behavior,  grain  boundary  cohesion  characteristics,  crack  path 
connectivity/tortuosity,  and/or  H-metal  interactions.  The  micromechanical  models  detailed  in  Eqn.  2  and 
3  provide  a  quantitative  framework  to  link  these  metallurgy  variations  with  the  HEAC  behavior  by 
speculating  how  lot-to-lot  differences  vary  the  pertinent  input  parameters  (kIG,  a,  aY s,  x^t,  CH,Diff3  CH.a3 
Dcff,  and  CHa-crit)*  However,  a  systematic  characterization  and  analysis  of  the  relative  influence  of  changes 
in  these  metallurgical  features  on  the  HEAC  behavior  is  lacking  for  Monel  K-500. 

The  objective  of  this  paper  is  to  evaluate  the  extent  to  which  the  HEAC  behavior  varies  for  Monel 
K-500  lots  that  meet  the  QQ-N-286G  specification.  First,  HEAC  susceptibility  is  assessed  through 
fracture  mechanics  testing  in  both  an  inert  environment  and  in  full  immersion  NaCl  solution  under  mild 
and  aggressive  cathodic  polarizations.  Second,  a  detailed  characterization  of  the  microstructure  is 
performed  to  identify  the  variation  in  metallurgical  features  for  each  tested  lot  of  Monel  K-500.  Third, 
several  aspects  of  the  H-metal  interactions  are  investigated.  Finally,  the  results  of  these  analyses  are 
coupled  with  the  micromechanical  model  for  KTH  (Eqns  1  and  2)  to  quantitatively  evaluate  how  observed 
changes  in  the  microstructure  may  relate  to  the  resulting  HEAC  behavior. 

Experimental  Methods 

2.1  Material 

Monel  K-500  is  composed  of  an  austenitic  Ni-Cu  solid  solution  matrix  (y)  that  is  precipitation- 
hardened  by  the  intermetallic  y’  (Ni3(Al,Ti)  phase.  The  y’  precipitates  are  homogeneously  distributed  and 
highly  coherent  (<0.1%  misfit  strain)  with  an  ordered  Ll2  structure  composed  ofNi  atoms  at  the  faces  and 
A1  (or  Ti)  atoms  on  the  corners  of  the  unit  cell  [42-44].  Due  to  this  low  misfit,  these  precipitates  form  as 
spherical  particles  -  a  morphology  that  is  maintained  even  after  prolonged  aging  [43].  Dislocation 
interactions  with  y5  are  generally  controlled  by  the  size  of  the  y5  precipitates,  with  a  noted  shift  from 
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dislocation  cutting  to  Orowan  looping  for  larger  particle  radii  [43,44].  A  heterogeneous  distribution  of 
MC-type  carbides  (typically  TiC)  are  also  found  in  this  alloy,  but  this  distribution  does  not  change  with 
aging  and  their  contribution  to  the  strength  of  Monel  K-500  is  considered  to  be  minimal  [43]. 

Five  lots  of  Monel  K-500  were  evaluated  in  this  work.  The  first  lot  was  obtained  from  ATI  Allvac 
(termed  Allvac)  as  a  15.9  mm  diameter  bar  that  was  hot  worked  and  then  age  hardened  as  follows:  16  h  at 
866K,  cooled  at  14K/h  to  755K,  and  then  air  cooled  to  ambient  temperature.  Three  lots  of  Monel  K-500 
were  harvested  from  retired  engineering  components  (intact  components  that  were  roughly  305  x  51  x  38 
mm  in  size).  The  exact  service  conditions  for  the  components  are  not  known,  but  are  estimated  as 
roughly  10-15  years  of  full  immersion  in  seawater  with  exposure  to  applied  potentials  likely  ranging  from 
-1000  to  -750  mVscc.  These  lots  are  termed  Component  1  (TR  1),  Component  2  (TR  2),  and  Component 
3  (TR  3)  and  have  an  unknown  processing  history.12  The  final  lot  of  Monel  K-500  was  received  from  the 
U.S.  Naval  Research  Laboratory  (NRL)  with  thermal  processing  designed  to  capture  the  upper-bound 
strength  requirement  of  QQ-N-286.  This  material  (termed  NRL-High  Strength  (NRL-HS)),  was 
processed  from  a  1 14  mm  diameter  bar  procured  from  Special  Metals  that  was  direct-aged  at  866  K  for  16 
h,  furnace  cooled  at  14  K/h  to  81 1  K,  held  for  1  h,  furnace  cooled  at  14  K/h  to  755  K,  held  for  1  h,  and  air 
cooled  to  ambient  temperature.  Table  1  presents  both  the  QQ-N-286  composition  requirements  and  the 
measured  bulk  chemical  compositions  for  each  lot;  concentrations  of  trace  elements  S,  P,  Sn,  Zr,  Mg,  Hf, 
and  Pb  were  determined  using  glow  discharge  mass  spectrometry.  Measured  mechanical  properties  are 
reported  in  Table  2;  the  yield  strength  and  Ramberg-Osgood  flow  rule  constants  were  determined  from 
compression  loading  of  cylindrical  specimens  to  2%  total  strain. 

2.2  Fracture  Mechanics-Based  HEAC  Testing 

2.2,1  Mechanical  Loading  and  Crack  Monitoring 

Single  edge  notch  tensile  (SEN(T))  specimens  were  machined  from  each  of  the  material  lots  with 
gauge  sections  of  12.5  mm  (10,6  mm  for  Allvac)  in  width  (W),  and  2.68  mm  in  thickness  (B).  Specimens 
were  loaded  parallel  to  the  bar/component  longitudinal  axis  with  the  Mode  I  crack  growth  occurring  in  the 
radial  direction.  An  initial  200  ±  10  pm  notch  was  placed  in  the  center  of  the  gauge  section  using 


12  HEAC  testing  of  the  component  specimens  were  performed  in  the  as-received 
conditions.  A  companion  study  suggests  that  H-charging  that  occurred  during 
field  service  may  lead  to  a  total  internal  H-concentration,  CH_totai,  on  the 
order  of  5  wppm  which  is  similar  to  the  6  wppm  at  the  crack  tip  at  -850  mVScE 
for  the  Allvac  material  lot  [18] .  The  effects  of  the  precharged  hydrogen  is 
discussed  in  more  detail  elsewhere  [124].  However,  HEAC  testing  of  the 
component  specimens  in  an  inert  environment  did  not  show  H-enhanced  IG 
cracking  due  to  the  precharged  H-content. 
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electrical  discharge  machining  (EDM).  Specimens  were  fatigue  precracked  in  laboratory  air  using  the 
following  protocol:  (1)  constant  maximum  load  (12.5  kN)  at  R  =  0.1  from  the  initial  notch  to  a  crack 
length  of  ~0.4  mm  (ending  Kmax~15  MPaVm)  and  (2)  constant  R  =  0.1  with  Kmax  decreasing  from  15 
MPaVm  to  a  final  value  of  13  MPaVm  at  a  crack  length  of  approximately  1  mm.  All  specimens  were 
loaded  in  clevis  grips  to  allow  for  free  rotation  in  compliance  with  K-solution  boundary  conditions  [45]. 

Active  crack  length  feedback  coupled  with  software-controlled  servo-hydraulic  actuator 
displacement  enabled  fracture  mechanics  testing  at  a  constant  increasing  elastic  K  (dK/dt).  Testing  in 
NaCl  solution  was  conducted  at  a  constant  dK/dt  =  0.33  MPaVm/hr,  which  corresponds  to  an  approximate 
grip  velocity  of  2xl0'6  mm/s,  while  inert  environment  testing  was  completed  at  a  constant  dK/dt  of  either 
0.33  or  1  MPaVm/hr  to  limit  test  duration.  Crack  progression  was  monitored  via  dcPD  and  the  grip 
velocity  was  adjusted  every  5  seconds  to  maintain  the  desired  dK/dt.  This  procedure  included  current- 
polarity'  reversal  to  eliminate  thermally-induced  voltages,  crack  potential  difference  normalization  by  a 
remote  reference,  and  active-normalized  potential  values  averaged  in  groups  of  50.  Each  potential 
measurement  was  converted  to  a  crack  length  via  Johnson’s  equation  [46],  where  the  potential  associated 
with  the  initial  notch  depth  (ao)  is  the  adjustable  constant  (V0)  [47].  The  resolvable  average  crack 
extension,  Aa,  for  this  technique  is  0.5  pm,  based  on  the  ability  to  resolve  0.1  pV  changes  in  measured 
potential  for  a  constant  applied  current  of  4.000  ^0.005  A.  Post-test  initial  crack  length  measured  via 
electron  microscopy  typically  varied  <5%  from  the  dcPD-based  calculated  value.  If  necessary,  the  initial 
crack  length  was  corrected  to  match  the  measured  value  from  fractography  and  all  subsequent  crack 
measurements  were  recalculated  using  the  measured  voltage  at  each  time  step.  The  crack  growth  rate 
(da/dt)  at  each  crack  length  and  time  combination  was  determined  using  the  incremental  (n=3) 
polynomial  method  described  in  ASTM  E647-13  Appendix  XI  [48].  High-resolution  crack  growth 
monitoring  for  environment-material  combinations  that  are  resistant  to  HEAC  is  complicated  by 
plasticity-based  resistivity  increases,  which  establish  a  functional  resolution  limit  of  0.2  to  1  nm/s, 
depending  on  the  applied  K  and  dK/dt  [18].  False  da/dt  were  observed  at  low  K  due  to  crack  surface 
electrical  contact  enabled  by  the  cathodic  polarization-induced  destabilization  of  the  surface  oxide  film;  a 
post -test  analysis  protocol  (as  detailed  elsewhere  [18])  was  used  for  data  correction.  The  effects  of 
specimen  plasticity  are  accounted  for  using  the  J-integral  procedure  to  calculate  a  total  stress  intensity 
(Kj).  Remaining  ligament  plasticity  was  evaluated  via  analytical  solutions  for  a  stationary  crack  in  a 
material  that  deforms  according  to  the  Ramberg-Osgood  flow  rule  and  the  effective  crack  length  (aefr, 
taken  as  the  physical  crack  length  plus  the  plain  strain  plastic  zone  correction).  Full  details  of  this 
approach  are  presented  elsewhere  [49]. 

Testing  and  material  response  complexities  can  cause  deviations  from  the  classic  da/dt  vs.  K 
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behavior,  i.e.  an  initial  absence  of  crack  growth  below  KTh  that  is  followed  by  a  steep  rise  in  da/dt  with 
increasing  K,  ultimately  reaching  a  K-independent  (Stage  II)  growth  rate  (da/dtn).  As  such,  it  is  necessary 
to  assign  conditions  for  determining  KTh  and  da/dtn.  For  the  current  study,  KTH  is  defined  as  the  K  at 
which  the  growth  rate  exceeds  the  plasticity-based  resolution  limit  of  the  dcPD  system  (as  discussed 
in  [18])  and  the  Stage  II  crack  growth  rate  is  defined  as  the  da/dt  at  K  =  60  MPaVm. 

2.2.2  Testing  Environment 

Three  test  environments  were  investigated  at  23°C:  full  immersion  in  an  aqueous  chloride 
electrolyte  with  an  applied  potential  of  either  -950  or  -850  mVscE  and  dry  N2  gas  (RH  <  5%).  Prior 
efforts  suggest  that  the  HEAC  susceptibility  in  these  environments  is  aggressive,  mild,  and  inert, 
respectively  [18].  To  contain  the  environment,  the  gauge  section  of  the  SEN(T)  specimens  was  placed 
inside  a  340  mL  cylindrical  Plexiglass  cell.  Immersion  testing  was  completed  by  circulating  0.6  M  non- 
deaerated  NaCl  solution  from  a  4L  reservoir  at  ~20  mL/min  in  ambient  conditions.  Testing  was 
completed  at  constant  applied  potentials  of  -850  and  -950  mV,  referenced  to  a  standard  calomel 
electrode  (SCE)  controlled  by  a  potentiostat  operated  in  floating  mode  (to  avoid  a  ground  loop).  The 
SEN(T)  specimen  was  grounded  through  grip  attachment  to  the  testing  machine.  A  platinum-coated  Nb 
mesh  counter-electrode  and  the  SCE  were  placed  in  the  cell,  with  the  counter-electrode  surrounding  the 
specimen.  Prior  to  onset  of  mechanical  testing,  the  open  circuit  potential  of  the  specimen  was  monitored 
for  1  hour,  followed  by  the  specimen  being  polarized  to  the  selected  test  potential.  The  bulk  solution  pH 
was  found  to  rise  from  6.0  to  between  6.2-6. 5  during  an  experiment.  Onl>  a  I  mm  window  surrounding 
the  Mode  I  crack  path  from  the  notch  to  an  a/W  of  0.7  was  exposed  to  solution;  all  other  surfaces  within 
the  cell  were  masked  with  3M  Electroplating  Tape  470  and  a  butyl  rubber-based  lacquer  (Tolber  Miccro 
XP-2000  Stop-off  Lacquer).  Gas  testing  was  completed  by  flowing  N2  into  the  environmental  cell  at  a 
rate  that  maintained  a  measured  relative  humidity  of  less  than  5%  for  the  duration  of  loading.  Upon  test 
completion,  each  specimen  was  rinsed  with  deionized  water,  sonicated  in  acetone  and  methanol  for  10 
minutes  each,  and  then  stored  in  a  desiccator. 

2.3  Characterization 

Scanning  electron  microscopy  (SEM)  was  performed  to  identify  the  final  precrack  length  and  to 
examine  the  fracture  morphology.  Orientation  imaging  microscopy  (OIM)  was  completed  for  each  lot 
using  electron  backscatter  diffraction  (EBSD)  with  a  30  kV  accelerating  voltage  and  a  7.6  nA  probe 
current.  EBSD  samples  were  extracted  from  the  non-gauge  portion  of  the  SEN(T)  •  specimens, 
mechanically  polished  to  0.10  pm,  and  then  vibratory  polished  using  0.05  pm  alumina  suspension.  A  0.5 
pm  step  size  was  employed  to  ensure  >10  steps  for  the  smallest  grains  observed.  The  resulting  OIM  maps 
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were  then  analyzed  using  the  Tango  post-processing  software  (Oxford  Instruments)  to  evaluate  grain 
boundary  misorientation  character  and  average  grain  size. 

Hydrogen  interactions  were  characterized  through  determination  of  hydrogen  evolution  reaction 
kinetics,  CH>diff,  CH,tot>  and  trap  binding  energy  using  well-developed  procedures  [25,50,51].  Specifically, 
the  hydrogen  evolution  reaction  kinetics  (HER)  were  studied  through  electrochemical 
interrogation  [25,50-52]  and  the  total  and  diffusible  hydrogen  concentrations,  hydrogen  diffusion  rates 
(Deff),  and  trap  binding  energies  were  determined  using  electrochemical,  thermal  emission  (LECO)  and 
thermal  desorption  methods,  respectively  [51]. 

Results 

3.1  Quantification  of  Cracking  Behavior 

3.1.1  Inert  Environment  Testing 

The  da/dt  versus  Kj  results  for  Allvac,  TR1,  TR2,  TR3,  and  duplicate  NRL-HS  tests  in  a  dry 
nitrogen  gas  environment  (RH  <  5%)  are  shown  in  Figure  1.  Consistent  with  prior  results  for  Allvac 
material  [18],  each  lot  demonstrated  two  regimes  of  behavior.  First,  da/dt  was  found  to  increase  linearly 
with  increasing  Kj  over  a  range  of  low  Kj  (10  <  Kj  <  35  MPaVm).  This  behavior  has  been  attributed  to  a 
crack  tip  plasticity-induced  change  in  resistivity  that  scales  with  applied  K  and  increases  the  dcPD- 
measured  potential  [18].  Thus,  it  is  possible  to  establish  a  false  daMt  '‘resolution  limit”,  as  discussed 
elsewhere  [18],  which  is  described  by  the  following  relationship  for  Allvac  tested  at  a  dK/dt  of  1 
MPaVm/hr  (solid  line  in  Figure  1):  log[da/dtFALSE  (mm/s)]  =  -6.3009  +  0.0127  Kj  (MPaVm).  Since  the 
plastic  zone  size  depends  on  both  K  and  yield  strength  [53],  it  is  expected  that  the  magnitude  of  false 
dcPD  rise  due  to  crack  tip  plasticity  effects  will  be  independent  of  the  applied  loading  rate.  Therefore,  the 
resolution  limit  should  scale  linearly  with  applied  dK/dt,  with  a  dK/dt  of  0.33  MPaVm/hr  simply  being  a 
three-fold  reduction  of  the  1  MPaVm/hr  resolution  limit.  Despite  lot-to-lot  differences  in  yield  strength, 
this  is  shown  to  be  reasonably  accurate  by  comparing  tests  conducted  at  a  constant  dK/dt  of  0.33 
MPaVm/hr  in  dry  nitrogen  with  the  predicted  0.33  MPaVm/hr  resolution  limit  (dashed  line  in  Figure  1). 
Subsequent  plots  will  report  the  resolution  limit  specific  to  the  dK/dt  used  for  data  generation. 

The  second  regime  of  behavior  is  typified  by  crack  growth  above  the  plasticity-induced  resolution 
limit  at  K,  values  greater  than  35  MPaVm/hr.  For  all  material  lots,  SEM  fractography  confirmed  100-200 
pm  of  real  sub-critical  transgranular/ductile  crack  extension  isolated  to  the  center,  high-constraint  portion 
of  the  specimens,  as  shown  for  inert  environment-tested  TR2  in  Figure  2a.  This  extension  occurs  well 
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below  the  reported  ranges  of  K1C  for  Monel  K  -500  (198-340  MPaVm  [18])  and  significantly  exceeds  the 
~25  pm  that  would  be  expected  for  crack  tip  blunting  (calculated  at  a  Kj  of  80  MPaVm  [18]).  The 
dimpled,  ductile  crack  progression  in  the  center  of  the  specimen  (Figure  2b,  right)  was  distinct  from  the 
morphology  in  the  precrack  region  (left).  High  magnification  imaging  of  the  dimple  features  associated 
with  the  slow  sub-critical  cracking  on  the  Allvac  and  TR2  fracture  surfaces  (Figure  2c  and  2d, 
respectively)  demonstrates  that  the  microvoid-based  failure  initiated  at  precipitates.  Though 
compositional  analysis  of  the  precise  particle  composition  via  energy-dispersive  X-ray  spectroscopy 
(EDS)  is  complicated  by  the  particle  size  being  less  than  the  probe  interaction  volume,  comparing  EDS 
data  from  the  particles  to  that  of  the  base  material  demonstrated  a  clear  increase  in  Ti  and  C,  suggesting 
that  these  particles  are  TiC-type  carbides  [43,44].  Dey  et  al.  observed  a  roughly  300  nm  y’-precipitate 
free  zone  (PFZ)  proximate  to  the  TiC  precipitates  that  is  facilitated  by  local  Ti  and  vacancy 
depletion  [43].  It  has  been  postulated  that  the  observed  sub-critical  ductile  crack  advance  is  governed  by 
a  thermally-activated  dislocation  glide  mechanism  in  the  PFZ  [54].  Analysis  of  a  deformation- 
mechanism  map  for  pure  Ni  suggests  that,  at  ambient  temperature  (23 °C)  and  high  stresses  (as  pertinent 
to  a  crack  tip),  enhanced  dislocation  glide  (which  is  not  true  time-dependent  creep  behavior)  is 
active  [55].  This  enhanced  dislocation  mobility  becomes  increasingly  important  as  the  flow  stress 
decreases  with  loading  rate.  Furthermore,  this  mechanism  may  be  augmented  by  an  H-contribution  from 
either  residual  internal  H  or  low  levels  of  crack  tip  production  due  to  residual  moisture  in  the  dry  N2 
and/or  H  production  at  mild  cathodic  polarizations  [54]. 

3.1.2  Aqueous  Chloride  Electrolyte  Testing 

3. 1.2.1  Applied  Potential  of -950  mVScE 

The  da/dt  versus  K;  data  for  Allvac,  TR1,  TR2,  TR3,  and  duplicate  NRL-HS  tests  in  0.6  M  NaCl 
at  an  applied  potential  of -950  mVScE  are  presented  in  Figure  3.  The  Allvac,  TR  1,  TR  2,  TR  3,  and 
NRL-HS  material  lots  show  measured  KTH  values  of  29,  14,  18,  19,  and  15/17  MPaVm,  respectively.  All 
Kth  values  are  closely  clustered  together  for  this  testing  condition,  with  the  exception  of  Allvac,  whose 
higher  threshold  value  suggests  an  increased  resistance  to  HEAC.  This  threshold  stress  intensity  is  similar 
to  a  previous  reported  value  of  30  MPaVm  at  -934  mVscE  for  the  Allvac  lot  [18].  Stage  II  crack  growth 
rates  (at  K  =  60  MPaVm)  of  7.2xl0'6,  3.8xl0'6,  2.7xl0'6,  4.6xl0'6,  and  2.6xl0’6/2.8xl0'6  mm/s  were 
observed  for  Allvac,  TR  1,  TR  2,  TR  3,  and  NRL-HS,  respectively.  The  variation  by  a  factor  of  3 
between  lots  is  within  the  expected  growth  rate  scatter  established  via  five  replicate  tests  of  a  single  lot  of 
Monel  K-500  in  0.6  M  NaCl  at  -1000  mV;  specifically,  a  4-fold  variation  in  da/dtn  was  observed  [56]. 
The  current  growth  rates  are  lower  than  the  previously  reported  da/dtn  of  1.9  to  2.3xl0'5  mm/s  at  -1000 
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mVscEand  l.lxl O'5  mm/s  at  -900  mVScE  for  Allvac  [18].  However,  the  current  da/dtn  data  is  reported  at 
Kj  of  60  MPaVm  for  a  dK/dt  of  0.33  MPaVm/hr,  whereas  prior  work  reported  values  for  Kj  of  50  MPaVm 
for  a  dK/dt  of  1.1  MPaVm/hr  [18].  While  observed  slower  growth  rates  for  Allvac  in  the  current  work 

t 

may  be  attributed  to  test-to-test  scatter,  they  are  directionally  consistent  with  prior  results  in  H-producing 
environments  that  show  increasing  dK/dt  can  enhance  HEAC  susceptibility  in  Ni-based  alloys  [18], 
steels  [57],  and  Ti-6A1-4V  [58]. 

Fractography  (each  at  roughly  a  K  of  40  MPaVm)  shown  in  Figure  4  demonstrates  that  flat  IG 
cracking  is  the  dominant  fracture  morphology  for  each  material  during  slow-rising  displacement  testing  in 
0.6M  NaCl  at  an  applied  potential  of -950  mVScE»  However,  NRL-HS  also  showed  distinct  transgranular 
(TG)  features,  along  with  dimpling  on  some  IG  facets  (Figure  4e-f).  Each  image  was  taken  at  the  same 
magnification,  therefore  different  feature  sizes  reflect  differences  in  grain  size  between  material  lots. 
High  magnification  images  (Figure  5)  reveal  slip  traces  on  the  fracture  surface  IG  features  for  each 
material.  A  detailed  study  of  the  density,  spacing,  and  surface  step  height  was  not  performed  for  each 
specimen,  but  a  coarse  observation  suggests  that  the  features  were  generally  similar  between  materials. 
The  slip  behavior  for  Monel  K-500  (and  the  resulting  strength  and  hardening  behavior)  is  largely 
governed  by  the  size  of  the  homogeneously-distributed  y’  precipitates  and  is  critically  dependent  on  the 
processing  history  [43,44].  Furthermore,  prior  efforts  have  shown  a  correlation  between  the  spacing  of 
the  slip  traces  on  the  IG  fracture  feature  with  TEM  characterization  of  the  dislocation  cell  structure 
underneath  the  fracture  surface,  suggesting  a  potentially  critical  effect  of  bulk  slip  behavior  on  the 
controlling  HEAC  damage  mechanism  [31,59].  Despite  the  differences  in  the  yield  strength  reported  in 
Table  2,  this  coarse  observation  of  similarity  in  the  IG  facet  slip  traces  does  not  suggest  a  strong  variation 
in  the  y’  and  slip  interaction  (i.e.  the  relative  amount  of  shearing  and  looping).  A  detailed  study  of  the 
effect  of  y’  morphology  on  the  HEAC  behavior  is  ongoing,  but  in  the  context  of  the  current  effort,  similar 
bulk  slip  behavior  amongst  the  material  lots  is  inferred  from  the  morphology  in  Figure  5  and  similar 
hardening  behavior  (Table  2). 

3. 1.2.2  Applied  Potential  of -850  mVScE 

Severe  H-producing  environments  will  give  rise  to  sufficiently  high  crack  tip  H-concentrations 
(CH-o)  such  that  any  effect  of  lot-to-lot  variation  in  metallurgy  on  HEAC  susceptibility  may  be  masked  by 
the  strong  propensity  for  IG-HEAC.  Therefore,  it  is  advantageous  to  investigate  the  effect  of  metallurgy- 
based  variations  in  the  HEAC  behavior  by  testing  at  applied  potentials  that  produce  a  lower  level  of  CH. 
diff.  Prior  work  established  an  empirical  relationship  between  diffusible  hydrogen  concentration  and  the 
applied  potential  for  Monel  K-500  by  utilizing  artificial  crevice  geometries  pertinent  to  current  LEFM 
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experiments.  These  experiments  demonstrated  that,  for  an  applied  potential  of -850  mVSCE,  the  value  of 
CH-Difr generated  in  the  crack  tip  environment  is  <10  wppm  [18,25].  As  such,  the  HEAC  kinetics  were 
established  at  this  less  negative  potential  (-850  itiVsce)  in  0.6  M  NaCl  for  each  material  lot. 

The  da/dt  versus  Kj  data  for  duplicate  Allvac,  TR1,  TR2,  TR3,  and  duplicate  NRL-HS  tests  in  0.6 
M  NaCl  at  an  applied  potential  of -850  mVScE  are  presented  in  Figure  6.  The  TR1,  TR  2,  and  TR  3 
material  lots  all  had  similar  da/dtn  crack  growth  rates  of  2.2X10-6  mm/s,  4.0xl0’6  mm/s,  and  3.8xl0'6  mm/s 
at  a  K  of  60  MPaVm  with  KTH  values  of  1 7,  42,  and  24  MPaVm,  respectively.  Macroscale  fractography  of 
TR  2  shows  the  smooth  EDM  notch  (Figure  7a,  far  left),  the  in-air  fatigue  precrack  region  that  consisted 
of  TG-slip  based  cracking  ending  at  the  dotted  line,  transitioning  to  mixed  IG/TG  crack  progression  along 
the  entire  crack  front  during  the  slow-rising  displacement  testing  to  the  solid  line,  followed  by  final 
ductile  failure.  The  mixed  IG/TG  morphology  observed  during  loading  at  an  applied  potential  of -850 
mVScE  in  0.6  M  NaCl  is  shown  at  a  K  of  roughly  40  MPaVm  in  Figure  7b,  7c,  and  7d,  for  TR1,  TR  2,  and 
TR  3,  respectively.  This  morphology  is  consistent  with  slow-rising  displacement  testing  of  Monel  K-500 
where  mild  cathodic  polarizations  resulted  in  identical  IG/TG  crack  growth  [18].  The  da/dt  vs.  K  data  for 
Allvac  (replicate  tests)  and  NRL-HS  (replicate  tests)  specimens  show  growth  rate  behavior  that  is  similar 
to  the  inert  dry  N2  results,  as  demonstrated  in  Figure  8a  for  NRL-HS.  Furthermore,  fractography 
confirmed  the  presence  of  the  same  macroscale  morphology  for  NRL-HS  and  Allvac  at  -850  mVScE  (e.g. 
Figure  8b)  as  was  observed  for  inert  dry  N2  testing  (e.g.  Figure  2a).  A  slight  variations  in  morphology 
was  observed  between  Allvac  (small  dimpled  features)  and  NRL-HS  (dimpled  ridge-like  structure,  Figure 
8c)  at  higher  magnifications.  Speculatively,  this  variation  is  attributed  to  an  increased  level  of  hydrostatic 
stress  in  NRL-HS  causing  secondary  cracks  to  form  perpendicular  to  the  crack  plane.  Regardless,  these 
findings  suggest  that  no  IG-HEAC  occurred  in  these  specimens,  therefore  KTh  and  da/dtn  are  not  reported 
for  the  Allvac  and  NRL-HS  lots.  The  inconsistency  in  the  HEAC  susceptibility  amongst  material  lots  can 
be  interpreted  as  either  (1)  variability  that  is  inherent  to  testing  in  a  potential  regime  where  minor  changes 
in  local  crack  environment  result  in  orders  of  magnitude  changes  in  HEAC  susceptibility  [18],  or  (2) 
metallurgy-specific  characteristics  that  shift  the  critical  potential  where  HEAC  susceptibility  initiates. 
The  former  is  treated  in  a  companion  paper  [52]  and  the  latter  will  be  developed  in  the  Discussion. 

3.2  Characterization  of  the  Grain  Structure 

Orientation  maps  colored  according  to  the  indicated  inverse  pole  figure  are  presented  for  each 
material  in  Figure  9a-e  with  grain  boundaries  highlighted  as  black  lines.  The  results  of  an  ASTM  El  12- 
113  grain  size  analysis  are  reported  in  Table  3,  along  with  the  calculated  grain  boundary  surface  area  to 
volume  ratio  (Sv;  mm2/mm3)  [60,61].  For  polycrystalline  materials,  the  IG-HEAC  behavior  is  influenced 
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by  the  prevalence  of  certain  special  boundaries  (low  angle  and  coincident  site  lattice  boundaries  (CSL)) 
and  the  degree  of  high  angle  boundary  connectivity  [62-64].  The  OIM  of  a  representative  0.05  mm2  area 
of  each  material  is  used  to  quantify  the  distribution  of  high  angle  (>15°),  low  angle  (>5°  and  <15°),  and 
CSL  boundaries.  The  resolution  of  the  analysis  enabled  characterization  of  orientation  variations  within 
grains  (attributed  to  sub-grain  crystal  distortions  and  dislocation  structures  [65]),  though  for  clarity  in  the 
intergrain  misorientation  analysis  and  grain  size  calculations,  the  orientation  of  each  grain  was  averaged 
using  the  Tango  software  and  misorientation  angles  below  5°  were  not  included  in  the  analysis.  The 
fraction  of  (1)  the  total  number  of  boundaries  and  (2)  the  total  Sv  that  are  low  energy  (including  random 
boundaries  <15°  and  CSL  boundaries)  range  from  53-61%  and  53-64%,  respectively  (Table  3),  with  the 
vast  majority  of  low  energy  boundaries  being  S3  boundaries.  Despite  this  seemingly  small  range,  the 
differences  in  grain  size  give  rise  to  a  drastic  variation  in  the  total  Sv  of  low  energy  boundaries  between 
material  lots,  which  may  be  important  in  the  context  of  the  susceptible  path  connectivity  [62,63].  For 
example,  43  and  44%  of  the  Sv  are  high  energy  boundaries  in  Allvac  and  TR  2,  but  grain  size  differences 
(13.8  and  35.3  pm,  respectively)  result  in  the  total  Sv  of  high  energy  boundaries  being  162  mm2/mm3  for 
Allvac  and  70  mm2/mm3  for  TR2. 

3.3  Quantification  of  Hydrogen  Interaction  Parameters 

Evaluating  the  effect  of  hydrogen  interactions  requires  the  examination  of  many  parameters. 
These  include  hydrogen  evolution  kinetics  to  address  differences  in  H  production  and  adsorption, 
hydrogen  concentration,  hydrogen  diffusion  rates,  and  trap  binding  energies.  The  segregation  of  trace 
elements  to  the  grain  boundary  is  also  of  interest  due  to  their  potential  to  promote  increased  hydrogen 
activity  and/or  lower  CH.crit  [66],  as  well  as  lower  intrinsic  fracture  strength  or  toughness  [25,67,68].  The 
kinetics  of  the  hydrogen  evolution  reaction  were  studied  by  measuring  the  Tafel  slope  and  exchange 
current  densities  for  water  reduction  (hydrogen  evolution)  in  0.6  M  NaCl  at  pH  8.0.  Production  rates  were 
similar  among  the  different  materials;  lot-to-lot  variation  does  not  affect  this  parameter.  Figure  10 
displays  the  effect  of  hydrogen  overpotential  on  hydrogen  concentration  where  both  concentration 
parameters  were  found  to  increase  with  increasing  hydrogen  overpotential.  As  expected,  the  total 
hydrogen  concentration  CHj0tai  was  larger  than  C^Difr  [25].  In  terms  of  lot-to-lot  variation,  there  was 
limited  difference  between  both  Cnjotai  and  CH,Difr,  except  for  the  Allvac  material,  which  consistently 
exhibited  hydrogen  concentration  near  the  lower  bound  of  all  other  heats  or  statistically  lower  hydrogen 
levels  based  on  repeat  measurements  at  some  potentials.  Variation  in  diffusible  hydrogen  concentrations 
could  arise  from  hydrogen  interacting  with  low  and  high  binding  energy  traps,  such  as: 
dislocations  [69,70],  grain  boundaries,  carbides  [71],  microvoids,  age  hardening  precipitates  (e.g. 
y’)  [25,69,72],  among  others.  Two  trap  binding  states  were  found  in  all  the  heats  tested;  these  results  are 
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in  agreement  with  previous  information  reported  on  solution  heat  treated  (SHT),  aged  and  SHT,  cold 
worked,  aged  Monel  K-500  materials  [25].  The  average  values  for  the  binding  energies  for  y’  were  found 
in  the  range  7.8  to  17.7  kJ/mol  in  the  following  ascending  order:  NRL-HS  <  TR  2  <  Allvac.  Changes  in 
trap  binding  energy  would  result  in  lot-to-lot  differences  in  the  hydrogen  diffusion  coefficient,  which 
would  critically  affect  the  proposed  HEAC  behavior  (as  quantified  in  Eqn  3).  However,  prior  analysis 
found  similar  hydrogen  diffusion  coefficients  (in  the  range  0.9-3.9  x  10’14  m2/s)  for  all  materials, 
suggesting  a  secondary  role  of  trapping/diffusion  on  lot-dependent  HEAC  behavior  [25,52]. 

Discussion 

The  crack  growth  kinetics  (Figure  1,  3,  6)  and  cracking  morphologies  (Figure  4,  8)  demonstrate  that, 
in  both  inert  (dry  N2  gas)  and  aggressive  (0.6  M  NaCl  polarized  to  -950  itiVsce)  environments,  the  sub- 
critical  cracking  behavior  shows  little  variation  in  HEAC  susceptibility'  between  different  lots  of  Monel 
K-500.  However,  at  mildly  aggressive  environments  (0.6  M  NaCl  polarized  to  -850  itiVsce)  where  low 
crack  tip  H-uptake  is  expected  [1 8],  a  noticeable  change  in  fracture  morphology  is  observed.  Specifically, 
the  Allvac  and  NRL-HS  lots  demonstrate  a  ductile  transgranular  (TG)  cracking  morphology  in  the  high 
constraint  center  of  the  specimen,  while  TR1,  TR2,  and  TR2  show  both  a  mixed  TG/IG  morphology  and 
crack  growth  kinetics  consistent  with  a  HEAC-based  mechanism.  This  discussion  will  augment  the 
fracture  mechanics  testing  results  with  detailed  materials  characterization  data  and  mechanism-based 
modeling  to  provide  insight  into  three  primary  questions:  (1)  Can  the  potential  dependence  of  the 
observed  lot-to-lot  variation  be  rationalized  in  the  context  of  a  hydrogen  embrittlement  mechanism?  (2) 
To  what  extent  do  various  material  lot-dependent  characteristics  (yield  strength,  composition,  grain 
boundary  character,  H-metal  interactions,  and  impurity  segregation)  correspond  with  the  observed 
changes  in  the  HEAC  susceptibility  for  mildly  aggressive  environments?  (3)  Can  the  differences  in 
microstructure  be  incorporated  into  micromechanical  modeling  to  capture  the  observed  variation  in 
HEAC  susceptibility? 

4.1  Effect  of  Microstructure  Variations  at  Mild  and  Aggressive  Applied  Potentials 

Decohesion-based  micromechanical  models  (Eqn.  2  and  3)  provide  a  quantitative  framework  for 
understanding  the  extent  to  which  variations  in  material  properties  will  influence  HEAC  metrics  (KTh  and 
da/dtn)  at  different  applied  potentials.  Specifically,  this  framework  can  be  utilized  to  inform  why  all 
material  lots  demonstrate  similar  HEAC  susceptibility  at  -950  mVScE,  while  significant  differences  are 
observed  for  materials  tested  at  -850  mVSCE-  The  justification  and  validation  of  the  micromechanical 
modeling  approach  was  established  in  prior  work  [18]  and  has  been  recently  refined  [54].  For  current 
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modeling,  the  input  parameters  for  KTh  modeling  are  oY s  =  786  MPa  (Allvac),  a”  =  0.0002  MPaVm,  and 
P’  =  0.20  (MPaVm)'1  (based  on  dislocation-shielding  simulations  of  the  near  crack  tip  stress  field  [32,73] 
and  validated  by  experimental  data  on  IN718  [37,54]).  CH-CT  is  calculated  from  Eqn.  1,  where  oH  =  9oYs 
and  an  empirical  relationship  relates  the  applied  potential  to  CH-Difffor  the  Allvac  material  [18].  kiG  is  the 
Griffith  toughness  of  the  IG  interface,  calculated  to  be  0.88  MPaVm  using  the  Griffith  relation  and 
theoretical  values  of  the  surface  energy  of  Ni  [54].  a  represents  the  potency  of  the  local  H  in  decreasing 
kiG  and  is  set  to  4.078  (MPaVm  per  atom  fraction  H)  to  fit  the  model  to  the  average  of  three  Allvac  KTh 
values  (established  from  previous  work  as  18  MPaVm  [18])  at  an  applied  potential  of -1000  mVSCE.  The 
model-predicted  effect  of  applied  potential  on  the  K^of  Allvac  is  reported  in  Figure  1  la-c  (solid  black 
lines)  using  the  parameters  detailed  above. 

The  relative  influence  of  microstructural  variation  at  different  applied  potentials  can  be  evaluated 
by  systematically  varying  input  parameters  that  may  exhibit  lot-to-lot  differences.  Though  details 
regarding  the  link  between  the  modeling  parameters  and  microstructure  features  will  be  established  later 
in  the  discussion,  a  sensitivity  analysis  is  useful  to  demonstrate  the  effect  of  such  changes  on  HEAC 
behavior.  Three  modeling  parameters  are  varied  in  the  current  KTh  modeling.  First,  the  sensitivity  of  the 
models  to  variations  in  H  uptake  behavior  is  evaluated  by  multiplying  and  dividing  the  calculated  CH-Difr 
by  2.5  (corresponding  to  a  250%  increase  and  a  60%  decrease).  This  250%  increase  in  uptake  behavior  is 
physically  realistic  for  Monel  K-500  [18]  and  is  supported  by  the  present  analysis  (Figure  10)  which 
indicates  lower  CH-Diff  in  Allvac  for  a  given  applied  potential  [52].  Figure  11a  clearly  demonstrates  that 
increasing/decreasing  CH-Diff  results  in  a  more  pronounced  variation  in  KTH  as  the  applied  potential  is 
increased  from  -950  to  -850  mVSCE.  For  example,  increasing  CH-Difr  by  250%  (Figure  1  la,  thick  solid 
blue  line)  results  in  a  decrease  of  KTh  by  61  MPaVm  at  -850  mVScE  whereas  a  decrease  of  only  15 
MPaVm  is  calculated  for  -950  mVScE-  Similarly,  the  changes  in  KrH  are  again  more  severe  at  -850  mVScE 
when  the  intrinsic  grain  boundary  toughness  (k1G)  is  varied  ±25%  with  all  other  inputs  at  the  baseline 
levels,  as  shown  in  Figure  1  lb.  Physically,  such  variations  in  grain  boundary'  toughness  could  result  from 
differences  in  the  fraction  of  high  energy  boundaries  or  boundaries  containing  significant  concentration  of 
deleterious  impurity  atoms.  Decreasing  kIG  from  0.88  to  0.66  MPaVm  (Figure  1  lb,  thick  solid  blue  line) 
results  in  a  KTH  decrease  of  62  and  14  MPaVm  for  -850  and  -950  mVScE,  respectively.  Finally,  the  -850 
mVscE  condition  is  significantly  more  sensitive  to  a  ±15%  change  in  oYs  (Figure  11c,  thick  solid  blue 
line);  an  increase  in  oYS  from  786  to  904  MPa  results  in  a  65  and  16  MPa  change  in  KTH  for  -850  and 
-950  mVscEj  respectively.  Critically,  this  variation  in  yield  strength  represents  the  difference  between  the 
Allvac  and  NRL  HS  lots,  as  shown  in  Table  2,  illustrating  the  large  change  in  HEAC  resistance  that  could 
be  expected  in  this  alloy. 
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A  similar  micromechanical  modeling-based  sensitivity  analysis  was  performed  to  understand  how 
changes  in  microstructure  may  influence  da/dtn  at  different  applied  potentials.  The  baseline  parameters 
for  current  modeling  mimic  those  in  prior  work  [18,54],  where  Defr=  lxlO'10  cm2/s,  Xcrit=  1  pm,  CH-0is 
calculated  in  the  same  manner  as  in  the  KTh  modeling,  and  CH-critis  set  to  647  wppm  [54].  The  baseline 
model  (Figure  12a-b;  black  solid  line)  demonstrates  HE  AC  growth  rates  increase  several  orders  of 
magnitude  as  the  potential  is  decreased  from  approximately  -840  to  900  mVScE;  however,  further 
reduction  of  the  potential  from  -900  to  -1 100  mVScE  results  in  a  minimal  increase  in  da/dt.  This  behavior 
is  understood  in  the  context  of  the  abundant  supply  of  crack-tip  H  at  the  low  potentials,  whereas  at  higher 
potentials  there  is  insufficient  supply  of  H  to  enable  embrittlement.  The  sensitivity  of  this  behavior  to 
varying  metallurgical  features  can  be  quantitatively  evaluated  by  varying  the  hydrogen  uptake  behavior 
and  yield  strength  inputs  in  Eqn  3;  the  Deffand  x^t  parameters  are  considered  invariant  with  regards  to 
microstructure  [18,25].  To  evaluate  the  sensitivity  of  the  model  to  variations  in  hydrogen  uptake,  the 
diffusible  hydrogen  concentration  was  again  increased  250%  and  decreased  60%  (Figure  12a-b,  thick 
solid  blue  line  and  dashed  line,  respectively).  As  demonstrated  by  Figure  12a,  the  difference  in  da/dtn 
significantly  increases  amongst  the  three  calculated  curves  when  moving  from  -950  to  -850  mVScE> 
suggesting  that  the  effect  of  metallurgical  variations  is  more  discernible  at  higher  potentials.  This  same 
behavior  is  noted  when  varying  the  yield  strength  ±15%;  specifically,  the  effect  of  yield  strength  variation 
on  da/dtn  is  more  pronounced  at  less  negative  applied  potentials.  These  observations  can  be  rationalized 
in  the  context  of  strong  cathodic  polarizations  leading  to  a  level  of  crack-tip  H  that  greatly  exceeds  the 
critical  value  required  for  cracking  (i.e.  CHo-cnt)*  Therefore,  any  variation  in  material  HEAC  resistance 
will  only  marginally  influence  the  cracking  behavior  at  -950  mVscE,  while  HEAC  susceptibility  becomes 
more  material-dependent  as  the  applied  potential  increases  to  -850  mVScE- 

In  toto ,  a  sensitivity  analysis  using  previously  validated  models  for  KTHand  da/dtn  (Eqn  2  and  3, 
respectively)  clearly  establishes  that  microstructure-based  changes  in  material  properties  are  expected  to 
have  an  increasing  effect  on  HEAC  susceptibility  at  applied  potentials  that  produce  lower  crack-tip  H 
concentrations.  These  results  are  consistent  with  experimental  observations  of  minimal  lot-to-lot  variation 
in  the  HEAC  behavior  at  -950  mVScE  (Figure  3-5)  and  significant  lot-dependent  HEAC  behavior  at  -850 
mVscE  (Figure  6-8).  Taken  together,  these  experimental  and  micromechanical  modeling  results 
demonstrate  that  microstructural  influences  can  alter  HEAC  susceptibility  with  the  caveat  that  these 
deviations  are  attenuated  when  environmental  conditions  become  severe.  Furthermore,  these  finding 
reinforce  the  need  to  accurately  capture  environmental  conditions  when  performing  engineering-scale 
modeling  so  as  to  ensure  reliable  predictions. 

4.2  Correlation  of  Microstructure  Features  with  HEAC  Susceptibility 
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The  susceptibility  of  Mi-based  alloys  to  IG  H-embrittlement  has  been  found  to  scale  with  yield 
strength  [74],  grain  size  [75-77],  grain  boundary  (GB)  character  [30,62,78],  and  composition  [79-81]. 
Isolating  the  relative  influence  of  each  feature  is  challenging  for  complex  engineering  alloys,  especially 
when  the  material  processing  history  is  uncertain.  Quantitative  characterization  of  each  characteristic 
enables  correlation  of  measured  properties  with  observed  HEAC  behavior,  thereby  providing  an 
opportunity  to  elucidate  dominant  features.  The  following  section  seeks  to  connect  the  measured  variation 
in  material  properties  with  the  observed  HEAC  susceptibility  found  for  each  material  lot. 

4.2.1  Material  Composition 

The  basic  composition  of  Monel  K-500  consists  of  roughly  66  wt%  Ni  and  30  wt%  Cu, 
augmented  by  the  addition  of  -3  wt%  A1  and  0.6  wt%  Ti  to  enable  further  strengthening  via  precipitation 
of  the  y’-phase  [82].  Trace  impurity  elements  can  have  significant  effects  on  creep  resistance, 
weldability,  corrosion  resistance,  and  mechanical  properties  in  Ni-based  alloys  and  must  be  controlled  to 
optimize  material  performance  [79,80,83].  Specifically,  impurity  elements  (e.g.  O,  H,  N,  Ar,  He,  S,  P, 
Pb,  Bi,  Sb,  As,  Se,  Ag,  Cu,  Tl,  Te)  have  been  broadly  characterized  as  deleterious  to  the  material 
properties  of  Ni  and/or  Ni-based  alloys  [79,84].  For  the  current  work,  no  clear  trend  of  trace  S,  P,  Pb,  Ag, 
or  Sb  content  and  HEAC  resistance  was  observed,  though  higher  levels  of  As  (and  marginally  Bi  and  Te) 
were  found  in  the  more  susceptible  TR  1,  TR2,  and  TR3  materials.  T  hese  subtle  differences  are  likely  not 
sufficient  to  deleteriously  influence  properties,  especially  considering  the  minute  bulk  impurity 
levels  [79,84].  However,  there  is  a  broad  literature  which  demonstrates  that  elements  such  as  Bi,  Te,  Pb, 
S,  B,  C,  Zr,  Hf,  P,  Mg,  and  Ta  have  an  affinity  for  segregating  to  GBs,  giving  rise  to  local  concentrations 
well  above  what  is  observed  in  the  bulk  [80,85,86];  this  segregation  can  magnify  the  influence  of 
variations  in  the  trace  element  content.  The  degree  of  segregation  is  dependent  on  the  element  and  the 
grain  boundary  misorientation  [63,87,88],  as  well  as  the  concentration  of  intentionally  added  getter 
elements  (e.g.  Mg,  Zr,  Hf,  C)  which  retard  the  segregation  of  deleterious  elements  to  the  GB.  Of  critical 
importance  to  the  IG-HEAC  behavior  is  the  (direct  or  indirect)  effect  of  these  impurity  and  trace  elements 
on  the  grain  boundary  susceptibility  to  (chemo-)mechanical  failure. 

While  elemental  C  segregation  to  GBs  has  been  shown  to  be  innocuous  [89]  or  even 
beneficial  [90]  in  Ni  alloys,  Natishan  et  al  suggest  that  a  semi-continuous,  thin  film  of  unalloyed  C  may 
contribute  to  the  IG  failure  of  Monel  K-500  in  hydrogen  [91,92].  This  film  was  sporadically  observed 
along  grain  boundaries  and  is  likely  in  the  form  of  graphite,  similar  to  what  has  been  reported  for  Nickel 
200  after  prolonged  exposures  to  425-650°C  [93].  In  the  context  of  present  results,  Table  1  demonstrates 
that  the  more  susceptible  TR1  and  TR2  have  increased  levels  of  C,  but  this  is  conflicted  by  HEAC- 
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susceptible  TR3  having  less  C  content  than  the  immune  Allvac  lot.  However,  prior  researchers  found  that 
the  prevalence  of  the  film  does  not  scale  with  the  bulk  carbon  content,  but  instead  is  dependent  on  applied 
thermal  processing  [92].  In  this  paradigm,  a  composition/processing  variation  that  failed  to  produce  a 
graphite  film  would  then  be  responsible  for  the  lack  of  IG  cracking  in  Allvac  and  NRL-HS.  Therefore, 
the  mixed  IG/TG  features  observed  for  TR1,  TR2,  and  TR3  (Figure  7)  at  -850  mVScE  could  be  due  to  this 
graphite  film-induced  IG  failure  as  opposed  to  HEAC.  However,  this  conjecture  is  not  consistent  with  the 
ductile/transgranular  morphology  that  is  observed  for  all  alloys  in  the  dry  N2  environment  or  other  low  H 
environment  testing  presented  in  companion  work  [54].  According  to  Natishan  et  al ,  the  graphite  film- 
enhanced  IG  susceptibility  of  Monel  K-500  was  most  pronounced  at  low  H  conditions,  while  inhibiting 
HEAC  susceptibility  at  higher  H  concentrations  [91,92],  which  directly  contradicts  the  current 
experimental  results  (Figures  3-6).  As  such,  this  graphite  film-based  mechanism  is  not  considered  to  be  a 
source  of  the  current  lot-to-lot  variation. 

Sulfur  (where  segregation  can  give  rise  to  GB  concentrations  over  2000-fold  higher  than  bulk 
levels  [79,94,95])  is  widely  considered  to  be  the  dominant  deleterious  segregant  for  Ni-based 
alloys  [64,79,96].  Similar  deleterious  behavior  was  postulated  for  P,  which  also  strongly  segregates  to 
GBs  [79],  however  experimental  [97]  and  theoretical  [89]  results  establish  that  P  is  either  inert  or  even 
beneficial  in  Ni-based  alloys  [80].  Using  either  quantum  mechanical  cluster  calculations  [89,98]  or  a 
local  spin  density  atomic  cluster  model  [99],  researchers  have  postulated  that  the  electronegativity  of  S 
changes  the  local  chemical  bonding  structure  to  decrease  the  GB  cohesive  properties  [64].  Additionally, 
the  embrittling  effect  of  sulfur  has  been  found  to  be  synergistic  with  that  of  H  [66,90,96,100],  which 
suggests  that  lot-to-lot  variations  in  sulfur  segregati on/getter ing  may  impact  HEAC  susceptibility.  There 
is  a  modest  correlation  between  increased  bulk  S  content  (Table  1)  and  HEAC  susceptibility  (Figures  6- 
8),  though  it  is  complicated  by  the  HEAC-resistant  NRL-HS  material  having  5  times  the  S  content  (17 
wppm)  of  the  susceptible  TR1  material  (3.1  wppm).  However,  the  rigor  of  correlating  bulk  S  content  to 
the  GB  cracking  susceptibility  is  compromised  by:  (1)  the  presence  of  getter  elements  that  will  strongly 
influence  the  concentration  of  S  available  for  segregation,  (2)  the  grain  boundary  surface  area  available 
for  segregation,  and  (3)  the  presence  of  trace  elements  that  are  known  to  increase  the  cohesive  strength  of 
the  GB  interface. 

Since  getter  elements  have  a  high  affinity  to  react  with  S  [101],  they  will  preferentially  form 
sulfides  (Mg,  Ca,  Zr  [79,102])  and/or  sulfo-carbides  (Zr,  Hf  [86,103]),  which  reduce  the  S  available  for 
segregation  to  the  grain  boundaries.  Thus,  to  quantify  the  influence  of  these  getter  elements  on  the  HEAC 
resistance,  it  is  critical  to  understand  the  degree  to  which  they  reduce  the  sulfur  available  for  segregation. 
Yamaguchi  et  al  developed  a  relationship  to  describe  the  reduction  of  available  sulfur  for  minor  additions 
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of  getter  elements  (Ca,  Mg,  Y,  and  Zr)  [104];  this  relationship  was  scaled  to  experimental  measurements 
of  hot-workability  and  Auger  analysis  of  grain  boundary  S  for  a  variety  of  Ni-based  superalloys. 
Rigorously  applying  this  model  to  Monel  K-500  to  quantify  the  reduction  of  tramp  residual  S 
concentration  is  complicated  by  alloy  differences  and  focus  on  the  hot-workability  in  model  development. 
However,  the  Yamaguchi  approach  can  be  used  to  qualitatively  evaluate  the  influence  of  the  variation  in 
getter  element  composition  (Table  1)  on  HEAC  resistance.  The  original  relationship  proposed  by 
Yamaguchi  is  modified  in  two  ways  to  calculate  the  effective  gettering  power  (via  a  total  weighted 
composition,  CGet)  of  each  alloy.  First,  Hf  is  included  with  the  same  weighting  factor  as  Zr  based  on  the 
proposed  similarity  of  solubility  product  for  Hf  and  Zr  [101].  Second,  the  Ca  term  is  eliminated  since  the 
Yamaguchi  study  (and  others  that  report  a  gettering  effect  of  Ca  [105])  explicitly  doped  the  .Ni-based 
alloys  with  elemental  Ca,  whereas  the  Ca  concentration  reported  in  Table  1  is  likely  in  the  form  of 
residual  CaO  (or  potentially  CaS)  not  removed  with  the  slag  during  the  desulphurization  process  [106]. 
While  at  molten  conditions,  CaO  will  effectively  remove  elemental  S,  but  CaO  is  chemically  stable  at  all 
known  processing  temperatures  and  is  not  expected  to  act  as  a  getter  of  tramp  residual  S  [107].  The 
resulting  relationship  is: 

C0et=  0.3Mg  T  O.lZr  +  O.lHf  (Eqn4) 

where  elemental  designations  represent  the  respective  bulk  concentrations  in  a  consistent  unit  (e.g.  wt%, 
wppm,  etc).  The  resulting  CGet  values  are  49,  75,  73,  72,  and  77  wppm  for  Allvac,  TR1,  TR2,  TR3,  and 
NRL-HS,  respectively.  Higher  values  of  CGetare  found  to  not  scale  with  decreased  HEAC  susceptibility, 
thereby  suggesting  that  an  additional  influence  outside  of  getter  element  composition  exists. 

For  a  given  bulk  impurity  level,  the  concentration  of  GB  segregants  will  decrease  with  shrinking 
grain  size  due  to  the  increase  in  available  surface  area  per  volume  of  grain  boundary;  this  behavior  has 
been  experimental  established  for  Ni-based  alloys  [76,85].  Normalizing  the  bulk  S  values  listed  in  Table 
1  by  the  lot  specific  GB  surface  area  per  volume  values  (Table  3)  results  in  nominally  lower  values  for 
materials  that  are  more  HEAC  resistant  (0.0046  (Allvac)  and  0.042  (NRL-HS)  wppm/(mm7mm3))  and 
higher  values  for  more  susceptible  materials  (0.074  and  0.075  wppm/(mm2/mm3)  for  TR2  and  TR3, 
respectively).  This  suggests  that  grain  size  may  play  a  tangible  role  in  determining  HEAC  susceptibility, 
even  though  TR1  (0.012  wppm/(mm2/mm3))  does  diverge  from  that  trend.  One  plausible  explanation  for 
this  deviation  is  that  the  effect  of  the  high  S  content  of  NRL-HS  may  be  offset  by  having  the  highest 
effective  gettering  power;  specifically,  the  high  Cs/Sv  would  be  mitigated  by  gettering  elements. 
Regardless,  these  tentative  observations  are  supported  by  preliminary  Auger  analysis  of  fracture  surface 
GB  features  from  H-charged  slow  strain  rate  testing  of  notched  bar  samples,  where  S  enhancement  was 


31 


observed  for  all  materials,  but  grain  size  per  volume  and  gettering  element  concentrations  vary  [52].  A 
more  thorough  Auger  analysis  is  ongoing  to  quantitatively  evaluate  the  interaction  of  bulk  S  composition, 
grain  size,  and  gettering  elements  on  the  deleterious  GB  segregation  of  S. 

GB  segregants  can  also  have  the  beneficial  effect  of  raising  HEAC  resistance.  In  particular,  B  has 
been  experimentally  [108]  and  theoretically  [64,85,89,99,109]  shown  to  increase  the  cohesive  strength 
of  the  GB.  It  has  also  been  postulated  that  Zr  will  have  a  beneficial  intrinsic  effect  on  GB 
cohesion  [103],  however  such  a  beneficial  effect  (above  the  gettering  efficiency  of  Zr)  is  not  widely 
accepted  [79,80].  The  materials  that  have  the  highest  bulk  B  levels  (Allvac;  3.8  wppm,  NRL-HS;  2.0 
wppm)  in  Table  1,  correlate  with  higher  observed  HEAC  resistance  (Figure  6);  this  trend  also  holds  true 
for  bulk  Zr  concentrations  (Allvac;  370  wppm,  NRL-HS;  650  wppm).  However,  no  correlation  is  found 
with  observed  HEAC  resistance  when  the  bulk  B  concentration  is  normalized  by  Sv,  as  the  susceptible 
TR2  and  TR3  lots  have  a  higher  value  than  that  observed  for  NRL  HS.  The  unknown  potency  and  the 
relative  concentration  of  B  that  segregates  to  the  GB  (rather  than  forming  borides  [108])  precludes 
conclusive  and  quantitative  understanding  of  the  extent  to  which  such  small  changes  (0.5-1. 5  wppm)  in 
bulk  B  content  will  impact  the  HEAC  resistance. 

In  toto,  the  compositional  analysis  above  demonstrates  that  the  materials  that  exhibit  enhanced 
HEAC  resistance  (Allvac  and  NRL-HS)  have  the  lowest  bulk  S  content  per  surface  area/volume  of  GB 
and  the  highest  bulk  B  (and  Zr)  levels,  however  there  is  no  clear  trend  with  the  combined  calculated 
gettering  power  of  Mg,  Zr,  and  Hf.  These  qualitative  trend  comparisons  reasonably  suggest  a  potentially 
important  role  of  S  segregation  differences  and/or  segregation  of  GB-strengthening  elements  in 
determining  HEAC  susceptibility.  A  full  compositional  analysis  of  the  observed  GB  particles,  complete 
Auger  analysis  of  the  GB  composition,  and  establishing  the  potency  of  each  trace  element  to  either  getter, 
embrittle,  or  enhance  properties  are  needed,  but  outside  the  scope  of  the  current  effort.  Regardless,  this 
analysis  provides  a  reasonable  basis  to  conclude  that  changes  in  impurity  and  trace  element 
concentrations  likely  contribute  to  (but  do  not  dominate)  the  varying  HEAC  susceptibility  between 
different  lots  of  Monel  K-500. 

4.2.2  Grain  Boundary  Characteristics 

4.2.2. 1  Grain  Size 

The  current  data  (Table  3  and  Figure  9)  demonstrate  that  materials  with  larger  grains  (TRI,  TR2, 
TR3)  are  more  susceptible  to  HEAC  cracking  than  smaller  grained  materials  (Allvac  and  NRL-HS). 
There  are  varying  reports  of  the  effect  of  grain  size  on  HEAC  behavior;  increasing  grain  size  has  been 
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reported  to  decrease  susceptibility  [110,111],  while  others  report  HEAC  behavior  being  independent  of 
size  variation  [112].  Isolation  of  the  intrinsic  effect  of  grain  size  is  complicated  by  the  interdependence 
of  other  properties  with  grain  dimensions;  specifically,  strength,  H-trapping/solubility  behavior  [77,1 13], 
GB  impurity  segregation  [76],  crack  path  dependent  mechanics  [1 14],  and/or  susceptible  crack  path 
connectivity  [111]. 


Stress  intensity  reductions  due  to  increasing  crack  deflection  away  from  the  Mode  I  path  with 
increasing  grain  size  have  been  shown  to  increase  fatigue  cracking  thresholds  and  decreased  crack  growth 
rates  in  the  slow  cracking  (less  than  10"6  mm/  cycle)  [1 14-116].  The  following  relationship  was  proposed 
to  describe  the  effective  stress  intensity  of  a  deflected  crack  (KDef): 


^def 


=  K 


D  cos2(^)+5 
D+5 


(Eqn  5) 


where  D  is  the  distance  over  which  the  tilted  crack  advances,  S  is  the  span  of  the  non-deflected  Mode  1 
crack,  and  0  is  the  angle  of  deflection  [114].  Large  grain  sizes  would  reasonably  lead  to  increased  D  and 
S  values,  resulting  in  a  Kdef/K  ratio  of  less  than  1.  Average  values  of  D  and  S  were  measured  from  the 
random,  high  angle  grain  boundary  maps  of  each  material  lot  (except  NRL  HS)  and  are  listed  in  Table  3. 
The  maximum  angle  of  deflection  was  found  to  be  60  degrees  and  was  similar  for  all  materials,  which 
results  in  a  measured  Kdef/K  of  0.82,  0.82,  0.77,  and  0.81  for  Allvac,  TR1,  TR2,  and  TR3,  respectively.  As 
was  expected,  the  largest  grain  size  material  (TR2,  35.3  pm)  had  the  largest  amount  of  deflection  and  the 
lowest  Kde/K  ratio.  When  the  magnitude  of  decrease  in  K  is  evaluated  amongst  the  material  lots,  the  lot¬ 
to-lot  variation  in  KTH  is  too  large  to  be  solely  due  to  crack  deflection.  Specifically,  accounting  for 
deflection  results  in  an  approximate  20%  reduction  of  KTH  in  each  lot,  which  is  insignificant  given  the 
large  difference  in  measured  KTh  for  Allvac  and  TR1  at  -950  mVSCE  and  the  total  lack  of  1G-HEAC  for 
Allvac  and  NRL-HS  at  -850  mVSCE-  Therefore,  crack  deflection-induced  differences  in  K  are  not 
considered  to  be  the  dominant  source  of  lot-to-lot  variation  in  HEAC  susceptibility. 


4. 2. 2. 2.  Grain  Boundary  Connectivity 

Researchers  have  established  that  the  grain  boundary  character  (misorientation,  interface  energy, 
and  structure)  will  significantly  influence  the  1G  cracking  susceptibility'  in  a  wide  range  of  structural 
metals  by  influencing  the  cohesive  strength,  segregation/depletion  behavior,  and  slip 
behavior  [62,63,1 13].  Critically,  certain  low-energy  “special15  boundaries  are  more  resistant  to  the  failure 
mechanisms  pertinent  to  HEAC;  these  boundaries  include  grain  boundaries  of  low  misorientation  (15°  or 
less)  and  coincident  site  lattice  (CSL)  boundaries  of  1<E<29  [75,1 17].  Both  the  number  and  Sv  fraction 

of  low  energy  (CSL  and  low  angle  <15°)  GBs  are  similar  for  all  materials,  ranging  from  53-62%  with  no 
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discemable  trend  regarding  HEAC.  Percolation  theory  has  been  used  in  the  context  of  IG  stress  corrosion 
cracking  (SCC)  to  correlate  the  prevalence  of  susceptible  boundaries  to  the  probability  of  having  a  fully 
connected,  continuous  susceptible  crack  path  [62,1 18J.  Such  an  approach  assumes  a  binary  distribution 
of  either  fully  susceptible  (high  angle)  or  immune  GBs  (low  angle  or  CSL),  based  on  the  misorientation 
angle,  thereby  ignoring  any  gradients  in  susceptibility  due  to  other  microstructural  considerations. 
Regardless,  percolation  analysis  estimates  that  the  threshold  fraction  of  resistant  grains  necessary  to 
disrupt  a  connected  path  are  >1 1%  for  IG/TG  mix  and  >77%  for  fully  ductile  TG  (via  a  3D  analysis  of  a 
tetrakai decahedron  geometry  [118]),  35%  (unconstrained  2D  analysis  of  a  hexagonal  array  [119,120]), 
and  50%  (2D  analysis  of  both  a  microstructurally  constrained  hexagonal  array  [120]).  Based  on  these 
percolation  theory  thresholds  and  the  limited  range  of  low  energy  boundaries  for  all  materials  (53-62%), 
there  is  no  clear  evidence  to  support  the  variation  in  HEAC  resistance  being  due  to  differing  low  energy 
boundary  interruptions  of  the  continuous  susceptible  crack  path. 

While  the  percolation  theory  approach  provides  insightful  correlations  between  the  percentage  of 
susceptible  boundaries  and  the  probability  of  a  connected  path,  sole  reliance  on  the  percent  of  resistant 
boundaries  fails  to  account  for  potential  grain  size  influences.  For  example,  despite  similar  percentages 
(43  and  44%)  of  the  Sv  being  high  energy  in  Allvac  and  TR  2,  grain  size  differences  (13.8  and  35.3  pm, 
respectively)  result  in  the  total  Sv  of  high  energy  boundaries  of  the  two  lots  being  drastically  different: 
162  mm2/mm3  for  Allvac  and  70  mm2/mm3  for  TR2.  A  fractal  analysis  of  connectivity  provides  a  means 
to  understand  the  influence  of  these  grain  size  effects  on  the  availability  of  a  susceptible  crack  path  by 
quantifying  the  fractal  dimension  (DR)  of  the  path  of  maximum  connectivity  network  (MCN)  using  a  box¬ 
counting  method  [75]. 

The  boundary  misorientation  angles  are  characterized  via  OIM  analysis  of  each  material  (Figure 
9);  these  data  can  be  used  to  trace  the  maximum  connectivity  network  (MCN)  of  random,  susceptible, 
high  angle  boundaries.  This  process  is  illustrated  for  the  Allvac  material  in  Figure  13.  Figure  13a 
presents  the  OIM  map  (with  averaged  grain  orientations)  and  an  inverse  pole  figure.  From  these  data, 
traces  of  the  random  (black;  misorientation  angle  >5°)  and  CSL  (various  colors)  are  then  identified 
(Figure  13b).  The  high  angle  boundaries  (>15°)  are  extracted  (Figure  13c)  and  a  computer  program  is 
used  to  identify  the  MCN  (red  lines  in  Figure  13d)  for  each  material  lot.  For  this  MCN  analysis,  two  data 
points  were  considered  connected  if  they  were  adjacent  or  <2  pixels  apart  (thus  <1  pm),  so  that  short  non- 
indexed  paths  along  the  grain  boundary  do  not  affect  the  connectivity  analysis.  This  1  pm  threshold  is 
physically  justified  based  on  comparisons  to  x^t,  (proposed  to  be  ~1  pm  [40]);  specifically,  an  advancing 
crack  along  a  random,  high  angle  grain  boundary  would  not  likely  be  arrested  by  an  anomalous  segment 

of  “special”  boundary  shorter  than  xcrit.  The  MCN  is  then  overlaid  with  a  square  grid,  each  box  of  side 
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length  r|  (=12.5,  25,  or  50  pm).  For  each  box  size  (r|),  the  number  of  boxes  containing  part  of  the  MCN 
(N(rj))  is  quantified.  The  fractal  dimension  of  the  MCN  (DR)  is  the  metric  by  which  the  connectivity  is 
evaluated  [75]  and  is  defined  as: 


A 


(Eqn  6) 


Higher  values  of  DR  reflect  greater  connectivity  of  the  random,  susceptible  high  angle  grain 
boundaries  [75].  As  such,  this  approach  enables  a  quantitative  comparison  of  random,  high  angle  grain 


boundary  connectivity  for  each  Monel  K-500  material  lot.  To  provide  the  greatest  accuracy,  the  range  of 


q-values  are  on  the  scale  of  the  grain  sizes  observed  for  Monel  K-500,  consistent  with  prior  arguments 
and  application  of  this  fractal  analysis  approach  [75,121,122]. 

The  fractal  dimensions  are  reported  in  Table  3  and  show  no  clear  correlation  between  increased 
connectivity  (high  DR)  and  increased  HEAC  susceptibility,  thereby  suggesting  that  lot-to-lot  differences 
in  crack  path  connectivity  do  not  govern  the  observed  variability  in  HEAC  behavior.  For  example,  Allvac 
and  NRL-HS  shows  the  highest  and  third  highest  DR,  yet  neither  exhibit  IG-HEAC  at  -850  mVScE  (Figure 
8).  The  role  of  grain  size  on  the  crack  path  connectivity  is  illustrated  in  Figure  14  by  plotting  the  DR 
versus  both  the  Sv  fraction  of  low  energy  boundaries  and  the  Sv  value  of  high  angle  boundaries.  The  data 
do  not  show  a  correlation  between  increasing  fraction  of  low  angle  boundaries  and  DR  (thus  connectivity) 
for  materials  with  different  grain  sizes.  However,  there  is  a  reasonable  correlation  of  increasing 
connectivity  with  increasing  Sv  of  high  angle  boundaries.  This  demonstrates  that,  for  varying  grain  size 
materials,  the  fractional  values  of  resistant  low  energy  boundaries  (nor  high  energy  susceptible 
boundaries)  are  not  an  accurate  proxy  for  the  degree  of  crack  path  connectivity,  rather  the  true  value  of  Sv 
is  governing.  It  has  been  controversially  postulated  that  such  behavior  may  be  caused  by  2 3  boundaries 
playing  a  marginal  role  in  disrupting  the  network  of  susceptible  boundaries  due  to  their  formation  in  the 
grain  interior,  which  would  fail  to  disrupt  a  susceptible  path  [120]. 

4.3  Hydrogen  Interactions 

Previous  work  found  different  diffusible  hydrogen  concentrations  as  a  function  of  applied 
potential  in  two  Monel  K-500  lots,  thereby  suggesting  that  the  relationship  between  the  diffusible 
hydrogen  concentration  (thus  the  stress-enhanced  concentration,  CH-a«»  via  Eqn  1)  and  applied  potential 
may  vary  amongst  the  material  lots  [18].  This  variation  is  confirmed  by  the  differences  in  measured  Cn-Diff 
found  for  three  lots  (Allvac,  NRL-HS,  and  TR2)  at  hydrogen  overpotentials  ranging  from  -0.85  to  -1.1 
VSce  (Figure  10).  A  lower  value  of  CH-Diff  leads  to  a  concomitant  reduction  of  CH-o  for  a  given  hydrogen 
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overpotential.  Moreover,  the  hydrogen  coverage  at  grain  boundaries  in  the  fracture  process  zone  for  a 
heat-invariant  grain  boundary  trap  binding  energy  would  be  reduced  by  the  lower  Cn,Diff;  though  it  should 
be  noted  that  grain  boundary  trap  states  are  difficult  to  distinguish  from  the  y*  trap  state,  which  is  the 
subject  of  future  work  [52].  A  decrease  in  CH-Difr  would  increase  HEAC  resistance  (as  demonstrated  in 
Figure  10)  by  slowing  accumulation  of  the  H  concentration  required  to  nucleate  crack  growth,  CHa-crit- 
This  change  in  H-material  interaction  would  interact  with  the  previously  discussed  lot-to-lot  differences 
(e.g.  S  segregation,  gettering  element  content,  and/or  high  grain  boundary  surface  to  volume  ratio  altering 
the  kio  and  CHa-crit)  to  influence  the  material  specific  HEAC  behavior.  These  speculations  could  be 
reinforced  by  better  assessment  of  elemental  sulfur  and  gettering  levels  and  trap  binding  energies  for 
high-  and  low-angle  boundaries. 

4.4  Modeling  KTH  for  Lot-to-Lot  Variations  in  Yield  Strength  and  Composition 

The  results  of  the  sensitivity  analysis  shown  in  Figure  lla-c  demonstrates  that  significant 
variations  in  KTH  can  arise  from  realistic  lot-to-lot  changes  in  yield  strength,  grain  boundary  toughness, 
and  hydrogen  uptake  behavior.  Furthermore,  the  discussion  presented  herein  has  suggested  a  potentially 
important  role  for  impurity  segregation  to  grain  boundaries  in  determining  a  material  lot’s  HEAC 
susceptibility,  specifically  with  regards  to  the  deleterious  effect  of  sulfur  and  the  beneficial  effect  of 
boron.  In  the  following  section,  the  KTh  model  (Eqn.  2)  is  leveraged  to  predict  the  liEAC  behavior  of  the 
five  material  lots;  this  effort  will  use  previously  developed  modeling  parameters  [18]  along  with  material 
lot  specific  inputs  as  a  basis  to  propose  model  improvements  that  incorporate  the  effects  of  impurity 
segregation  to  GBs. 

Predicted  KTH  versus  Eapp  results  for  Allvac,  TR1,  TR2,  TR3,  and  NRL-HS  are  shown  in  Figure 
15  using  previously  validated  parameters  for  the  Allvac  material  (a”  =  0.0002  MPaVm,  (3’  =  0.20 
(MPaVm)'1,  kio  =  0.88  MPaVm,  and  a  =  4.078  MPaVm/(atm  frac  H)  [18,54])  and  each  material’s 
respective  yield  strength  (Table  2).  As  expected,  the  variation  in  yield  strength  amongst  the  material  lots 
(and  its  effect  on  the  local  hydrostatic  stress)  give  rise  to  drastic  differences  in  calculated  threshold  values 
as  the  applied  potential  is  increased.  However,  as  shown  by  the  experimental  results  plotted  alongside  the 
predicted  curves,  the  previously  validated  modeling  parameters  are  not  able  to  fully  capture  the  effect  of 
microstructure  changes  on  HEAC  susceptibility.  In  particular,  the  NRL-HS  is  predicted  to  be  the  most 
susceptible  of  the  five  lots  (blue  dashed  line),  but  was  found  to  be  immune  to  HEAC  when  experimentally 
tested  at  -850  mVSCE,  as  shown  by  the  black  arrow  extending  upward.  This  deviation  from  experimental 
results  could  be  rooted  in  the  influence  of  sulfur  and  boron  on  the  grain  boundary  toughness,  kiG; 
especially  with  regards  to  NRL-HS  which  has  the  highest  bulk  sulfur  content.  To  further  explore  the 
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effect  of  impurity  segregation,  we  make  two  assumptions:  (1)  any  effect  of  impurity  elements  on  kiG 
occurs  prior  to  H  exposure  and  (2)  sulfur  and  boron  are  the  only  two  elements  influencing  the  GB 
toughness.  This  leads  to  the  following  modification  of  Eqn.  2: 


^th 


_ !  ((kip  f3Cs+yCB )  flCHq) 

P*  aM^Ys 


(Eqn  7) 


Where  (3  (MPaVm/(atm  frac  S))  and  y  (MPaVm/(atm  frac  B))  are  the  weighting  factors  that  define  how 
effectively  sulfur  and  boron  alter  kIG,  Cs  and  CB  are  the  concentrations  of  sulfur  and  boron  at  the  grain 
boundary,  and  all  other  variables  are  defined  as  in  Eqn  2.  To  predict  the  concentrations  of  boron  and 
sulfur  at  an  arbitrary  grain  boundary,  the  Auger  analysis  of  Ladna  and  Bimbaum  were  used  to  estimate 
the  amplification  of  the  bulk  concentration  [108].  These  amplifications  (48,000x  and  ~l,270x  for  S  and 
B,  respectively)  were  linearly  scaled  to  reflect  differences  in  grain  size  between  the  original  and  current 
study,  resulting  in  final  amplifications  of  approximately  180x  and  3x  in  NRL-HS  (for  S  and  B, 
respectively).  Cs  and  CB  were  calculated  for  each  material  lot,  along  with  the  previously  calculated  CHa, 
and  then  a,  |3,  and  y  were  fit  to  minimize  the  square  of  the  difference  between  the  calculated  KTH  and  the 
experimentally  measured  KTHfor  each  lot  at  -950  mVScE,  where  lot-to-lot  variations  in  metallurgy  should 
have  minimal  effect.  From  this  analysis,  a  =  4.167  (MPaVm/(atm  frac  H)),  P  =  2.753  (MPaVm/(atm  frac 
S)),  and  y  =  4. 12x1  O'4  (MPaVm/(atm  frac  B)),  which  suggests  that  boron  has  minimal  effect  on  grain 
boundary  toughness. 

Calculated  KTh  versus  Eapp  using  the  newly  fit  parameters  are  shown  in  Figure  16  for  Allvac, 
TR1,  TR2,  TR3,  and  NRL-HS,  along  with  experimental  results  for  -950  and  -850  mVScE-  Despite  a 
noticeable  shift  in  the  results,  especially  at  less  cathodic  potentials,  a  significant  difference  in  theoretical 
and  experimental  results  still  exists.  Critically,  several  deficiencies  in  the  fidelity  of  the  data  inputs  need 
to  be  addressed  before  the  efficacy  of  this  modeling  approach  can  be  determined,  specifically:  (1) 
inaccurate  estimates  of  sulfur  and  boron  concentrations  at  the  grain  boundary,  (2)  lack  of  specific  data  on 
the  effect  of  gettering  elements  on  sequestering  sulfur  in  Monel  K-500,  and  (3)  lack  of  specific  hydrogen 
uptake  relationships  for  non-Allvac  lots.Physically,  the  proposed  alterations  to  Eqn  2  are  reasonable,  as 
sulfur  and  boron  have  been  shown,  both  experimentally  and  theoretically,  to  impact  grain  boundary 
cohesion  [79,108,109,123].  Improvements  to  the  modeling  would  therefore  be  based  on  gathering  data 
specific  to  the  alloy  of  interest  so  as  to  ensure  that  unique  factors  are  accurately  captured.  Specifically, 
Auger  analysis  of  GBs  in  each  lot  as  a  function  of  applied  heat  treatment  would  provide  an  accurate  GB 
amplification  factor  for  boron  and  sulfur,  while  also  addressing  questions  regarding  the  efficiency  of 

gettering  elements.  Further  still,  understanding  the  specific  hydrogen  uptake  relationship  for  each  material 
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is  critical  to  accurately  estimating  the  stress-enhanced  hydrogen  concentration  in  the  FPZ.  As 
demonstrated  in  Figure  11a,  realistic  changes  in  uptake  behavior  can  have  drastic  influences  on  the 
calculated  threshold  values. 

Finally,  a  variation  in  modeling  parameter  dependencies,  especially  regarding  the  dependence  on 
yield  strength,  has  been  found  for  different  materials  and  should  be  evaluated  in  detail  for  the  Ni-based 
alloys.  Specifically,  a  less  severe  dependence  of  KTh  on  yield  strength  was  found  for  Mo  (exp(a_1 2))  as 
compared  to  the  original  Fe-3%Si  solution  (exp(a’1))  [73].  However,  previous  efforts  have  shown  that 
values  predicted  using  Eqn  2  were  in  excellent  agreement  with  experimental  results  [18,54],  suggesting 
that  other  influences  may  be  dominating. 

Conclusions 

The  influence  of  microstructure  variation  on  hydrogen  environment-assisted  cracking  of  Monel 
K-500  are  evaluated  using  five  nominally  peak-aged  lots  of  material  (all  of  which  satisfy  the  QQ-N-286 
specification)  tested  under  slow-rising  stress  intensity  loading  while  immersed  in  NaCl  solution  under 
cathodic  polarizations.  Quantitative  HEAC  parameters  were  established  using  fracture  mechanics  testing 
and  were  coupled  with  a  detailed  characterization  of  the  microstructure  (e.g.  composition,  grain  character 
analysis,  and  H-metal  interaction  parameters.)  The  results  of  these  analyses  are  coupled  with 
micromechanical  models  to  quantitatively  evaluate  how  observed  changes  in  the  microstructure  relate  to 
the  resulting  HEAC  behavior.  The  following  conclusions  are  established. 

•  Limited  lot-to-lot  variation  in  HEAC  susceptibility  was  observed  in  five  peak-aged  Monel  K- 
500  lots  when  immersed  in  0.6M  NaCl  under  an  applied  potential  of -950  mVSCE  during 
slow-rising  K  testing,  suggesting  minimal  influence  of  metallurgical  variation  in  severe  H- 
producing  environments. 

•  Extensive  lot-to-lot  variation  in  both  HEAC  metrics  (KTH  and  da/dtn)  and  fractography  was 
observed  amongst  the  material  lots  when  immersed  in  0.6  M  NaCl  under  an  applied  potential 
of -850  mVSCE  during  slow-rising  K  testing. 

•  The  influence  of  metallurgical  features,  including  grain  size,  impurity  segregation,  grain 
boundary  character  and  connectivity,  gettering  element  composition,  and  crack  path 
tortuosity  were  evaluated.  Of  these,  yield  strength  and  impurity  segregation  are  considered  to 
be  the  most  influential  in  determining  HEAC  resistance. 
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o  The  degree  of  impurity  segregation  is  critically  influenced  by  the  grain  size,  inclusion 
of  gettering  elements,  and  the  bulk  concentration  of  both  beneficial  (e.g.  B)  and 
deleterious  (e.g.  S)  elements. 

•  Sensitivity  analysis  of  the  decohesion-based,  micromechanical  models  suggests  that 
physically  realistic  variations  in  hydrogen  uptake  behavior,  yield  strength,  and  intrinsic  grain 
boundary  toughness  give  rise  to  large  changes  in  HEAC  susceptibility. 

•  Predicted  values  of  KTH  did  not  align  with  experimental  results  when  previously  validated 
modeling  parameters  were  employed,  suggesting  that  the  current  models  need  to  be  .extended 
to  more  accurately  capture  the  effect  of  lot-to-lot  variations  in  microstructure. 

•  A  modification  to  KTh  model  was  proposed  to  incorporate  the  effect  of  impurity  segregants, 
resulting  in  new  modeling  parameters.  Correlation  between  experimental  and  predicted 
values  were  compromised  by  uncertain  inputs  regarding  the  degree  of  grain  boundary 
segregation.  Ongoing  efforts  will  establish  these  parameters  to  provide  a  more  rigorous 
evaluation. 
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Table  1  -  Chemical  Composition  of  Monel  K-500  Lots  (wt.  %,  except  S,  P,  Sn,  Pb,  Mg,  Zr,  Hf,  B  are  wppm) 


Lot 

Ni 

Cu 

Al 

Fe 

Mn 

Si 

Ti 

C 

Co 

S 

P 

Sn 

Pb 

Mg 

Zr 

Hf 

B 

Allvac 

66  12 

28  57 

2  89 

0.80 

0.81 

0.08 

0.45 

0.17 

0.01 

1.60 

92 

2.4 

2.1 

39 

370 

26 

3.8 

NRL-HS 

63.44 

30  74 

3.20 

0.91 

0.85 

0.10 

0.57 

0.14 

0.02 

17  00 

40 

2.2 

3.5 

40 

650 

0.8 

2.0 

TR1 

65.55 

29  18 

2  86 

0  93 

0  67 

0  04 

0.47 

0.22 

0  01 

3.10 

55 

1  2 

2  4 

170 

230 

4.6 

0.3 

TR2 

64.66 

30.15 

2  73 

0.69 

0.73 

0.09 

0.45 

0.20 

0.04 

11.00 

71 

6.9 

2.5 

130 

330 

5.6 

1.5 

TR3 

64  17 

30.34 

3.21 

0.71 

0.81 

009 

0.45 

0  13 

004 

12  00 

74 

6.8 

24 

130 

320 

5.4 

1,5 

QQ-N- 

286 

>  63.0 

27-33 

2.3- 

3.15 

< 

2.0 

< 

1.5 

< 

0.50 

0.35- 

0.85 

< 

0.18 

< 

0.25 

<60 

< 

200 

<60 

<60 

N/A 

N/A 

N/A 

N/A 

Table  2  -  Mechanical  Properties  of  Monel  K-500  Lots 


Lot 

ctvs  (MPa) 

E  (GPa) 

n 

a 

Allvac 

786 

180 

20 

0.39 

NRL-HS 

910 

191 

18 

0.405 

TR1 

898 

183 

18 

0.39 

TR2 

795 

202 

20 

0  39 

TR3 

792 

194 

20 

0.39 

Table  3  -  Grain  Boundary  Size  and  Characteristics  of  Monel  K-500  Lots 


Lot 

d  (Mm) 

Sv 

(mm2/mm3) 

Low  Angle  Boundary 
Fraction  (%) 

Fractal 

Dimension,  Dr 

D 

(pm) 

S 

(pm) 

Allvac 

13.8 

341 

57 

1.59 

10.5 

10.8 

NRL-HS 

11.2 

407 

62 

1.49 

N/A 

N/A 

TR1 

25  8 

252 

60 

1.53 

17.6 

18.3 

TR2 

35.3 

149 

56 

1.38 

23.2 

14.0 

TR3 

31.8 

160 

53 

1.35 

24.1 

23.4 
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Figure  1  -  Crack  growth  rate  versus  total  K  for  five  lots  of  peak-aged  Monel  K-500  stressed  under  two 
different  dK/dt  protocols  in  a  dry  nitrogen  (RH  <  5%)  environment.  Closed  symbols  represents  specimens 
tested  at  dK/dt  =  1  MPaVm/hr,  while  open  symbols  were  tested  at  dK/dt  =  0.33  MPaVm/hr.  Expected 
false  da/dt  produced  by  crack  tip  plastic  deformation  for  each  dK/dt  are  represented  by  the  solid  (1 
MPaVm/hr)  and  dashed  (0.33  MPaVm/hr)  lines. 
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Figure  2  -  Fractographs  illustrating  the  (a)  dimpled,  ductile  fracture  morphology  observed  in  the  center 
(b)  of  inert  environment-tested  specimens.  Additional  images  of  (c)  Allvac  and  (d)  TR2  show  the  carbide 
precipitates  (red  boxes)  found  at  the  base  of  many  micrcvcic  features. 
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Figure  3  -  Crack  growth  rate  versus  total  K  for  five  lots  of  peak-aged  Monel  K-500  stressed  in  NaCI 
solution  at  a  constant  applied  potential  of  -950mVSCE-  The  dashed  line  represents  the  expected  false 
da/dt  produced  by  crack  tip  plastic  deformation  for  dK/dt  =  0.33  MPaV m/hr. 
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Figure  4  -  F'actcgraphy  of  the  intergranular  (16)  fracture  morphology  observed  for  (a)  Allvac,  (b)  TR1,  (c) 
TR2,  (d)  TF3  and  (e-f)  NRL  HS  when  tested  in  D.6M  NaCI  solution  with  an  app  ied  potential  0"  -950mVScE- 
In  additior  to  the  IG  morphology,  RL  HS  (e-f)  also  exhibited  some  transgranu  ar  features. 
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Figure  5  -  rractographs  of  facet  slip  ".races  for  (a)  Allvac,  (b)  TR1,  (c)  TR2,  (d)  TR3,  and  (e-f)  NRL  HS  when 
tested  in  0.6M  NaCI  solution  with  an  applied  potential  of  -S50mVSCE. 
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Figure  6  -  Crack  growth  rate  versus  total  K  for  five  lots  of  peak-aged  Monel  K-500  stressed  in  NaCI 
solution  at  a  constant  applied  potential  of  -850mVSCE.  The  dashed  line  represents  the  expected  false 
da/dt  produced  by  crack  tip  plastic  deformation  for  dK/dt  =  0.33  MPaVm/hr. 
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Figure  7  -  Fractcgraphy  of  samp  es  tested  at  an  applied  potential  of  -850  mVSCE  in  0.6M  NaCI:  (a) 
overview  of  TR2  fracture  surface  and  the  fracture  morphology  observed  at  the  crack  length  where  K  =  40 
MPa'/m  for  (b)  TR1,  (c)  TF.2.  and  (d)  TR3. 
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Figure  8  -  Comparison  of  (a)  crack  growth  kinetics  data  for  NRL  HS  tested  in  dry  N2  3nd  C.6M  NaCI  at  -850 
mVscE  and  fractographs  of  the  (b)  dimpled,  ductile  fracture  morphology  observed  in  the  center  (c)  of  the 
0.6M  NaCI  at  -850  mVScE  specimen. 
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Figure  9  -  Orientation  maps  for  (a)  Allvac,  (b)  TR1,  (c)  TR2,  (d)  TR3,  and  (e)  NRL  HS  collected  using 
electron  backscatter  diffraction.  Grain  boundaries  (black  lines)  were  defined  by  misorientations  greater 
than  15°. 
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Figure  10  -  Effect  of  hydrogen  overpotential  on  total  and  diffusible  hydrogen  concentrations  for  the  Allvac,  TR2, 
and  NRL-HS.  Error  bars  represent  95%  confidence  interval. 
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Figure  11  -  Sensitivity  analysis  conducted  on  the  modeling  to  determine  the  effect  of  (a)  a  ±  250% 
change  in  diffusible  hydrogen  concentration  relative  to  the  baseline  value  found  for  Allvac,  (b)  a  ±  25% 
change  in  the  Griffith  toughness  of  the  grain  boundary  (ktG)  relative  to  the  baseline  value  of  0.88 
MPaVm,  and  (c)  a  ±  15%  change  in  the  yield  strength  from  the  Allvac  yield  strength  (786  MPa)  on  the 
calculated  threshold  stress  intensity. 
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Figure  12  -  Sensitivity  analysis  conducted  on  the  modeling  to  determine  the  effect  of  (a)  a  ±  250% 
change  in  diffusible  hydrogen  concentration  relative  to  the  baseline  value  found  for  Allvac  and  (b)  a  ± 
15%  change  in  the  yield  strength  from  the  Allvac  yield  strength  (786  MPa)  on  the  calculated  H-diffusion 
limited  crack  growth  rate,  da/dtn. 
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Figure  13  -The  process  for  establishing  the  fractal  dimension  included:  (a)  obtaining  an  orientation  map 
via  electron  backscatter  diffraction,  (b)  identification  of  high  (black  lines)  angle  grain  boundaries  and  low 
(red  lines)  angle  grain  boundaries,  (c)  isolating  high  angle  grain  boundaries  by  removing  all  low  angle 
grain  boundaries,  (d)  measuring  the  maximum  connected  network  (highlighted  red  lines)  and  overlaying 
a  square  grid  with  a  known  side  length. 
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Figure  14  -  Fraction  of  low  energy  boundaries  (open  symbols)  and  the  Sv  of  high  energy  boundaries 
(closed  symbols)  versus  the  measured  fractal  dimension  for  each  lot  of  Monel  K-500. 
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Figure  15  -  Plot  of  the  calculated  threshold  stress  intensity  versus  applied  potential  curves  for  each  lot  of 
Monel  K-500  using  previously  validated  parameters  with  experimental  data  points  included  for  :950  and 
-850  mVSCE.  Black  arrows  extending  upward  represent  experimental  observations  of  HEAC  immunity  at 
the  tested  potential. 
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Figure  16  -  Plot  of  the  calculated  threshold  stress  intensity  versus  applied  potential  curves  for  each  lot  of 
Monel  K-500  using  new  parameters  with  experimental  data  points  included  for  -950  and  -850  mVSCE. 

New  parameters  include  the  effects  of  grain  boundary  impurity  segregation  of  sulfur  (P)  and  boron  (y). 
Black  arrows  extending  upward  represent  experimental  observations  of  HEAC  immunity  at  the  tested 
potential. 
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Chapter  2 


Effects  of  Heat  to  Heat  Variations  in  Hydrogen  Interactions  on  the  Hydrogen  Embrittlement 

Susceptibility  of  Monel  K-500 

B.C.  Rincon  Troconis,  H.  Ha,  Z.  Harris,  J.  Burns  and  J.R.  Scully 
Center  for  Electrochemical  Science  and  Engineering 
Department  of  Materials  Science  and  Engineering 
University  of  Virginia 
Charlottesville,  VA  22904 


Abstract 

The  influence  of  heat  to  heat  variation  on  hydrogen  environment-assisted  cracking  (HEAC)  of 
aged  hardened  Monel  K-500  (UNS  N05500,  63  %Ni,  33  %Cu,  3.15  %A1,  wt.  %)  was  evaluated.  Four 
peak-aged  lots  of  material  laboratory  were  tested  under  slow  strain  rate  tensile  testing  (SSRT)  while 
immersed  in  NaCl  solution  under  various  cathodic  polarizations.  All  materials  met  a  US  federal 
specification  based  on  mechanical  properties  and  composition.  One  heat  was  harvested  from  seawater 
deployment  where  HEAC  was  observed  after  exposure  at  unknown  levels  of  cathodic  protection.  Another 
heat  cracked  at  less  severe  cathodic  potentials  in  field  ocean  exposure  when  pre-cracked  despite  lower 
strength.  Upon  removal,  hydrogen  was  detected  by  thermal  desorption  in  both  cases.  In  addition  to  SSRT, 
selected  microstructural,  compositional,  as  well  as  hydrogen-material  interactions  were  investigated  and 
compared.  These  included  assessments  of  diffusible  and  total  hydrogen  concentrations,  trap  binding 
energies,  diffusion  rate,  grain  size,  hardness,  and  compositional  analysis  including  sulfur,  key  trace  sulfur 
gettering  elements  and  other  beneficial  elements  which  alter  sulfur  solubility.  Moreover,  grain  boundary 
chemistry  was  assessed  after  aging  and  hydrogen  fracture  by  Auger  Electron  Spectroscopy.  Sulfur 
enrichment  was  observed.  All  materials  were  susceptible  to  hydrogen  embrittlement  with  differences 
between  field  susceptible  heats  and  the  others.  The  heats  identified  in  seawater  and  field  testing  exhibited 
more  facile  intergranular  cracking,  absorbed  more  hydrogen  relative  to  other  heats  as  a  function  of 
cathodic  polarization,  and  contained  unfavorable  concentrations  of  trace  detrimental  grain  boundary 
segregating  elements  relative  to  the  beneficial  gettering  elements  in  combination  with  a  low  density  of 
grain  boundary  trap  or  segregation  sites  per  unit  volume.  Greater  susceptibility  is  traced  to  enhance 
hydrogen  uptake  and  speculated  grain  boundary  sulfur  enrichment  based  on  alloy  trace  chemistry  coupled 
with  grain  boundary  trapping  of  hydrogen  supplied  from  charging  during  cathodic  polarization.  Material 
hardness  within  the  limits  from  Rockwell  C  28  to  36  was  not  found  to  be  the  dominant  factor. 
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Introduction 


Monel  K-500  is  a  precipitation  hardenable  nickel  based  alloy  with  composition  63%Ni,  33  %Cu, 
3.15  %A1,  1.5  %Mn,  2  %Fe,  0.85  %Ti  (in  wt%)  [1].  This  alloy  is  utilized  for  different  applications  such 
as  marine  service,  chemical  processing,  pulp  industry,  oil  and  gas  production  and  electronic  components. 
The  matrix  is  characterized  by  a  Ni-Cu  face  centered  cubic  structure  [2].  The  peak  aged  materials  present 
an  uniform  distribution  of  spherical  Ni3(Al,X)  precipitates  known  as  y\  where  X  can  be  Cu,  Mn,  Ti  or  Si. 
This  alloy  commonly  also  contains  incoherent  carbides  in  the  form  of  TiC  [2]  that  suffer  little  coarsening 
and  shape  change  during  aging.  But,  when  alloyed  for  long  period  of  times  other  carbides  are  also  present 
in  the  chemical  form  of  M23C6  (M:  Ni,  Mn  or  Fe).  Another  microstructure  feature  that  can  be  present  in 
this  alloy  are  the  annealing  twins  [2].  These  are  present  in  solution  heat  treated  (SHT)  aged,  SHT  water 
quenched  (WQ)  and  SHT  air  cooled  (AC)  materials.  All  these  microstructural  features  can  interact  with 
H.  Despite  of  the  alloy  immunity  in  some  environments,  structural  failures  related  to  hydrogen 
embrittlement  (HE)  have  been  reported  [3-6]. 

The  reported  failures  took  place  under  sacrificial  cathodic  polarization  in  seawater.  Reduction  in 
tensile  ductility  and  notch  strength  was  observed  for  SHT  aged  Monel  K-500  alloy  coupled  to  zinc 
producing  potentials  from  -1  to  -1.1  Vsce-13  [3,  5].  Intergranular  failures  of  bolts  subjected  to  cathodic 
polarization  have  also  been  reported  in  age  hardened  Monel  K-500  alloys  [4,  6].  These  failures  were 
attributed  to  the  hardening  of  the  material  through  age  hardening  after  threading,  achieving  a  hardness 
Rockwell  C  (HRC)  of  approximately  39.  Therefore,  a  modification  of  the  heat  treatment  [4],  involving 
annealing  after  threading  followed  by  age  hardening,  was  proposed  to  maintain  the  hardness  bellow  HRC 
35.  In  despite  of  this  processing  modification,  bolt  failures  continued  to  be  reported  at  lower  hardness 
levels  [6].  Threaded  bolts  annealed  at  a  temperature  of  980-1050  °C  during,  water  quenched  (WQ)  and 
precipitation  hardened  at  500-600  °C  during  16  hours  with  a  final  hardness  of  25  (Below  limit 
recommended  in  US  federal  specification  [3])  failed  [6].  Sporadic  failures  occurred  after  1  year  exposure 
under  sacrificial  cathodic  protection  using  A1  anodes  under  a  load  representing  60%  of  the  yield  strength 
of  the  material  [6].  Slow  strain  rate  testing  (SSRT)  experiments  performed  in  ASTM  artificial  seawater  on 
pre-charged  specimens  potentiostatically  or  using  A1  and  Fe  anodes,  (representing  high  and  low  hydrogen 
overpotential)  resulted  in  intergranular  cracking . 

HE  failures  have  been  indicated  to  result  from  grain  boundary  segregation  particularly  in 
combination  with  hydrogen  [7,  8].  Elements  present  in  the  bulk  alloys  in  the  order  of  ppm  (trace 

13  All  potentials  are  hereafter  referenced  to  Saturated  Calomel  Electrode 
(+0.241  V  vs .  a  standard  hydrogen  electrode) 
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elements)  can  segregate  to  the  grain  boundary,  producing  highly  localized  enriched  volume  of  material 
with  concentrations  with  orders  of  magnitude  higher  than  that  of  the  bulk  concentration.  That  is  the  case 
of  elements  such  as  S,  Pb,  P,  Zr,  Hf,  Ca  and  Mg  [9-12].  The  segregation  strength  will  depend  on  different 
parameters  such  as  the  specific  element,  bulk  solubility,  segregation  energy,  grain  boundary  nature  and 
the  presence  of  other  elements  that  can  interact  with  the  segregating  element,  decreasing  segregation  [13, 
14].  Hydrogen  is  thought  to  segregate  preferentially  to  the  grain  boundaries  facilitating  grain  boundary 
de-cohesion,  leading  to  HE  [15].  Nickel  hydriding  has  been  detected  [16-18].  But,  is  not  believed  to  be 
the  main  factor  controlling  failure  in  sea  water  under  typical  cathodic  protection  conditions  although,  it  is 
observed  at  extreme  hydrogen  uptake  conditions.  The  process  of  hydrogen  embrittlement  proposed  is  in 
good  agreement  with  experiments  and  involves  crack  tip  hydrogen  production,  uptake,  and  diffusion  into 
the  fracture  process  zone  and  at  grain  boundaries  where  fracture  occurs  [15,  19-21].  This  process  can  be 
enhanced  by  S  presence  in  the  material,  which  is  an  element  known  to  be  deleterious  for  Ni-based  alloys 
[8,  22,  23].  Traces  of  this  element  can  result  in  grain  boundary  concentrations  that  are  2000  times  the  bulk 
concentration  [7,  23,  24].  It  has  been  postulated  that  the  electronegativity  of  this  element  interacts  with 
the  bonding  strength  in  the  grain  boundary,  resulting  in  grain  boundary  de-cohesion  [22].  This  hypothesis 
has  been  rationalized  in  terms  of  the  local  spin  density'  atomic  cluster  [25]  and  quantum-mechanics  cluster 
calculations.  Another  element  known  to  segregate  to  grain  boundary  is  P  [23].  The  literature  offers 
information  regarding  the  inert  and  even  beneficial  effect  of  this  element  [10,  26,  27].  Sb  is  another 
element  that  has  been  found  to  be  deleterious  for  alloy  performance.  Sb  segregates  to  the  grain  boundary 
and  lowers  the  threshold  stress  intensity  [28].  H  and  S  have  been  found  to  act  synergistically  towards  HE 
[1,  29-32].  A  relationship  taking  into  consideration  the  average  grain  boundary  fracture  surface  coverage 
of  segregating  elements  (Cx)  and  their  fractional  relative  potency  with  respect  to  S  to  promote  hy  drogen 
cracking  has  been  defined  empirically  from  straining  electrode  data  [32]: 

Cesq  =  Cs  +  0.07  Cp  +  0.6  Csb  Eq.  1 

Also,  another  equation  has  been  formulated,  in  this  case  to  explain  the  effect  of  impurity 
segregations  on  the  fracture  stress  of  a  notched  specimen  [33], 

/t  —  vacuum  nfrx  pry  Fn  9 

Gfrac  -  afrac  ~  aLH  “P4  nCi-z 

where,  x  and  y  are  ty  pical  fractions,  cr^?ac^um  is  the  fracture  stress  in  vacuum  and  a  and  p  describe  a 
weighting  factor  of  the  deleterious  effect  of  impurities  on  fracture  stress. 
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S  is  known  to  interact  with  hydrogen  and  enhance  HE.  However,  there  are  certain  elements, 
called  getter  elements,  that  can  be  intentionally  added  to  the  alloy  to  react  with  S  and  avoid  its  segregation 
to  the  grain  boundary.  These  elements  can  form  sulfo-carbides  (Zr,  Hf)  [12,  34]  and  sulfides  (Mg,  Ca,  Zr) 
[23].  Yamaguchi  developed  a  model  to  evaluate  the  S  gettering  power  of  these  elements  based  an 
empirical  analysis  considering  experiments  related  to  hot  workability  and  chemical  characterization  of 
grain  boundary  through  Auger  electron  spectroscopy  [35].  A  modified  expression  is  available  considering 
the  addition  of  Hf  effect  with  a  similar  weight  as  Zr  (similar  solubility  product)  and  neglecting  the 
gettering  power  of  Ca  as  it  is  considered  to  be  present  mainly  as  stable  CaO  [36]. 


CGet=  0.3Mg  +  0.1Zr +  0.1Hf  Eq.  3 

Hardness  has  been  explored  as  a  parameter  that  could  affect  the  HE  of  Ni-based  alloys  [37,  38]. 
The  hardness  of  an  Inconel  material  was  reduced  using  laser  driven  surface  modification.  The  laser 
treatment  modified  the  surface  layer  of  the  material,  up  to  1  mm  in  thickness,  dissolving  y’  precipitates 
and  enhancing  the  material  ductility  of  the  alloy  under  H2  environment  when  compared  to  a  non-treated 
material  exposed  to  the  same  conditions  [37].  Hardness  is  a  key  parameter  in  decohesion  models  applied 
to  Monel  K-500.  It  can  be  shown  that  hardness  alone  can  produce  a  drop  in  Kth  and  rise  in  stage  II  crack 
growth  with  hydrogen;  but  significantly,  has  a  strong  effect  on  the  cathodic  potential  threshold  [15,  36]. 

A  relationship  between  the  critical  stress  for  cracking,  a* ,  and  material  dependent  variables  such 
as  the  uniaxial  tensile  yield  strength,  oy,  in  the  presence  of  gaseous  hydrogen  for  steel  has  been 
established  [38], 


a'  =  ex*  -  A(PhJ  /2exp 


V  h  (2.6<7h  —  1.3  <7y) 


3  RT 


Eq.  4 


where  a0*,  A,  VH  are  materials  variables  and  represent  the  material-dependent  constants  (o0*  and  A)  and 
the  partial  molar  volume  of  hydrogen  inside  the  metal,  PH2  is  the  hydrogen  pressure,  ou  is  the  principle 
local  stress  perpendicular  to  the  crack,  R  is  the  gas  constant  and  T  is  temperature.  Modelling  of  the 
threshold  stress  intensity  utilizing  Eq.  4  shows  that  for  the  same  PH2  an  increase  on  ay  represents  a 
decrease  in  the  threshold  stress.  However,  it  is  important  to  point  out  that  oy  cannot  be  varied 
independently  and  that  the  values  o0*  and  A  cannot  be  measured,  which  makes  impossible  the  comparison 
of  Eq.  4  with  real  data.  Another  study  [39]  based  on  the  effect  of  aging  temperature  and  therefore 
hardness  on  hydrogen  permeation  parameters  elucidates  the  increase  in  absorbed  hydrogen  with  hardness 
for  18Ni(250)  steel,  which  can  result  in  a  decrease  in  critical  stress  intensity  [40]  and  the  threshold  stress 
intensity  for  cracking  [41].  In  addition,  the  diffusivity  was  found  [39]  to  decrease  with  hardness  and  then 
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achieve  steady  state  for  over  aged  materials,  which  can  result  in  slow  crack  growth  rates  for  diffusion 
controlled  cracking  [40].  In  a  secondary  hardened  martensitic  steel  greater  hardness  brought  about  by 
M2C  carbides  provided  traps  sites  for  hydrogen  which  could  repartition  to  crack  tips  [42]. 

Hydrogen  uptake,  diffusivity  and  trap  binding  energy  were  determined  [43]  for  single  phase  solid 
solution,  aged,  and  cold  worked  then  aged  Monel  K-500  alloy.  The  total  hydrogen  concentration,  CH,totin 
alkaline  3.5  %  NaCl  was  found  to  be  in  the  range  0-132  ppm,  while  the  diffusible  hydrogen  concentration, 
Cn.diff,  was  lower  and  in  the  range  0-36  ppm,  approximately  [43].  Both  hydrogen  concentrations  were 
found  to  be  dependent  on  hydrogen  overpotential.  The  diffusible  hydrogen  concentration  was  found  to  be 
dominated  by  the  hydrogen  in  the  lattice  and  to  a  lesser  extent  some  hydrogen  in  reversible  traps, 
excluding  the  y’  precipitates,  while  the  total  hydrogen  concentration  included  the  diffusible  hydrogen,  the 
hydrogen  present  in  the  reversible  y’  precipitate  traps  and  the  hydrogen  in  irreversible  trap  sites  [43].  The 
reversible  traps  could  involve  microstructural  features  such  as  copper  solute  in  the  solid  solution  matrix, 
dislocations,  twin  boundaries,  low  energy  grain  boundaries  and  some  sites  on  y’  precipitates  [44,  45]. 
Grain  boundaries  and  dislocations  have  been  reported  as  hydrogen  trap  sites  [46];  however  they  have  a 
low  site  density  in  the  order  of  1017  sites/cm\  except  in  ultra-fine  grain  size  microstructures  [47,  48].  The 
room  temperature  hydrogen  diffusion  coefficient  was  found  to  range  from  0.9  to  3.9  x  10‘14  m2/s  [43].  The 
apparent  activation  energy  for  diffusion  was  found  to  be  approximately  29-41  ±  1.5  kJ/mol.  The  trap 
binding  energy  for  y’  precipitates  (10.2  ±  4.6  kJ/mol)  was  found  to  agree  with  previous  studies  [44,  46] 
and  lower  than  the  activation  energy  of  perfect  lattice  diffusion  (25.6  ±  0.5  kJ/mol)  [43].  Despite  all  of  the 
work  performed  previously  [43,  44],  pertaining  to  hydrogen  interaction  in  Monel  K-500  alloys,  there  is  a 
heat  to  heat  variation  that  needs  to  be  addressed  since  it  affects  the  HE  susceptibility  of  this  alloy. 

This  manuscript  aims  to  address  some  of  the  possible  effects  of  heat  to  heat  variations  on  the 
susceptibility  of  Monel  K-500  to  HEAC  with  emphasis  on  hydrogen  material  interactions.  A  previous 
study  addressed  possible  origins  of  differences  from  a  fracture  mechanics  perspective  [36].  In  this  study, 
slow  strain  rate  tests  (SSRT)  of  notched  samples  pre-charged  with  hydrogen  during  48  h  in  0.6  M  NaCl 
solution  at  different  hydrogen  overpotentials  elucidated  differences  in  terms  of  fractography  morphology, 
local  breaking  stress  and  time  to  failure  (TTF)  among  the  heats.  In  air  and  at  a  hydrogen  overpotential  of  - 
0.137  V  the  fractography  was  characterized  by  micro-void  ductile  fracture,  while  during  polarization  to 
more  negative  potentials  the  failure  mode  changed  to  a  mixture  of  ductile/intergranular  (IG)  fracture  at  - 
0.237  V  and  finally  to  IG  fracture  at  hydrogen  overpotentials  of  -0.287  and  -0.387  V.  To  explain  this 
behavior,  the  heats  were  characterized  in  terms  of  mechanical  properties,  microstructure,  bulk  and  grain 
boundary  composition,  and  H-interactions.  Parameters  such  as  hardness,  grain  size,  grain  boundary 
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character,  impurity  segregation,  hydrogen  production  rate,  total  and  diffusible  hydrogen  concentration, 
hydrogen  diffusion  and  trap  binding  energies  were  determined  for  the  different  heats  and  will  be 
discussed  in  terms  of  their  possible  effects  on  hydrogen  embrittlement  (HE)  susceptibility.  The  discussion 
presented  herein  will  provide  insightful  information  pertaining  to  the  effect  of  heat  to  heat  variations  on 
the  susceptibility  of  the  Monel  K-500  alloy  to  HE. 


Methods 

Materials  and  Hardness 

Four  different  Monel  K-500  alloy  (UNS  N05500)  heats  with  compositions  shown  on  Table  1 
were  tested.  One  of  them  was  supplied  by  Allegheny  Technologies  Incorporated.  This  Monel  K-500 
material  was  hot-finished  and  aged,  specifically,  it  was  aged  at  593  °C  (1 100  F)  for  16  hours,  then  furnace 
cooled  at  25  F/hour  to  482  °C  (900  F)  followed  by  air  cooling.  This  condition  is  labeled  as  Allvac. 
Another  heat  was  harvested  from  the  field  after  exposure  to  a  marine  environment.  There  is  no 
information  about  exact  heat  treatment  nor  exposure  time  but  the  material  was  aged.  This  heat  is  labeled 
as  TR  2.  Two  other  heats  were  supplied  by  the  U.S.  National  Research  Laboratory.  These  heats  were  heat 
treated  [49]  to  obtain  upper  and  lower  bounds  of  QQ-N-286G  federal  specification  [3].  These  heats  were 
labeled  as  low  strength,  NRL  LS  (Special  Metals  heat  M5437KG12),  and  high  strength,  NRL  HS  (Special 
Metals  heat  M60N4KG14),  specimens.  The  heat  treated  specimens  were  of  4  inch  and  4.5  inch  diameter 
bars  for  NRL  LS  and  NRL  HS,  respectively.  NRL  LS  was  hot  rolled,  continuously  annealed  at  982  °C 
(1800  F)  followed  by  water  quench,  rotatory  straightened,  rough  turned  (0.250  inch  from  diameter)  then 
3-point  straightened  and  aged.  The  aging  procedure  corresponds  to  593  °C  (1 100  F)  during  2  h,  furnace 
cooled  to  482  °C  (900  F)  at  100  F/h  and  then  air  cooled.  The  NRL  HS  was  direct-aged.  The  aging 
procedure  w 'as  carried  at  593  °C  (1 100  F)  during  16  h,  furnace  cooled  at  25  F/h  up  to  538  °C  (1000F) 
during  6  hours  then  furnace  cooled  at  25  F/h  to  900F  for  9  h  and  air  cooled.  All  NRL  specimens  were 
initially  exposed  as  0.750  inch  (1.905  cm)  thick  compact  tension  discs  to  seawater  under  cathodic 
polarization  from  -0.750  to  -1.1  V  vs  Ag/AgCl  before  retrieval.  Specimens  were  baked  for  48  hours  at 
450  °C  prior  to  laboratory  tests  presented  herein.  The  hardness  of  each  heat  was  measured  at  least  5  times 
on  Rockwell  C  or  B  Scale,  the  average  value  is  reported  on  Table  2. 

Microstructure 
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The  microstructure,  grain  boundary  nature  and  texture  were  evaluated  using  a  FEI  Quanta  650 
scanning  electron  microscope,  SEM,  based  electron  backscattered  diffraction,  EBSD  with  a  beam 
accelerating  voltage  of  30  kV  and  a  probe  current  of  5  nA.  These  settings  allow  an  angular  resolution  less 
than  1°  and  a  spatial  resolution  less  than  30  nm  in  any  direction  and  less  than  10  nm  axially.  The 
maximum  theoretical  error  in  any  measured  feature  size  is  less  than  10%.  This  error  was  calculated  using 
the  SEM  parameters.  Representative  areas  of  each  heat  were  chosen  to  acquire  the  Orientation  imaging 
maps  (OIM).  For  this  purpose,  specimens  were  cut  and  grinded  with  SiC  paper  up  to  1200  grit,  then 
polished  using  3  and  1  pm  glycerol  polycrystalline  diamond  suspensions,  followed  by  vibratory  polishing 
in  propylene  glycol  alumina  suspension  of  0.05  pm  to  reduce  surface  deformation.  The  surfaces  were 
cleaned  using  cotton  balls  soaked  with  ivory  soap  to  remove  alumina  residue  and  then  rinsed  with  d.i. 
water  followed  by  ethanol.  The  samples  were  dried  with  compressed  air.  EBSD  results  are  reported 
elsewhere  [36].  The  average  grain  size  is  reported  in  Table  2. 

Slow  Strain  Rate  Tensile ,  SSRT,  Test 

SSRT  tests  were  performed  on  all  the  Monel  K-500  heats  using  a  Cortest,  Inc  frame  at  a  cross¬ 
head  speed  on  the  order  of  10"6  mm/s  until  failure.  The  SSRT  specimens  were  machined  into  a  notched 
cylinder  geometry  with  dimensions  shown  on  Figure  17.  The  tests  were  performed  in  deaerated  0.6  M 
NaCl  at  25  °C.  Prior  SSRT  testing  the  specimen  was  baked  at  450  °C  during  48  h  to  desorb  the  hydrogen 
present  inside  the  material.  Then,  it  was  charged  using  a  three  electrode  electrochemical  cell  setting  with  a 
Pt  mesh  as  counter  electrode  and  a  Saturated  Calomel  Electrode,  SCE,  as  reference  electrode.  This  cell 
was  designed  to  allow  charging  in  situ  prior  and  during  SSRT  testing.  The  grips  were  masked  utilizing 
Micro  Super  XP-2000  Stop-off  Lacquer  to  avoid  galvanic  coupling  with  the  specimen.  The  potential  was 
controlled  using  a  research  grade  Gamry  Reference  600  potentiostat  during  48  hours  at  -1.1  VSce  prior 
and  during  SSRT  testing.  After  specimen  fracture,  both  fracture  surfaces  were  thoroughly  rinsed  with  tap 
water,  ultrasonically  cleaned  in  D.I.  water  for  20  minutes,  followed  by  ultrasonic  cleaning  in  ethanol  for 
20  mins  and  dried  with  compressed  air.  One  fracture  surface  was  analyzed  using  the  SEM,  while  the  pair 
was  send  for  surface  analysis  of  the  grain  boundary  chemistry  using  Auger  Electron  Spectroscopy,  AES. 

AES 


AES  depth  profile  analyses  were  carried  on  using  a  Physical  Electronics  USA,  Inc.  (PHI)  model 
680  scanning  auger  microprobe,  with  a  Schottky  field  emission  electron  source  and  a  Cylindrical  Mirror 
Analyzer  (CMA).  Auger  spectra  were  collected  using  a  10  kV  and  20  nA  primary  electron  beam.  The 
depth  profile  on  the  grain  boundary  was  achieved  thru  calibrated  sputtering  the  specimen  surface  in  a 
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localized  grain  boundary  with  a  rastered  (3mm  x  3mm)  argon  beam  at  3kV.  The  atomic  fraction  of  each 
element  x,  AFX,  was  calculated  according  to, 


AF  —  lx/Sx 
x  Z  On/Sn) 


Eq.  5 


where  Ix  is  the  intensity  measured  for  element  x  in  the  counts  per  second  versus  kinetic  energy  spectrum 
and  Sx  is  the  sensitivity  factor  of  element  x. 


Electrochemical  Techniques 

A  three  electrode  cell  setting  with  a  Pt  mesh  as  counter  electrode  and  a  SCE  as  reference 
electrode  was  used  for  all  electrochemical  tests.  The  potential  was  controlled  using  a  potentiostat.  The 
electrolyte  used  for  all  the  experiments  was  0.6  M  NaCl  pH  8.0  by  adjustment  with  NaOH.  Cathodic 
potentiodynamic  polarizations  were  performed  on  all  the  Monel  K-500  heats  at  a  scan  rate  of  0.167 
mV/s. 


Electrochemical  Hydrogen  Charging 

To  evaluate  the  hydrogen  uptake  of  the  different  Monel  K-500  heats,  120  ±  10  pm  planar  thick 
specimens  were  ground  to  600  grit  surface  finish,  rinsed  with  D.I.  water,  degreased  with  ethanol 
ultrasonically,  dried  with  compressed  air  and  charged  with  hydrogen  via  cathodic  potentiostatic 
polarization  at  different  potentials  (-0.850,  -0.950,  -1.0  and  -1.1  VSCE)  in  0.6  M  NaCl  (pH  8.0  by 
adjustment  with  NaOH)  during  7  days.  In  the  case  of  determination  of  binding  energies  a  potential  of -1.0 
Vsce  was  used  during  charging.  The  charging  time  and  specimen  thickness  to  achieve  approximately  90% 
hydrogen  saturation  of  the  specimen  were  chosen  based  on  diffusion  calculations  assuming  that  surface 
concentration  and  hydrogen  diffusion  coefficient  (0.9  -  3.9  x  10'14  m2/s  [43])  were  constants,  semi-infinite 
solid  and  slab  geometry  [50].  Prior  hydrogen  charging,  specimens  were  grinded  to  600  grit,  followed  by 
ultrasonic  cleaning  in  ethanol  for  one  minute  and  finally  dried  with  compressed  air. 

Determination  of  Diffusible  and  Total  Hydrogen  Concentration  through  Electrochemical  Extraction  and 
LECO ,  respectively 

The  diffusive  hydrogen  concentration,  CH  Dim  was  determined  using  an  electrochemical  extraction 
method  [43,  51].  This  technique  is  based  on  the  electrochemical  method  of  hydrogen  permeation 
developed  by  Devanathan  et  al.  [52]  After  charging,  the  specimen  was  slightly  abraded  with  600  grit  SiC 
paper  and  then  rinsed  with  D.I.  Water  and  ethanol  followed  by  drying  with  compressed  air.  The  hydrogen 
diffusing  out  from  the  specimen  was  oxidized  at  -0.85  VMmse  in  de-aerated  Borate  Buffer  (pH:  10)  during 
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one  hour,  even  though  only  the  first  half  hour  of  the  data  was  used  for  calculations  (-0.85  VMmse  is  more 
positive  than  -1.21  VMmse,  which  is  the  Nernst  potential  for  hydrogen  at  pH:  10).  The  Barnacle  Cell  test 
was  performed  at  least  three  times  on  each  Monel  K-500  heat.  CHlDifr  was  calculated  following  [51], 


-'H.Diff 


_  Jt 

zF\nt) 


-1/2 


Eq.6 


where  Jt  is  the  hydrogen  oxidation  current  density  at  time  t  (a/cm2);  z  is  the  number  of  electrons  involved 
in  the  reaction  (z:  1  equivalent  per  mole  for  hydrogen  evolution  reaction);  F  is  the  Faraday  constant 
(96,498  C/  equivalent)  and  Deff  is  the  effective  hydrogen  diffusivity  (a  value  of  1.3  x  10'14  m2/s  was 
utilized  according  to  those  found  in  the  literature  for  Monel  K-500  alloys  [43]).  Equation  6  is  only  valid 
when  L2/Dt  >4  [51],  where  L  is  the  depth  from  which  the  hydrogen  is  extracted.  In  addition,  was 
assumed  to  be  independent  of  CH.Dtfr,  fact  that  was  verified  previously  [43].  The  CH/r0tai  was  determined 
from  LECO  hydrogen  measurements  heating  to  a  maximum  temperature  of  650  °C  according  to  ASTM 
Standard  E1447-09/CTP  3008-1/IG  [53]. 


Determination  of  Total  Hydrogen  Concentration  and  Binding  Energy  through  Thermal  Desorption 
Spectroscopy 


A  Thermal  Desorption  Spectroscopy  system  described  previously  [54,  55]  was  utilized  to 
determine  the  trap  binging  energy,  EB,  and  the  effective  hydrogen  diffusivity,  Dcff.  Desorbed  hydrogen 
was  measured,  under  high  vacuum,  as  hydrogen  partial  pressure,  PH2,  during  programmed  ramped  heat 
experiments  to  calculate  the  parameters  previously  mentioned.  Backgrounds  were  measured  using  the 
residual  gas  in  the  system  without  any  specimen  at  different  heating  rates. 

The  Deff  was  calculated  taking  into  consideration  the  dependence  of  the  increase  in  integrated- 
absorbed  hydrogen  level  in  a  plate  of  known  thickness,  21,  on  the  Deff.  The  solution  of  the  diffusion 
equation  by  a  membrane  under  non-steady  state  conditions  with  equal  surface  concentrations  and  uniform 
initial  concentration  gives  [56]: 

£=  ^-Zn^^^exp{-D(.2n  +  l)2n2tm2)}  Eq.7 

where,  Mt  denotes  the  total  mass  of  H  absorbed  by  the  membrane  at  time,  t,  and  M*  indicates  the  total 
mass  of  H  absorbed  after  infinite  time.  Experimental  data  was  gathered  for  specimens  charged  at  different 
times  and  tested  for  absorbed  hydrogen  content  by  thermal  desorption  using  a  heating  rate  of  10  °C/min 
up  to  a  temperature  of  550  °C.  The  data  was  background  corrected,  the  hydrogen  desorbed  was  integrated 
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with  respect  to  heating  time  and  compared  to  theoretical  curves  ofMt/M00  vs  t12  using  different  Deff 
values  to  choose  the  curve  that  best  fitted  the  experimental  data. 

For  the  case  of  a  cylinder  geometry  for  NRL  heats,  the  solution  of  the  diffusion  equation  for  non¬ 
steady  state  is  given  by  [56], 


Eq.  8 


where,  a  is  the  cylinder  radius  and  a  is  the  positive  root  of  J0(aan)  =  0  being  J0  the  Bessel  function  of  the 
first  kind  of  order  zero.  In  this  case,  the  NRL  samples  were  pre-charged  at  different  time  intervals  at  - 
0.800,  -0.850,  -0.890  and  -1.000  VSCe  in  seawter.  Theoretical  calculations  were  performed  for  hydrogen 
diffusivities  of  5  x  10'10  cm2/s,  1  x  10‘10  cm2/s  and  2  x  10'10  cm2/s. 

Experiments  at  different  heating  rates  of  3,  5,  7.5  and  10  °C/min  were  performed  to  determine  the 
Eb.  For  this  purpose,  120  pm  ±  10  pm  thick  flat  plate  specimens  were  electrochemically  charged  with 
hydrogen  in  0.6  M  NaCl  (pH:  8.0  by  adjustment  with  NaOH)  during  7  days  at  -1.0  VSce-  This  parameter 
can  be  estimated  from  the  H  desorption  energies,  Ed.  It  is  known  that  the  Ed  is  equal  to  the  sum  of  the 
activation  energy  for  H  migration  from  an  interstitial  lattice  site  to  a  trap,  Em,  and  the  EB.  Now,  assuming 
that  Em  is  equal  to  the  activation  energy  for  H  diffusion  by  an  interstitial  jump  mechanism  in  the  trap-free 
or  perfect  lattice,  Em,  then  Ed  =  Em  +  Eb.  This  relationship  between  the  energies  is  only  valid  if  detrapping 
and  lattice  diffusion  are  the  dominant  processes  during  H  desorption.  The  Ed  for  a  specimen  with  a 
homogenous  H  concentration  across  its  thickness  is  defined  by  [57]: 


Eq.  9 


where  x  is  the  fraction  of  the  desorbed  H  quantity  integrated  from  a  single  desorption  site  and  therefore  a 


single  deconvoluted  TDS  peak,  A  is  a  constant,  n  is  the  reaction  order,  R  is  the  gas  constant  (8.31  J/mol 
K)  and  T  is  the  temperature  in  K.  Equation  9  provides  a  reasonable  approximation  of  Ed  when  the 
following  conditions  are  met:  fast  diffusion  or  slow  diffusion  taken  into  account,  existence  of  a  single 
dominant  trap,  minimal  retrapping  and  low  trap  occupancy  [58].  When  measuring  the  desorbed  hydrogen 
utilizing  TDS,  different  desorption  peak  temperatures  Tm  are  measured.  These  values  depend  on  the 
heating  rate  [59,  60].  These  Tmand  the  utilized  heating  rate  values  were  used  to  calculate  Ed  using  the 
following  equation  [57],  which  is  valid  for  detrapping  and  diffusion-controlled  egress  under  certain 
conditions  [61]. 
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Ed 


R 


Eq.  10 


The  data  was  plotted  as  In  vs  1  /Tm  and  then  the  Ed  was  calculated  from  the  slope  of  the 

*m 

linear  regression  fitted  to  the  experimental  data  points  with  a  confidence  interval  of  95%.  The 
characteristic  desorption  spectra  was  always  composed  of  two  peaks,  which  were  deconvoluted  based  on 
Gaussian  distribution  model  using  the  OriginPro  7.5  software  to  determine  Tm.  The  lowest  Tm  value  was 
associated  with  the  hydrogen  trapped  in  the  lattice  in  accordance  with  trap  theory  [62]. 


Determination  of  Grain  Boundary  Trap  Binding  Energy > 

The  calculation  procedure,  in  the  case  of  limited  experiments,  for  the  Ed  was  performed  utilizing 
the  following  equation  [63], 

Ed  =  RTmax(BTmax  -  1)  Eq.  1 1 

where,  R  is  the  gas  constant,  Tmax  is  the  maximum  peak  temperature,  and  B  is  a  constant  that  is  material- 
dependent.  B  was  calculated  from  data  presented  in  a  previous  publication  for  SHT  +  Aged  Monel  K-500 
alloy  [43].  Then,  the  subtration  of  the  low  temperature,  Edi,  from  the  high  temperature  desorption  peak 
energy,  E^  represent  the  Eb  for  the  grain  boundary.  However,  to  detect  the  binding  energy  of  the  grain 
boundary  without  interference  from  other  traps,  the  y*  particles  were  solutionized  in  accordance  with  the 
Special  Metals  data  sheet  [1]  at  a  selected  temperature  of  1010  °C  during  forty  minutes  and  WQed  to 
impede  precipitation.  The  specimen  was  extracted  from  the  TR  2  heat  and  prepared  to  a  thickness  of 
300  pm.  Heat  treatments  performed  during  shorter  periods  of  time  were  tested  to  make  sure  that  the 
hardness  of  the  specimen  decreased  due  to  precipitate  dissolution  and  that  the  grain  boundaries  did  not 
grow.  Then,  the  specimens  were  electrochemically  charged  with  hydrogen  in  0.6  M  NaCl  (pH:  8.0  by 
adjustment  with  NaOH)  during  7  days  applying  a  potential  of  -1.0  VSce-  One  of  these  specimens  was 
analyzed  using  the  TDS  right  after  charging  and  another  one  was  tested  after  hydrogen  desorption  at  25 
°C  during  6  days. 
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Results 


Time  to  Failure  and  Maximum  Remote  Stress  after  SSRT 

Figure  18  displays  the  time  to  failure  (TTF)  and  the  maximum  remote  stress  of  the  different  heats 
during  concurrent  charging/SSRT  test  after  electrochemical  FI  charging  in  deaerated  0.6  M  NaCl  (pH:  8.0 
by  adjustment  with  NaOH)  during  48  hours  at  -0.850,  -0.950,  -1.0  and  -1.1  VSce-  The  average  TTF  ranged 
between  14.30  and  34.14  h.  The  effect  of  cathodic  polarization  on  TTF  was  evaluated  testing  TR  2  and 
Allvac  heats.  TTF  was  found  to  decrease  as  the  cathodic  polarization  potential  decreases  for  both  heats, 
which  is  consistent  with  the  expectation  of  intergranular  HEAC  at  high  hydrogen  overpotentials.  In 
addition,  analysis  of  this  parameter  shows  a  higher  susceptibility  of  TR  2  to  HEAC  than  Allvac  at  every 
potential  tested,  except  at  -0.850  Vsce-  At  this  potential  the  TTF  is  the  same  for  both  heats.  The  TTF 
values  for  the  four  heats  are  displayed  in  Figure  2  after  the  test  was  performed  at  -1.1  VSce-  The  results 
elucidate  the  following  order  in  terms  of  TTF,  TR2  <  Allvac  ~  NRL  LS  <  NRL  HS.  This  order  represents 
the  hydrogen  embrittlement  assisted  cracking  (HEAC)  resistance. 

In  terms  of  the  maximum  stress,  this  parameter  was  calculated  taking  into  consideration  the  cross 
sectional  area  at  the  ends  of  the  specimen  next  to  the  thread,  to  yield  the  maximum  remote  stress.  The 
calculated  values  were  in  the  order  of  129.50  to  387.07  MPa.  The  Allvac  lot  was  always  found  to  fail  at  a 
higher  maximum  remote  stress  than  TR  2  under  cathodic  polarization,  while  the  opposite  was  found  under 
benign  conditions  (in  air).  All  lots  were  tested  at  an  applied  potential  of -1.0  VSce  and  it  was  found  that 
the  HEAC  most  resistant  alloy  was  NRL  HS  followed  by  Allvac,  NRL  LS  and  TR  2,  which  correlates 
well  with  the  TTF  results  at  the  same  applied  potential.  It  is  interesting  to  note  that  the  TTF  did  not  trend 
with  hardness  as  discussed  below. 

SSRT  Fractography 

The  fractography  of  the  edge  of  the  notched  cylindrical  specimens  subjected  to  SSRT  for  the 
Allvac  and  TR  2  lots  is  presented  in  Figure  19.  The  failure  modes  for  the  benign  (in  air)  and  less 
aggressive  wet  conditions  (-0.850  VSce)  are  characterized  as  micro-void  dimple  fracture  for  both  lots, 
while  under  the  most  aggressive  condition  (-1.1  Vsce)  the  failure  mode  is  characterized  as  inter-granular 
(IG)  fracture  at  the  edge  of  the  sample  with  a  ductile  stress  overload  in  the  center  of  the  specimen  for  both 
lots.  However,  the  failure  mode  becomes  distinctive  at  an  applied  potential  of  -0.950  VSce-  Under  this 
condition  the  fractography  is  represented  as  micro-void  dimple  fracture  for  Allvac  lot,  while  it  is 
characterized  by  a  mixed  micro-void  dimple/IG  fracture  for  the  TR  2  lot.  This  observation  confirms  the 
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TTF  and  maximum  remote  stress  results  that  indicate  a  higher  HEAC  susceptibility  for  TR  2  when 
compared  to  Allvac. 

The  fractography  for  NRL  LS  and  HS  is  displayed  in  Figure  20.  The  fractography  revealed 
differences  in  the  fracture  behavior  of  these  lots.  The  failure  mode  was  characterized  by  IG  cracking  for 
NRL  LS  lot,  while  it  consisted  of  a  mixed  IG/trans-granular  (TG)  failure  mode  for  NRL  HS  lot. 

Local  Chemistry  Characterization  of  Grain  Boundaries  via  AES  after  SSRT  at  - 1.1  VSce  and  Trace 
Elements  present  in  the  Bulk 

After  SSRT  test  at  an  applied  potential  of  -1.1  VSce  at  least  two  grain  boundaries  from  each  heat 
were  chemically  characterized  using  AES.  Figure  21a  presents  an  example  of  the  spectra  taken  for  at  least 
two  grain  boundaries  on  the  fracture  surface  of  the  different  heats  after  the  test,  the  concentration  profile 
characterization  as  atomic  percent  and  the  normalized  concentration  profile  characterization  with  respect 
to  a  depth  of  21  nm  (near  the  bulk  value)  to  address  the  fluctuation  of  the  elemental  concentration  at  the 
grain  boundary.  Figure  21a-d  display  the  normalized  concentration  profiles  for  the  different  heats.  All  the 
grain  boundaries’  surfaces  presented  a  superficial  layer  composed  of  O,  C  and  S.  The  oxygen  and  carbon 
peaks  are  the  result  of  the  C  deposition  from  the  atmosphere  that  the  fracture  surfaces  experience  as  soon 
as  the  SSRT  setup  is  de-assembled  and  handled.  In  terms  of  the  oxygen,  this  element  is  present  as  result 
of  the  formation  of  an  air  oxide  as  soon  as  the  new  surfaces  are  depolarized  and  are  in  contact  with  the  air 
or  solution.  All  heats  indicated  sulfur  enrichment  represented  by  a  S  peak  at  depths  between  0.93  and  2 
nm  as  shown  in  Figure  21a-d  and  Table  3.  The  S  concentration  was  also  evaluated  as  a  raw  and 
normalized  S/Ni  ratio  to  obtain  a  more  representative  level  when  analyzing  the  sulfur  level  on  grain 
boundaries  where  different  elements  were  measured.  The  sulfur  composition  at  the  grain  boundary  for  the 
different  heats  was  found  to  be  in  between  1.71  and  4.12  at%.  The  maximum  value  was  found  on  the 
Allvac  heat,  which  does  not  correlate  with  HEAC  resistance  results.  When  the  S  composition  was 
normalized  to  a  depth  of  21  nm,  the  maximum  level  of  sulfur  was  found  on  NRL  LS,  which  it  is  not  the 
most  HEAC  susceptible  heat.  Another  approach  to  interpret  the  data  is  to  use  the  ratio  of  the  element  with 
respect  to  the  base  element  to  avoid  any  variations  attributed  to  the  presence  of  small  amounts  of  different 
elements  among  the  heats.  In  this  case,  the  maximum  S/Ni  values  did  not  correspond,  in  some  instances, 
to  the  maximum  S  concentration  depth.  However,  the  maximum  S/Ni  level  agreed  with  the  maximum  S 
concentration  and  was  found  on  the  Allvac  grain  boundary,  but  did  not  correlate  with  the  HEAC 
susceptibility.  This  lack  of  correlation  of  S  segregation  with  FtEAC  susceptibility  was  also  confirmed 
with  observation  of  the  S/Ni  ratio  normalized  to  a  depth  of  21  nm. 
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Error!  Reference  source  not  found,  presents  the  values  of  surface  area  of  grain  boundary  to 
volume  ratio  [36];  trace  elements  content  as  expressed  above  [36]  pertaining  to  the  deleterious  or 
enhanced  effect  of  these  elements  with  respect  to  HE  susceptibility.  The  deleterious  and  gettering  effects 
of  the  trace  elements  were  evaluated  utilizing  Equations  1  and  3.  The  equilibrium  sulfur  concentration  as 
described  elsewhere  was  found  to  be  in  the  range  5.2  to  20  ppm,  being  the  most  HEAC  resistant  heat, 
NRL  HS,  the  material  w  ith  the  highest  equilibrium  sulfur  concentration.  In  terms  of  the  gettering  factor,  a 
higher  value  represents  an  enhanced  circumstance  where  the  available  sulfur  in  the  microstructure  would 
have  a  higher  probability  of  being  gettered  and  therefore  can  noy  segregate  to  grain  boundaries.  The 
gettering  factor  was  found  within  49  and  86.3  ppm,  which  represent  the  values  for  Allvac  and  NRL  LS 
heats,  respectively.  These  results  alone  do  not  correlate  well  with  HEAC  susceptibility  evaluated  through 
SSRT  metrics.  Therefore,  taking  into  consideration  the  deleterious  or  beneficial  effects  of  certain  trace 
elements,  an  equation  was  formulated  to  try  to  explain  the  effect  of  trace  elements  on  the  segregation  and 
gettering  of  elements  and  its  interaction  with  HEAC  susceptibility.  The  equation  used  gives  the  same 
weight  to  S,  Sb,  Sn,  as  detrimental  elements,  and  to  Hf  as  beneficial  element  in  combination  with  a  third 
of  the  weight  of  Zr  (also  a  beneficial  element).  The  sum  of  the  deleterious  elements  was  given  by  (+) 
while  beneficial  elements  were  subtracted  (-).  This  yielded  a  single  bulk  compositional  parameter  in  an 
attempt  to  shed  light  on  the  combined  effect  of  trace  elements  on  HE  susceptibility.  The  S+Sn+Sb-Hf- 
Zr/3  values  were  TR-2  -97.0;  NRL  LS  -  76.7;  Allvac  -121.3;  and  NRL  HS  -197.9  when  expressed  in  ppm 
concentrations  given  as  a  sum.  Here  the  maximum  value  of  this  parameter  was  obtained  for  the  heat  more 
susceptible  to  HE.  The  use  of  this  factor  for  the  different  heats  evaluated  provided  the  following  order, 
NRL  HS  <  Allvac  <  NRL  LS  <  TR  2,  which  correlates  exactly  with  the  HE  susceptibility  trend  found 
previously  by  SSRT  (Figures  2  and  3).  In  addition,  the  surface  area  of  grain  boundary'  and  hence  trap  and 
segregation  site  density  on  boundaries  compared  to  the  volume  ratio  followed  the  trend  TR  2  <  NRL  LS  < 
Allvac  <  NRL  HS,  which  also  agreed  with  the  HE  susceptibility  order  assessed  by  TTF  metric  from 
SSRT  experiments.  Hence,  the  segregation  parameter  in  combination  with  the  density  of  grain  boundary 
trap  or  segregation  sites  per  unit  volume  appears  to  indicate  the  susceptibility  . 

Electrochemically  driven  Hydrogen  Production  on  Heats  Surfaces 

The  hydrogen  evolution  reaction  kinetic  was  addressed  by  calculation  of  the  Tafel  slope  from 
potentiodynamic  polarizations  performed  in  the  cathodic  regime  in  deaerated  0.6  M  NaCl  at  pH  8.0  by 
adjustment  with  NaOH  (Figure  23a).  The  cathodic  polarization  presented  two  distinct  kinetic  behaviors. 
At  the  higher  potentials  mass  transfer  control  prevails  until  a  potential  of  -0.8  Vsce>  approximately,  where 
there  is  a  transition  to  charge  transfer  control.  This  critical  potential  is  below  the  Nernst  potential  for 
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hydrogen  evolution  at  the  electrolyte  pH  (-0.713  VSce).  The  Tafel  slope  was  calculated  from  the  charge 
transfer  control  portion  of  the  curve  as  shown  in  Figure  23b.  The  experimental  hydrogen  production  rates 
were  found  to  be  similar  among  the  different  heats  and  between  -0.233  and  -0.222  V/decade.  The  heat  to 
heat  variations  did  not  affect  this  parameter. 

Diffusible  and  Total  Hydrogen  Concentration  via  Barnacle  Cell  and  LECO,  respectively 

Regarding  the  hydrogen  concentrations  and  trapping,  Figure  24  displays  the  effect  of  hydrogen 
overpotential  on  diffusible  and  total  hydrogen  concentration  for  the  different  heats  after  electrochemically 
charging  with  hydrogen  during  7  days  in  0.6  M  NaCl  (pH:8.0  by  adjustment  with  NaOH).  Both 
concentration  parameters  were  found  to  increase  with  hydrogen  overpotential.  The  Allvac  heat  was 
observed  to  exhibit  lower  total  and  diffusible  hydrogen  concentrations  over  a  range  of  cathodic  potentials. 
As  expected,  the  total  hydrogen  concentration  CHj0tai  was  larger  than  CH,oiff  as  typically  observed  [43]. 
This  behavior  is  related  to  the  fact  that  hydrogen  not  only  interacts  with  the  crystal  lattice,  but  also  with 
low  and  high  binding  traps,  such  as:  dislocations  [44,  45],  grain  boundaries,  carbides  [42],  microvoids, 
age  hardening  precipitates  like  f  [43,  44,  64],  among  others.  In  terms  of  the  heat  to  heat  variation,  for 
both,  CH,Totai  and  CH)Diff  there  is  limited  difference,  except  for  the  Allvac  heat.  This  heat  consistently 
exhibited  hydrogen  concentration  values  near  the  lower  bound  of  all  other  heats  or  statistically  lower 
hydrogen  levels  based  on  repeated  measurements  at  some  overpotentials. 

Trap  Binding  States  and  Desorption  and  Binding  Energies 

The  total,  lattice,  as  well  as  the  f  trap  binding  energy  of  the  different  heats  were  determined  from 
ramped  TDS  experiments  executed  at  different  heating  rates  (3,  5,  7.5  and  10  °C/min).  Figure  25  displays 
the  normalized  hydrogen  pressure  versus  time  for  as-received  and  baked  TR2  specimens.  Hydrogen  is 
indicated  by  peaks  with  one  or  more  maxima  observed.  In  this  case,  the  experiment  elucidates  the  intake 
of  hydrogen  by  the  TR  2  heat  while  it  was  exposed  in  seawater  under  cathodic  polarization,  while  the  TR 
2  specimen  baked  at  450  °C  indicates  no  peak.  This  baking  procedure  was  performed  in  all  the  samples 
during  48  hours  prior  electrochemical  charging  for  all  the  H-interaction  tests  mentioned  in  this 
manuscript,  except  for  the  experiments  pertaining  to  the  calculation  of  the  hydrogen  production  rate. 

The  hydrogen  desorption  rate  data  for  the  Allvac  (Figure  26a),  TR  2  (Figure  27a),  NRL  LS 
(Figure  28a)  and  NRL  HS  heats  (Figure  29a)  presented  the  same  characteristic  behavior  for  all  the  heats, 
that  is,  a  curve  with  two  maxima,  which  represent  two  hydrogen  states  or  trap  states.  The  peaks  identified 
on  each  heat  for  the  hydrogen  desorption  rate  as  function  of  temperature  experimental  data  were 
deconvoluted  using  OriginPro  7.5  software.  The  presence  of  these  two  states  in  the  Monel  K-500  heats 
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tested  is  in  agreement  with  previous  information  reported  on  solution  heat  treated  (SHT),  aged  and  SHT, 
cold  worked,  aged  Monel  K-500  materials  [43],  These  peaks  are  also  present  for  samples  pre-charged  at 
all  the  different  potentials,  mentioned  previously,  as  shown  in  Figure  30  for  example  when  all  the  heats 
were  pre-charged  at  a  potential  of  -1.1  VSce-  The  peak  at  higher  temperature  was  previously  identified  as 
the  y’  interface  trap  and  it  was  found  to  be  smaller  than  the  area  under  the  left  peak  at  lower  temperatures 
attributed  to  lattice  hydrogen  and  low  energy  reversible  traps  [43].  The  area  under  both  peaks  for  this  plot 
of  hydrogen  desorption  rate  vs  temperature  represents  the  total  hydrogen  concentration  for  each  heat, 
which  follows  this  order,  Allvac  <  NRL  HS  <  NRL  LS  <  TR  2.  This  technique  corroborates  the  finding 
that  Allvac  was  the  heat  with  less  diffusible  hydrogen  concentration,  with  respect  to  the  other  heats,  after 
pre-charging  at  the  same  condition  in  electrochemical  extraction  (Figure  24).  The  desorption  peak 
temperatures  measured  with  TDS  at  the  different  heating  rates  after  pre-charging  at  -1.0  VSce  were  found 
to  increase  with  the  heating  rate,  which  is  the  expected  behavior  for  a  thermally  controlled  desorption 
process  [65,  66].  This  behavior  has  been  reported  previously  [43].  The  different  Ed  values  were  calculated 
from  the  slope  of  the  linear  regression  of  the  heating  rate  -  desorption  peak  temperatures  pairs  for  Allvac 
(Figure  26b),  TR  2  (Figure  27b),  NRL  LS  (Figure  28b)  and  NRL  HS  heats  (Figure  29b)  and  were 
reported  at  the  95  %  confidence  interval  for  each  heat.  The  left  and  right  Ed  values  are  presented  in  Table 
4  and  show  that  Edj  average  values  range  in  between  23  and  34  kJ/mol  and  follow  Allvac  <  TR  2  ~  NRL- 
LS  <  NRL-HS.  The  value  range  presented  in  this  investigation  approximately  overlaps  with  that  reported 
previously  for  Monel  K-500  alloys  [43].  The  y‘  binding  energy  was  calculated  from  the  subtraction  of  Edj 
from  Ed  r.  The  average  values  for  the  binding  energies  for  y4  were  found  in  the  range  7.75  to  19.1  kJ/mol 
in  the  following  ascending  order:  NRL  HS  <  TR  2  <  Allvac  <  NRL  LS  (See  Table  4). 

Additional  experiments  focused  on  elucidation  of  grain  boundary  trap  binding  energy.  The  results 
and  a  comparison  with  the  original  TR  2  heat  are  presented  in  Figure  31.  The  resolutionized  materials 
present  two  trap  binding  states,  one  at  approximately  200  °C,  which  agrees  with  the  left  peak  assigned  to 
lattice  sites  for  the  non-resolutionized  material  and  another  one  at  350  °C,  which  should  correspond  to  the 
grain  boundary  trap  sites  since  the  peak  representing  the  y‘  trap  sites  is  no  longer  present  as  result  of  the 
SHT  +  WQ  heat  treatment  at  1010  °C  during  40  mins.  The  Ed  for  the  350  °C  peak  was  found  to  be  55.0 
kJ/mol  while  the  Ed  for  the  200  °C  peak  was  30.7  kJ/mol,  both  were  calculated  using  Eq.  11  and  the  B 
constant  (0.01863)  calculated  from  data  presented  in  a  previous  publication  for  SHT  +  Aged  Monel  K- 
500  alloy  [43].  Therefore,  subtraction  of  the  high  temperature  Ed  from  the  low  temperature  Ed  results  in  a 
Eb  of  approximately  24.3  kJ/mol  for  trap  sites  assigned  to  grain  boundaries,  which  is  above  the  Eb  for  the 
y‘  trap  for  the  TR  2  heat  (Table  4). 
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Effective  Hydrogen  Diffusivity  Calculation  via  Saturation  Curve  Analysis  utilizing  TDS 

Saturation  curves  were  constructed  measuring  the  total  mass  of  hydrogen  in  pre-charged  NRL 
heats  in  seawater  at  an  applied  potential  of  -0.800,  -0.850,  -0.890  and  -1.000  VSce  during  different  time 
intervals  (Mt)  until  saturation  (M*).  The  M/M*  ratios  are  plotted  against  time  in  Figure  32.  The  time 
needed  for  saturation  was  approximated  utilizing  diffusion  models  for  semi-infinite  cylindrical  geometry 
for  a  saturation  of  90%  and  a  D^efr  of  1.3  x  lO'10  cm2/s  [56].  In  addition,  theoretical  curves  were 
calculated  using  Eq.  8  and  are  presented  in  Figure  32.  The  experimental  data  was  found  to  oscillate 
between  the  theoretical  saturation  curves  for  a  hydrogen  diffusivity  of  5  x  10"n  and  2  x  10'10cm2/s,  except 
for  one  data  point,  and  presents  the  best  fit  for  a  hydrogen  diffusivity  of  1  x  10‘Iucm2/s.  A  few  additional 
experiments  were  performed  for  the  TR  2  and  Allvac  heats  on  a  plane  sheet  geometry  (not  shown  in  the 
plot).  The  data  points  suggest  that  the  hydrogen  diffusivity'  of  these  heats  does  not  exceed  the  one  for 
NRL  materials.  Hence,  differences  in  hydrogen  diffusion  rate  are  not  likely  the  factors  accounting  for 
differences  in  HE  susceptibility 


Discussion 

Factors  effecting  hydrogen  embrittlement  susceptibility 

The  susceptibility  of  Monel  K-500  to  HEAC  is  controlled  by  hydrogen  [15].  Hydrogen 
embrittlement  is  clearly  exacerbated  by  cathodic  potentials  which  increase  the  hydrogen  overpotential  and 
in  turn  hydrogen  uptake  as  discussed  below  [43].  Large  increases  in  trapped  and  lattice  H  as  a  function  of 
hydrogen  overpotential  were  indicated  in  0.6  M  NaCl  in  previous  TDS  investigations  [43].  Moreover, 
prior  publications  predicted  significant  effects  of  hydrogen  on  Kth  and  stage  11  crack  velocity,  [da/dt]n 
[15,  36].  There  is  a  threshold  based  on  CH.Cnt above  which  stage  II  crack  velocity  increases  markedly.  KTh 
decreases  monotonically  at  potentials  negative  to  about  -0.800  to  -0.900  VSce  depending  on  hydrogen 
concentration,  Griffith  fracture  toughness  and  yield  strength  [15,  36]  . 

Recall  that  TR  2  material  was  harvested  from  seawater  deployment  where  HEAC  was  observed 
and  hydrogen  was  detected  upon  retrieval  from  service  (Figure  25).  It  is  clear  from  the  data  indicated 
herein  that  TR  2  and  NRL  LS  exhibit  greater  susceptibility  to  hydrogen  embrittlement  as  evidenced  by 
lower  breaking  stresses,  times  to  failure  in  SSRT  (Figure  18)  as  well  as  by  fractography  (Figure  19c  and  d 
and  Figure  20)  which  indicates  a  greater  extent  of  intergranular  fracture  even  at  less  severe  cathodic 
potentials.  Allvac  exhibits  slightly  less  embrittlement  susceptibility  but  is  clearly  embrittled  for  instance 
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at  -1.1  VSce  while  remaining  somewhat  more  resistant  at  -0.95  and  -0.85  VSce  (Figure  18  and  Figure  19b 
and  c).  This  susceptibility  trend  and  lot-based  difference  was  also  seen  in  fracture  mechanics  based 
environment  assisted  cracking  under  identical  conditions  [36].  The  purpose  of  this  paper  is  to  attempt  to 
identify  factors  controlling  these  lot-to-lot  differences  based  on  metallurgical  factors.  These  issues  are 
more  relevant  at  intermediate  potentials14  as  clearly  all  Monel  K-500  lots  tested  are  embrittled  at  -1.1  VSce 
and  this  is  consistently  observed  in  past  studies  [4,  21,  67].  More  rapid  time  to  failure  and  lower  breaking 
stress  for  TR  2  and  NRL  LS  does  not  appear  to  correlate  with  hardness  as  indicated  in  Figure  33. 

Concerning  hydrogen  parameters,  all  materials  exhibited  similar  hydrogen  production  rates  as 
indicated  by  Figure  23.  However,  TR  2  and  NRL  LS  exhibited  enhanced  CH,diffand  CH, total  relative  to  the 
Allvac  material  but  not  NRL  HS  (Figure  24).  The  trend  in  CH,diff  follows  the  lot-to-lot  hardness  data 
shown  in  Figure  33  with  the  exception  of  the  Allvac  material  which  is  harder  than  NRL  LS  but  absorbed 
markedly  less  hydrogen.  The  trend  towards  an  increase  in  CH,diffwith  hardness  could  correspond  with  a 
higher  density  of  reversible  trap  sites  in  the  harder  material.  The  y’  trap  binding  energy  decreases  with 
increasing  hardness  (Figure  35)  suggesting  some  H  interchange  with  lattice  sites  and  participation  in 
CH.difffcr  these  weaker  sites  [43].  The  decrease  in  y9  binding  energy  with  hardness  could  be  explained  by 
slight  change  in  precipitate-lattice  coherency  or  in  the  density  of  strongly  bound  sites  on  these  particles 
[45].  Other  parameters  in  hydrogen  analysis  are  similar  such  as  hydrogen  diffusion  coefficient.  A 
difference  between  hydrogen  diffusion  in  NRL  heats,  Allvac  and  TR  2  materials  was  not  detected. 

This  trend  in  susceptibility  is  also  supported  by  near-grain  boundary  sulfur  analysis  from  Auger 
Electron  Spectroscopy  and  bulk  chemistry.  First,  it  should  be  noted  that  all  materials  contain  sulfur  (Table 
1)  at  concentrations  above  the  reported  solubility  limits  for  S  in  Ni  at  aging  temperatures  [67,  68].  Also, 
all  alloys  exhibit  near  surface  enhancements  in  sulfur  content  to  various  extents  which  falls  rapidly  with 
depth  as  seen  previously  for  Ni  (Figure  21a-d)  [68,  69].  What  is  interesting  to  note,  though  is  that  the 
alloys  tested  have  drastically  different  ratios  of  grain  boundary  embrittling  species  like  sulfur  and 
beneficial  gettering  species  such  as  Hf,  Zr,  Mg  and  Ca  or  Boron  [12,  23,  34].  These  elements  may  getter 
sulfur  by  forming  sulfo-carbides  [12,  34]  and  sulfides[23].  It  is  interesting  to  note  that  a  combination  of  a 
lower  S+Sb+Sn-Hf-Zr/3  parameter  with  a  low  ratio  of  surface  area  of  grain  boundary  to  the  volume  is 
found  for  the  less  susceptible  materials  such  as  Allvac  (Error!  Reference  source  not  found.).  However, 
it  should  be  noted  that  grain-to-grain  differences  in  near  grain  boundary  surface  peak  sulfur,  S/Ni  ratio 


14  Low  voltage  anodes  are  sought  to  provide  adequate  levels  of  cathodic 
protection  for  steel  structures  yet  not  be  too  severe  so  as  to  embrittle 
MK500  fasteners  in  seawater. 
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and  normalized  quantities  do  not  indicate  a  clear  trend  in  sulfur  segregation  (Table  3).  This  is  likely  due 
to  carbon  contamination  confounding  analysis  and  the  fact  that  sulfur  content  varies  by  grain. 

Each  of  these  three  factors  including  sulfur  segregation,  hydrogen  uptake  and  grain  boundary 
trapping  as  well  as  hardness  are  discussed  in  more  detail  below. 

H  Interactions  and  Uptake 

Given  similar  HER  rates  on  planar  electrodes  as  a  function  of  hydrogen  over  potential,  a 
difference  in  the  hydrogen  production  rate  is  not  likely  the  cause  of  differences  in  hydrogen  susceptibility. 
The  hydrogen  evolution  reaction  on  Ni  in  alkaline  environments  during  cathodic  polarization  has  been 
extensively  investigated  [70-73]  .  The  mechanism  involves  rate  determining  slow  discharge  from  a  water 
molecule  followed  by  rapid  Tafel  (chemical)  recombination  at  low  currents  as  shown  through  the 
following  steps  often  regarded  as  the  Volmer-Tafel  mechanism  of  the  hydrogen  evolution  reaction 
mechanism  [72].  At  high  current  densities,  reaction  (2)  may  be  too  slow  [72]  so  the  reaction  mechanism 
may  be  coupled.  Here,  Ni-Cu  is  taken  to  represent  the  Ni-Cu  fee  lattice  containing  the  elements  listed  in 
Table  1  and  is  abbreviated  as  the  solvent  Ni. 

1  H20  +  Ni-Cu  +e*  =  Ni-Cu-Hads  +OH‘  Eq.  12 

2  Ni-Hads  +  Ni-Hads  =  H2(g)  +  2Ni  Eq.  1 3 

Of  course,  hydrogen  absorption  and  uptake  can  occur  via  a  third  step.  This  third  step  is  regarding 
to  be  fast  relative  to  slow  diffusion  in  the  Ni-Cu  fee  lattice  and  as  such  is  regarded  to  be  at  equilibrium: 

3  Ni-Hads  =  Ni-Habs  Eq.  14 

Hence,  the  relationship  between  cathodic  potential  and  hydrogen  surface  coverage  as  well  as 
absorbed  hydrogen  concentration  and  potential  are  given  as: 

1  / — /?riF\ 

h  =  kt  F(  1  -  0Ni_H)exp  (-p-J 

2  i2  =  k2  F(QNi_H)2 

The  Tafel  slope  of  (1)  is  rate  controlling  and  P=0.5  exceeds  1 18  mV/dec  but  is  non-linear  with 
coverage.  Figure  23  is  consistent  with  this  prediction.  At  low  coverage  and  for  coupled  rates  or  similar 
rate  determining  reaction  rates  for  1  and  2,  it  may  be  shown  that  the  coverage  of  hydrogen  Ni-Hads  is 
given  by  the  following  reaction  and  the  uptake  of  hydrogen  may  be  combined  to  yield  the  lattice 
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Eq.  15 

Eq.  16 


hydrogen  content  [74]: 


Eq.  17 


r  _  ^abs  n 

lH, lattice  —  T  uNi-H 
Kdes 


Eq.  18 


f*  _  ^  a  j  _  ^abs 

^H, lattice  “  &H -Lattice W lattice  ~  7 

Kdes 


Eq.  19 


where,  kabs,  kdes,  ki,  and  k2are  all  rate  constants,  q  is  the  hydrogen  overpotential  while  R,T,  and  F  have 
their  usual  meanings.  From  this  equation  it  can  be  seen  that  CHjattice  is  a  strong  function  of  hydrogen 
overpotential,  and  four  rate  constants.  One  or  more  of  these  rate  constants  are  most  likely  altered  by  heat 
to  heat  variations.  CH,  lattice  is  related  to  CH,diff taking  into  consideration  weak  energy  reversible  traps  such 
as  dislocations,  copper  solute,  twin  boundaries,  low  energy  grain  boundaries  and  some  weak  sites  on  y* 
particles.  Finally,  CH)tot  represents  the  sum  of  the  lattice  hydrogen  and  all  traps  (reversible,  irreversible  and 
otherwise)  as  discussed  extensively  elsewhere  [43]. 

These  factors  can  all  produce  further  differences  due  to  heat-to-heat  variations  in  the  parameters 
effecting  these  equations. 


Eq.  20 


Eq.  21 


1  °  Hrevtrap 


Eq.  22 


where  0  is  the  fractional  trap  site  coverage  and  N  is  the  trap  site  density,  Gutticc  is  given  by  Cl/Nl. 
expressed  as  atoms/interstitial  site  or  as  a  fractional  coverage.  In  these  equations,  it  can  be  seen  that  the 
combination  of  reversible  trap  site  density  and  coverage  dictated  by  binding  energies,  EB,  of  each  specific 
“weak”  trap  site  where  hydrogen  is  mobile  and  in  equilibrium  with  lattice  sites  determines  CH,revtraP  from 
each  weak  trap.  The  y’  binding  energy  was  calculated  from  the  subtraction  of  Edtj  from  Ed>r.  The  average 
values  for  the  binding  energies  for  y‘  ranged  from  7.8  to  19.1  kJ/mol  in  the  order:  NRL  HS  <  TR  2  < 
Allvac  <  NRL  LS.  The  -8  kJ/mole  value  may  be  regarded  as  reversibly  trapped  hydrogen  [75]  at  y4 
precipitates  while  the  trapping  near  -19  kJ/mole  may  be  irreversible  at  short  times.  This  weakly  trapped 
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hydrogen  in  combination  with  CH,jattice  determine  CH,diff-  The  status  of  grain  boundaries  as  traps  is  unclear 
since  when  EB  is  high  such  as  might  be  seen  on  high  angle  boundaries;  grain  boundary  traps  may  be 
irreversible  but  likely  represent  a  combination  of  reversible  and  irreversible  sites  depending  on  boundary 
angle,  segregation  and  degree  of  misfit.  It  should  be  noted  that  previous  work  detected  no  difference  in 
the  fraction  of  low  angle  boundaries  amongst  the  heats  examined  [36]  suggesting  that  grain  boundary 
engineering  considerations  may  not  explain  susceptibility  differences  [13]. 

Given  a  likely  P  value  of  0.5,  then  the  CH  lattice  CH,diff,  and  CH,tot  of  hydrogen  may  be  predicted  as 
a  function  of  hydrogen  overpotential  r\  as  seen  in  Figure  24.  Based  on  previous  comparison  of  lattice  and 
diffusible  hydrogen  determinations  [43],  the  lattice  and  diffusible  hydrogen  concentration  are  similar 
because  of  the  high  lattice  solubility  for  hydrogen  in  fee  nickel  copper  and  given  by  the  lower  curve  of 
Figure  24  (and  left  hand  peaks  in  Figures  10-14)  assuming  low  energy  traps  are  in  equilibrium  with  lattice 
interstitial  sites.  It  can  be  seen  that  Cn.diff  and  Ch, Total  are  indeed  exponentially  related  to  hydrogen 
overpotential  (Figure  24).  In  the  case  of  Allvac,  it  is  likely  based  on  theoretical  grounds  established  above 
that  the  pre-exponential  factor  in  Equation  19  is  lower  while  it  is  higher  for  TR2  for  reasons  presently 
uncertain. 

Moreover,  the  hydrogen  coverage  at  grain  boundaries  0gb  in  the  fracture  process  zone  for  a  heat 
invariant  grain  boundary  trap  binding  energy  would  then  be  decreased  owing  to  the  lower  CH, lattice  present 
in  Allvac  (Figure  24)  and  this  is  shown  in  Figure  36  where  for  a  given  grain  boundary  trap  binding 
energy,  0gbis  lower  when  CH,iattice  and  equivalently  CH,diffare  lowered.  This  would  be  the  case  for  Allvac  as 
an  example.  Also  shown  in  Figure  36  is  the  effect  of  grain  boundary  trap  binding  energy  on  0gb.  Figure  3 1 
was  obtained  by  resolutionizing  to  eliminate  y’  traps.  Recall  that  a  value  of  approximately  24.3  kJ/mole 
was  determined  in  resolutionized  TR2  Monel  K-500  from  Figure  31.  A  reasonable  theoretical  range  from 
10-30  kJ/mole  is  considered  in  Figure  36. 15  A  higher  grain  boundary  trap  binding  energy  results  in  a 
higher  coverage.  Previous  results  by  EBSD  show  a  similar  fraction  of  low  and  high  angle  grain 
boundaries  across  all  heats  [36].  Therefore,  an  argument  cannot  be  made  that  EB  is  greater  in  any 
particular  heat  based  on  differences  in  grain  misorientations.  The  potential  for  low  S  segregation,  high 
sulfur  gettering  element  content,  and  high  grain  boundary  surface  to  volume  ratio  could  indeed  lower  the 
grain  boundary  trapping  and/or  hydrogen  activity  in  the  Allvac  heat  rendering  it  more  difficult  to  reach 

15  It  should  be  noted  that  the  possibility  exists  that  any  given  high  angle 
boundary  traps  hydrogen  with  a  greater  trap  binding  energy  that  is  undetected 
by  TDS  of  polycrystalline  material.  This  is  because  any  single  boundary  type 
represents  too  low  of  a  number  of  trap  sites  for  detection  by  TDS  and  only  a 
large  number  of  traps  sites  with  a  binding  energy  above  the  lattice 
desorption  energy  can  be  detected. 
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high  0gb  in  Allvac  and  comparatively  easier  in  TR2  and  NRL  LS.  The  sulfur  issue  is  complex  and  should 
be  discussed  in  further  detail  below. 

Grain  boundary  Segregation  of  Sulfur  and  Hydrogen  -  Trace  Elements  Effects 

Sulfur  is  observed  to  segregate  to  grain  boundaries  in  nickel  and  nickel  based  alloys  [1,  8,  9,  14, 
23,  24,  32,  68,  69,  76,  77].  Auger  compositional  analysis  is  difficult  on  grain  boundaries  owing  to  carbon 
contamination,  differences  in  grain  boundary  inclination  relative  to  incident  electron  beam  and  auger 
electron  take  off  angle,  machine  specific  sensitivity  factors  for  elements,  and  many  other  factors  [68,  69]. 
However,  the  equilibrium  sulfur  coverage  at  boundaries  may  be  estimated  in  Ni  as  a  function  of  the  sulfur 
solubility’  [22,  68].  Here  we  estimate  the  sulfur  grain  boundary  atomic  fraction  based  on  McLean  theory 
[78].  Use  of  this  theory  is  supported  by  Auger  data  which  shows  narrow  surface  enhancements  of  S  which 
fall  rapidly  with  sputter  depth.  The  sulfur  solubility  can  be  estimated  for  equilibrium  sulfur  segregation 
with  an  activation  energy  Es  of  72.8  kJ/mole  and  at  the  aging  temperature  of  866  K  [68]. 

=  CLikexp  (|f)  Eq’23 

gt>  1 4-  rs  exv(ME\ 

f  +  LbulkexP \RTJ. 

where,  CbuIk  is  the  bulk  equilibrium  sulfur  concentration  expressed  as  an  atomic  fraction  and  Csgb  is  the 
sulfur  content  of  grain  boundaries.  The  bulk  sulfur  concentration  in  Ni  and  similarly  Ni-lOCu  are  both 
about  30  x  10'°  at  fraction  at  866  K  which  is  lower  than  the  sulfur  contents  of  all  alloys  tested  (Table  1) 
[68].  Use  of  the  equilibrium  concentration  at  866  K,  e.g.,  the  aging  temperature,  yields  an  equilibrium 
sulfur  segregation  atomic  fraction  of  0.42. 16  This  implies  that  every  heat  of  Monel  K-500  would  be 
expected  to  exhibit  equal  sulfur  segregation  not  accounting  for  crystallographic  effects  on  Es.  However, 
Hf  concentration  as  well  as  A1  in  Ni  lower  the  S  bulk  solubility  such  that  the  bulk  content  may  only  be  30 
x  10'6  at  fraction  at  866  K  [68].  This  produces  a  dramatic  decrease  in  Csgbto  only  0.069  atomic  fraction.  It 
is  interesting  to  note  that  this  estimation  is  in  some  agreement  with  experimental  data.  Thus  the  joint 
effect  a  of  high  gettering  elements  that  also  effect  S  solubility  such  as  Hf  and  to  a  lesser  extent  Zr,  Ca,  and 
Mg,  may  exert  a  substantial  effect  on  equilibrium  segregation  tendency  explaining  the  great  resistance  of 
Allvac  and  NRL  HS  Monel  K-500. 


16  Mujlford  reports  that  equilibrium  segregation  at  the  aging  temperatures  for 
Monel  K-500  in  <40  hours.  It  should  also  be  noted  that  peak  equilibrium 
segregation  of  sulfur  is  both  predicted  and  measured  in  the  aging  range 
[68]  R . A.  Mulford,  Grain-Boundary  Segregation  in  Ni  and  Binary  Ni-Alloys 
Doped  with  Sulfur,  Metall  Trans  A,  14  (1983)  865-870.. 
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Prediction  of  Breaking  Stress  and  Influence  of  Parameters  discussed  above 


The  breaking  stress  can  be  calculated  as  a  function  of  the  hydrogen  content  and  based  on  the 
lower  hydrogen  levels  in  Allvac  and  possibility  lower  segregated  sulfur  levels  predict  the  trend  seen  in 
Figures  2,  3  and  4.  An  unresolved  issue  concerns  whether  or  not  co-segregation  of  atomic  hydrogen  is 
required  [1],  and  whether  co-segregation  of  hydrogen  and  sulfur  act  in  an  additive  manner  [32,  79-81]  or 
synergistically  [1,  8,  82]  to  embrittle  boundaries.  The  atomistic  process  by  which  interface  strength  is 
reduced  by  a  foreign  atom  is  debated.  The  decohesion  model  proposes  that  a  sufficient  enrichment  of 
hydrogen  or  sulfur  at  boundaries  causes  a  discernible  weakening  of  bonds  between  adjacent  atoms  [83].  If 
hydrogen  and  sulfur  accumulate  at  a  planar  defect,  then  the  decohesive  strength  is  lowered  [26] 
selectively  along  that  interface.  This  can  lead  to  preferential  breakage  of  bonds  given  sufficient  applied 
stress  [84,  85].  Embedded  atom  calculations  show  that  a  £9  (221)  tilt  boundary  in  nickel  experiences  a 
decrease  in  cohesive  strength  from  18  to  8  GPa  when  the  segregated  monolayer  coverage  increases  from  1 
x  10’17  to  6  x  10’16  H  atoms/cm2  [85].  Various  expressions  have  been  used  to  describe  how  grain  boundary 
fracture  stresses  are  each  reduced  by  such  impurity  segregation.  The  general  form  of  these  expressions,  as 
first  described  by  Akhurst  [38]  and  Briant  [16]  gives  a  linear  relationship  between  fracture  stress  and 
hydrogen,  Enos  reported  a  logarithmic  relationship  [86,  87].  Herein  the  following  relationship  describes 
the  behavior  seen  in  Figure  2  using  the  exponential  uptake  law  given  in  Eq.  19  and  the  fracture  criterion 
expressed  below. 


cr fracture 


=  g0  —  aHlog 
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Eq.  24 
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where  aH  is  a  potency  factor  for  hydrogen,  fis  is  a  potency  factor  for  sulfur  along  the  lines  of  elsewhere 
[15,  36],  g0  is  the  notch  breaking  stress  in  air  or  at  minimal  CH  di^0  and  CsuifUr0  below  the  levels  for 
intergranular  cracking.  Gfracture  is  the  fracture  stress  which  causes  intergranular  cracking  when  equaled 
by  the  applied  stress.  This  expression  may  be  modified  for  the  stress  field  of  the  notch  as  follows;  the 
Cn.diff  could  be  modified  by  the  tensile  triaxial  stress  field,  the  sulfur  content  could  be  modified  by  the 
boundary  content,  and  the  applied  stress  could  be  given  by  the  local  stress  concentration  instead  of  remote 
stress.  The  local  applied  stress  can  be  raised  to  Gjracture  as  a  function  of  the  geometric  stress 
concentration  (e.g.,  notch  or  crack  tip  radius),  plastic  flow  characteristics  of  the  material,  as  well  as  global 
and  local  plastic  constraint.  Equation  24  does  not  predict  a  critical  impurity  or  hydrogen  concentration 
threshold,  yet  one  may  still  be  observed.  This  is  because  it  is  reasonable  to  assume  that  some  grain 
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boundary  hydrogen  concentration  will  decrease  <Jfracture  below  some  other  critical  stress  or  strain 
associated  with  the  criteria  for  intergranular  fracture,  otherwise  a  transition  from  ductile  to  intergranular 
cracking  seen  in  Figure  3  is  not  seen.  Another  issue  is  whether  the  mean  or  the  extreme  value  of  local 
concentrations  determined  from  the  distribution  of  measured  values  should  be  used.  The  distribution  is 
expected  due  to  various  crystallographic  grain  boundary  mis-orientations  that  alter  segregation  tendencies 
on  certain  grain  boundaries 

Moreover,  all  of  the  factors  discussed  above  reflected  in  heat-to-heat  variations  are  either 
incorporated  in  Eq.  24-25,  including  hydrogen  sulfur  segregation  and  yield  strength,  or  captured  above  in 
the  relationships  for  sulfur  segregation  and  hydrogen  absorption  that  are  inputs  to  this  equation.  As 
absorbed  hydrogen  and  equilibrium  sulfur  content  increase  then  the  H  and  S  trapped  on  boundaries 
increases,  subsequently  (J/racture  decreases.  Yield  strength  increases  the  local  value  of  the  applied  stress, 
<7 applied  at  any  given  applied  remote  stress.  It  can  be  seen  by  exercising  this  relationship  that  Allvac  and 
NRL  HS  are  more  resistant  while  TR2  and  NRL  LS  are  less  resistant.  It  can  also  be  seen  that  yield 
strength  is  not  necessarily  a  dominant  factor  in  this  case  although  usually  an  important  factor  in  hydrogen 
embrittlement. 


Conclusions 

Monel  K-500  (63%Ni,  33  %Cu,  3.15  %AI,  wt.  %)  was  susceptible  to  hydrogen  embrittlement 
under  severe  cathodic  polarization  (e.g.,  -1.1  VSce)  in  seawater  in  all  heats  of  aged  material  that  were 
within  a  federal  specification.  Differences  between  field  susceptible  heats  and  others  explored  and 
compared.  Monel  K-500  harvested  from  field  failures  cracked  intergranularly  more  readily  at  modest 
levels  of  cathodic  applied  potentials.  Moreover,  the  extent  of  intergranular  fracture  increased  with  more 
negative  applied  potentials.  More  facile  intergranular  fracture  correlated  with  sulfur  segregation  and 
extent  of  applied  cathodic  potential.  H  absorption  increased  strongly  with  hydrogen  overpotential  and  was 
not  a  simple  function  of  hardness  for  material  within  the  specification.  Various  assessments  of  hydrogen 
concentration,  by  three  investigators  all  show  that  Allvac  (high  HEAC  resistance  at  -0.850  and  -.950 
Vsce)  absorbs  less  hydrogen  than  harvested  materials  from  field  failures  and  does  not  crack 
intergranularly  at  -.950  V  vs.  SCE..  Moreover,  grain  boundary  chemistry  was  assessed  after  aging  and 
hydrogen  fracture  by  Auger  Electron  Spectroscopy  .  Some  degree  of  sulfur  enrichment  was  observed  in  all 
heats.  The  heats  observed  to  crack  in  seawater  and  field  testing  contained  unfavorable  concentrations  of 
trace  detrimental  grain  boundary  segregating  elements  relative  to  beneficial  sulfur  gettering  elements  in 

combination  with  a  low  density  of  grain  boundary  trap  or  segregation  sites  per  unit  volume.  Other 

84 


hydrogen  metal  interactions  such  as  trapping  and  diffusion  rate  could  were  not  key  metrics  indicting 
susceptibility.  Greater  susceptibility  is  traced  to  enhanced  hydrogen  uptake  and  speculated  grain 
boundary  sulfur  enrichment  based  on  a  balance  of  alloy  trace  chemistries  coupled  with  grain  boundary 
trapping  of  hydrogen  supplied  from  the  lattice  and  possibly  lower  energy  y’  trap  sites  which  absorbed 
hydrogen.  Material  hardness  within  the  limits  from  Rockwell  C  28  to  36  were  not  found  to  be  the 
dominant  factors  in  HE  susceptibility. 
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Table  1.  Chemical  composition  in  wt%  or  ppm  for  Monel  K-500  heats  (UNS  N05500)  measured  in  accordance  with  ASTM  E2594  (ICP-OES) 
[88],  unless  otherwise  indicated. 
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Table  2.  Hardness  and  grain  size  values  for  Monel  K-500  heats  (UNS  N05500)  measured  in  accordance 
with  ASTM  E2594  (ICP-OES)  [88]  and  determined  from  EBSD  data  and  published  previously  [36]. 


Heats 

HRC 

Grain  size  (pm) 

All  vac 

31 

13.8 

TR  2 

34.1  ±0.6 

35.3 

NRL-LS 

28.50  ±  1.63 

22.5 

NRL-HS 

36.65  ±0.27 

11.2 

Table  3.  S,  normalized  S,  S/Ni  and  normalized  S/Ni  in  at  %  or  at%/at%  on  Monel  K-500  obtained  via 
AES  analysis  on  the  grain  boundaries  of  the  intergranular  fracture  surfaces  resulting  from  SSRT  at  -1.1 
Vsce  in  deaerated  0.6  M  NaCl  (pH:  8.0  by  adjustment  with  NaOH). 


at%  or 
at%/at% 

TR2 

A1 

TR2 

A3 

Allvac 

A1 

Allvac 

A2 

NRLHS 

A2 

NRLHS 

A3 

NRLLS 

A1 

NRLLS 

A2 

NRL¬ 
LS  A3 

Highest  S 
(at%) 

Depth  (nm) 

2.7 

1.5  nm 

1.7 

1.5  nm 

3.3 

2  nm 

4.1 

1  nm 

3.1 

1  nm 

3.1 

1  nm 

1.7 

1.85 

nm 

3.8 

0.93 

nm 

1.9 

1.85 

nm 

S 

Normalized 
at  2 1  nm 

3.3 

4.0 

3.5 

11.4 

2.2 

2.4 

2.6 

14.6 

3.1 

S/Ni  (max 
value  and 
location) 

0.1 

(0.1,0 

nm) 

0.0 

(0.1,0 

nm) 

0.2 

(0.2,  1 
nm) 

0.1 

0.2 

0.1 

0.1 

(0.1, 

0.9  nm) 

0.11 

0.1 

(0.1, 

0.9 

nm) 

S/Ni 

Normalized 
at  21  nm 

5.7 

(7.7,  0 
nm) 

4.1 

(14.3,0 

nm) 

11.4 
(12.6,  1 
nm) 

24.6 

6.3 

6.7 

6.71 

(7.6, 

0.9  nm) 

28.97 

4.9 

(6.1, 

0.9 

nm) 

Table  4.  Measured  hydrogen  desorption  and  binding  energies  for  the  aged  Monel  K-500  heats. 


Ed  (kJ/mol) 

SD  (kJ/mol) 

Eb  (kJ/mol) 

left 

right 

left 

right 

Allvac 

23.0 

40.6 

1.4 

0.5 

17.7 

TR  2 

26.4 

38.8 

0.7 

0.7 

12.4 

NRL  LS 

26.0 

45.0 

0.2 

0.4 

19.1 

NRLHS 

33.6 

41.3 

1.6 

0.9 

7.7 
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Figure  18.  Time  to  failure  of  TR2,  Allvac,  NRL-LS  and  NRL-HS  heats  of  Monel  K-500  during  concurrent  charging/SSRT  test 
after  electrochemical  H  pre-charging  in  deaerated  0.6  M  NaCl  during  48  hours  at  a  constant  potential  of -0.850,  -0.950,  -1.0  and  - 

1-1  vSCE. 
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Figure  19.  Fractography  ofTR  2  and  Allvac  heals  cf  Monel  K-500  after  concurrent  charging/SSRT  test  running  the  experiment  a) 
in  a  r  and  prechargmg/charging  at  a  constant  po:entiaI  of  b)  -0.850,  c)  -0.950  and  d)  -1.1  Vsce  ir.  deaerated  0  6  M  NaCl  (pH:  8.0 
by  adjustment  with  NaOH). 


Figure  20.  Fractography  ofNRL  LS  and  HS  lots  cf  Monel  K-500  after  concurrent  charging/SSRT  test  at  -1.1  VSce  in  deaera:ed 
0.6  M  NaCl  (pH:  8  0  by  adjustment  with  NaOH) 
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Figure  21.  Normalized  Auger  spectrum  on  Monel  K-500  with  respect  to  signal  at  a  depth  of  21  nm  from  a  grain  boundary  for  a) 
TR  2,  including  procedure  used  to  make  the  calculation,  b)  Allvac,  c)  NRL  LS  and  d)  NRL  HS  after  SSRT  performed  at  a  applied 
potential  of -1 .1  VSCE  in  deaerated  0  6  M  NaCl  (pH:  8.0  by  adjustment  with  NaOH). 


I  Surface  Area  of  GB  to  volume  ratio  (mm2/mm3)  (xIO) 
I  Cs  in  eq  (ppm) 

I  Cget  (ppm) 

I  S+Sn+$b-Hf-Zr/3  (ppm)  [/1 0] 


Figure  22.  Chemical  characterization  of  Monel  K-500  grain  boundaries  described  as  surface  area  of  grain  boundary  to  volume 
ratio;  and  trace  elements  content  pertaining  to  the  deleterious  or  enhanced  effect  of  these  elements  with  respect  to  HE 
susceptibility.  The  S+Sn+Sb-Hf-7,r/3  values  were  TR-2  (-97.0;  NRL  LS  -  76.7;  Allvac  -121.3;  NRL  HS  -197.9  when  expressed 
in  ppm  concentrations) 
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a) 


b) 

Figure  23.  Evaluation  of  the  cathodic  kinetics  of  hydrogen  evolution  on  TR  2  and  Allvac  Monel  K-500  materials  in  0.6  M  NaCl 

(pH  8  00  by  adjustment  with  NaOH). 
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Figure  24.  Diffusible  and  total  hydrogen  concentrations  determined  by  the  electrochemical  method,  barnacle  cell,  and  thermal 
desorption  spectroscopy,  respectively,  for  the  different  Monel  K-500  heats  after  charging  in  VSCe  in  deaerated  0.6  M  NaCl  (pH: 
8.0  by  adjustment  with  NaOFI)  at  -0.850,  -0.850,  -1.0  and  -1.1  VSce  during  7  days  (the  95%  confidence  intervals  are  shown  for 
the  mean  diffusible  hydrogen  concentration  values) 


Temp  (°C) 


Figure  25.  Normalized  hydrogen  pressure  egress  for  as-received  and  baked  TR  2  Monel  K-500  at  450  °C. 
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a) 


1000/Tm  (K'1) 


b) 

Figure  26.  Experimental  results  extracted  from  ramped  TDS  experiments  performed  at  different  heating  rates  on  120  pm  ±  10  pm 
thick  Allvac  heat  Monel  K-500  electrochemically  charged  with  hydrogen  in  0.6  M  NaCl  (pH:  8.0  by  adjustment  with  NaOH) 
during  7  days  applying  a  potential  of  -1.0  VSce-  a)  Hydrogen  desorption  rate  (dCH/dt)  as  a  fiinction  of  temperature  and  b) 
calculation  of  the  Ed  for  the  identified  hydrogen  desorption  peaks  and  their  95%confidence  intervals. 
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Figure  27. Experimental  results  extracted  from  ramped  TDS  experiments  performed  at  different  heating  rates  on  120  pm  ±  10  pm 
thick  TR  2  heat  Monel  K-500  electrochemical ly  charged  with  hydrogen  in  0.6  M  NaCl  (pH:  8.0  by  adjustment  with  NaOH) 
during  7  days  applying  a  potential  of  -1.0  Vsce-  a)  Hydrogen  desorption  rate  (dCn/dt)  as  a  function  of  temperature  and  b) 
calculation  of  the  Ed  for  the  identified  hydrogen  desorption  peaks  and  their  95%confidence  intervals. 
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1000/Tm  (K'1) 


b) 

Figure  28.  Experimental  results  extracted  from  ramped  TDS  experiments  performed  at  different  heating  rates  on  120  pm  ±  10  pm 
thick  NRL-LS  Monel  K-500  electrochemically  charged  with  hydrogen  in  0.6  M  NaCl  (pH:  8.0  by  adjustment  with  NaOH)  during 
7  days  applying  a  potential  of  -1  0  VSCe-  a)  Hydrogen  desorption  rate  (dCn/dt)  as  a  function  of  temperature  and  b)  calculation  of 
the  Ed  for  the  identified  hydrogen  desorption  peaks  and  their  95%confidence  intervals. 


102 


r 


a) 


1000/TJK'1) 


b) 

Figure  29.  Experimental  results  extracted  from  ramped  1 DS  experiments  performed  at  different  heating  rates  on  120  pm  ±  10  pm 
thick  NRL-HS  heat  Monel  K-500  electrochemically  charged  with  hydrogen  in  0.6  M  NaCl  (pH:  8.0  by  adjustment  with  NaOH) 
during  7  days  applying  a  potential  of  -1.0  VSce-  a)  Hydrogen  desorption  rate  (dCH/dt)  as  a  function  of  temperature  and  b) 
calculation  of  the  Ed  for  the  identified  hydrogen  desorption  peaks  and  their  95%confidence  intervals. 
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2,5x10 


Figure  30,  Hydrogen  desorption  rate  (dCH/dt)  with  respect  to  temperature  calculated  from  ramped  TDS  experiments  performed 
on  1 20  pm  ±  10  pm  thick  Monel  K-500  for  the  different  heats  at  a  heating  rate  of  3  °C/min.  The  flat  specimens  were 
electrochemically  charged  with  hydrogen  in  0.6  M  NaCl  (pH:  8  0  by  adjustment  with  NaOH)  during  7  days  at  -1 .1  Vsce- 


Temperature  (°C) 


Figure  31 .  Desorption  rate  versus  temperature  obtained  utilizing  thermal  desorption  spectroscopy  at  a  heating  rate  of  3  °C/min  for 
charged  solution  heat  treated  and  aged  Monel  K-500,  solution  heat  treated  and  water  quenched  material  and  solution  heat  treated 
and  water  quenched  hydrogen  desorbed  material.  For  this  last  material  hydrogen  desorption  took  place  at  25  °C  during  6  days 
The  specimens  were  electrochemically  charged  with  hydrogen  in  0.6  M  NaCl  (pH:  8.0  by  adjustment  with  NaOH)  during  7  days 
applying  a  potential  of -1.0  Vsce-  The  heat  used  for  these  experiments  was  the  TR  2  heat. 
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Figure  32.  Saturation  curve  analysis  represented  in  terms  of  the  H  mass  ratio  with  respect  to  a  hydrogen  saturated  state  versus  the 
suare  root  of  time  for  NRL  heats  of  Monel  K-500.  Samples  were  pre-charged  at  different  time  intervals  at  -0.800,  -0.850,  -0.890 
and  -1  000  Vsce  in  seawater.  Theoretical  calculations  using  the  solution  for  a  non-steady  diffusion  model  for  a  cylinder  with 
known  initial  hydrogen  profile  was  utilizd  to  calculate  the  saturation  curve  for  hydrogen  difTusivities  of  5  x  10‘10  cnr/s,  1  x  10'10 
cm2/s  and  2  x  10*10  cm2/s. 


Hardness  (HRC) 


Figure  33.  Comparison  between  difussible  hydrogen  concentration  of  the  different  heats  of  Monel  K-500  versus  hardness  of 
plane  sheet  specimens  electrochemically  charged  with  hydrogen  in  0.6  M  NaCl  (pH:  8.0  by  adjustment  with  NaOH)  during  7 
days  applying  a  potential  of -0.850  and  -1. 100  VSCE 
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Figure  34.  Effect  of  hardness  on  the  maximum  remote  stress  and  time  to  failure  of  the  different  heats  Monel  K-500  tested  during 
concurrent  charging/SSRT  test  after  electrochemical  H  pre-charging  in  deaerated  0.6  M  NaCl  (pH:  8.0  by  adjustment  with 
NaOH)  during  48  hours  at  a  constant  potential  of-1 .1  VSCE 
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Figure  35.  Effect  of  hardness  on  the  y'  binding  energy  of  the  different  heats  of  Monel  K-500. 
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Figure  36.  Effect  of  binding  energy  and  diffusible  hydrogen  concentration  on  the  theoretical  hydrogen  coverage  on  the  grain 
boundaries.  The  following  parameters  were  assumed  for  the  calculations*  pMond  k-5oo^.44  g/cm3,  NNj3A|:  8.29  x  1022  sites/  cm3 
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